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PREFACE 


NASA has an interest in turbines related priinaril}^ to aeronautics and 
space applications. Airbreathing turbine engimvs provide' j('t and turbo- 
shaft propulsion, as well as auxiliary power for aircraft. Prope'Ilant- 
driven turbines provide rocket propulsion and auxiliary pow(‘r for spac('- 
craft. Clos('d-cycle turbine engines using inert gases, organic fluids, and 
metal fluids have been studied for providing long-duration ('lectric powc'r 
for spacecraft. Other applications of current interest for turbine' engines 
include land-vehicle (cars, trucks, buses, trains, etc.) propulsion powcT 
and ground-based electrical pe)wer. 

In view of the turbine-system interest and efforts at Lewis Rc'search 
Center, a course entith'd “Turbine Design and Application’^ was ])re- 
sented during 1968-09 as part of the In-House Graduate' Study Program. 
The course was somewhat revised and again pres(‘iit('d in 1972-7‘L Various 
aspects of turbine technology were covered including thermodynamic and 
fluid-dynamic concepts, fundamental turbine concepts, vt'locity dia- 
grams, losses, blade aerodynamic design, blade cooling, mechanical de- 
sign, operation, and performance. 

The notes written and used for the course have been r<‘vised and edited 
for publication. Such a publication can serve as a foundation for an intro- 
ductory turbine course, a means for self-study, or a reh'n'iice for sc'lected 
topics. The first volume presented the material covering thermodynamic 
and fluid-dynamic concepts, fundamental turbine' conce'pts, and vc'locity 
diagram design. This second volume presents the material related to 
blade aerodynamic design and turbine' (‘iiergy losses. 

Any consistent set of units will satisfy the (*quations prese'iited. Two 
commonly used consistent sets of units and constant value's are given 
after the symbol definitions. These are the SI units and the U.S. cus- 
tomary units. A single set e>f equations covers be>th sets of units by in- 
cluding all cemstants required fe>r the U.S. customary units and defining 
as unity those not required for the SI units. 


Arthur J. Glassman 
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CHAPTER 4 


Blade Design 


By Warner L. Stewart and 
Arthur i. Glassman 


Tho (l(\sign of a tiirhino consists of throe major stops. Tho first is the 
d(^t('rminat ion of tho ov(‘rall roejuiromonts of flow, work, and spood. Those 
ar(‘ usually (‘stablishod l)y tho particular application. The second stop is 
tho (‘volution of velocity diagrams consistent with the desired efficiency 
and/or number of stages. This was discussed in chapter 3. The third step 
is th(‘ d(‘sign of tlu‘ blading that will produce the flow angles and velocities 
ro(iuir(‘d by tho velocity diagrams. This stop involves the determination of 
tho size, shape, and spacing of the blades. 

This chapt(‘r covers som(‘ of tho more important aspects of blade design. 
Th(‘ height of tho blad(‘ is s(‘t by the overall requirements of flow, speed, 
and inl(‘t stat<‘ conditions and the selected velocity diagram, which dictates 
the fluid state conditions throughout the turbine. The blade chord is 
usually s(‘l(‘(‘t(‘d to b(‘ a minimum valu<‘ consistent with mechanical con- 
sid(‘rations. The chord must b(‘ long c‘iiough to allow accurate fabrication 
and assur(‘ structural int('grity during operation. The selection of blade 
spacing, which can be ex])ress(‘d nondimensionally as solidity (ratio of 
chord to spacing) or axial solidity (ratio of axial chord to spacing), in- 
volv(‘s many coiisid(‘rations that will be discussed in the first part of this 
chapter. lflad<‘ profile design, which includes blade exit and inlet geome- 
tri(\s as well as th<‘ connecting surface profiles, is then discussed in the last 
part of this chapter, (,'hannel flow theory, which is the basis for the 
analytical procedures us(‘d to accomplish the jirofile design, is discussed 
in the next chapt<‘r. 
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SOLIDITY 

One of the important aspects of turbine blading design is the selection 
of the blade solidity, which is the ratio of chord or axial (‘hord to spacing. 
A minimum vuliu^ is usually desired from the stand[>oint of reducing 
weight, cooling flow, and cost. However, chord r(‘duction is limit(‘d by 
mechanical considerations, and incr(‘as(‘d spacing (‘V(*ntually results in 
decreased blade' ('fficic'iicy due to separate'd flow. This s(*ction will concern 
its(‘lf with th(^ aerodynamic factors affe'cting solidity seh'ction. Th(‘ dis- 
cussion will include' tlu' e'tTe'ct oi ve'locity diagram re'eiuiri'ments e>n solidity 
and the're'lation Ix'twe'e'n blade loading and solidity. Also included will be a 
description of advance'd blading concepts that are' be'ing studie'd for use* to 
suppre'ss se'paration and there'by re'duce' the' pe'rmissibie seelidity. 


Effect of Velocity Diagrams on Solidity 

Figure' 4-1 shows a typical s(*t ol blade' inle't and e'xit diagrams as well 
as the stati('-])re'ssure‘ elistribiitie>n aroimd a blade'. The' ve'lex-itie's in this 
figure are' shown as abse)]ute ve'locitie's. The* discussion in this chapte'r 
pertains to reel or blade' rows as we'll as to state)r })]ade* rows. Whe'ii re'fe'rring 
to a rotor, we' must use' re'lative rathe'r than absolute' velocitie'S in the 
('(juations and figure's. Since' in tliis chapter we' are' conce'rne'd with blade 
re)ws rathe'r than with stage's, the' angle' conve'ntioei will diffe'r slightly from 
that used in pre'vuous e'hapte'rs. 44ie‘ e'xit tange'ntial-ve'lex'ity compone'iit 
and flow angle are take'u as nc'gative' value's. The' inle't value's are' positive 
if the inle't and exit tatige'iitial-ve'le)city ceimpoiK'uts are in e>ppe)site dire'c- 
tions, and ne'gative if in the same direction. 

If e)iie conside'rs the' two-dime'nsie>naI flow through a passage of unit 
hc'ight betwe'en twe) blades, thc'u the tangential force' e'xe'rted by the 
fluid as it fle)ws fremi blade inle't (subscript 1) to (‘xit (subscript 2) is 


where 


»Sp2l X.2 ( T n ,i ~ T,/ , 2 ) 
9 


(4-1) 


F,t tange'utial force, X ; lb 

g conversion ce)ristant, 1 ; 32.17 ( Ibm ) (ft )/ (Ibf ) (se'c-) 

S' l)lade spae'ing, m; ft 

p density, kg/ni'^; Ib/ft* 

Ty axial ce)mpoiu'ut of velocity, m/se'c ; ft/ sec 
Vu tangential ce)mpe>nent of ve'lex'it v, m/sec ; ft se'c 

This tange'utial force' e'xe'rte'd by tlie' fluid must be' the* same* as the* force* 
due to the' static-pre*ssure distribution around the blade*, as was discusse'd 
in chapter 2. The* lowe'r part e>f figure 4-1 shows a typical static-pre'ssure* 
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Stations 1 2 



Fkjt'ui-: 4-1. Tyj)ical bhido-row velocity diagrams and surfa(*e static-pressure 

distribution. 


distribution around tlu‘ blade* row as a function axial distance*. The are*a 
be*twe*e*n tlie* twe> cur\'e*s re*prcse*nts the' teetal blade* fe>rce* acting e>n the 
fle>w in the* tange*ntial elire*ctie>n. Thus, 


F, = c. 




(4-2) 


whe*re* 

Cjt axial che>rd, m ; ft 

pp pr(*ssure'-surfa(‘e* static pre*ssure*, X/m-; U)/ft^ 
p„ suctiem-surface* static jire'ssure*, X/m’; lb ft^ 
X axial distance*, in ; ft 

The axial se)lidity, is 


Cx 


(4-8) 


3 



TURBIN K DKSIUN AND APPLICATION 


Substituting equations (4-1) and (4-2) into eciuation (4-3) then yields 


P2F.,2(F„.i-F..,) 



(4-4) 


At this point, we introduce two tangc^ntial loading (*()(*fiicients that have 
been used to relate the actual blade loading to an ideal l>iad(‘ loading. The 
first is the widely used cot^ffieient introduced by Zwcnfr^l (red. 1). This 
co(‘fficient is basc^d on an ideal loading that assunu's { 1 ) th(‘ static pn'ssun^ 
on the pn'ssurc' surface to b(‘ constant and etpial to th(^ inl(4 total pressun* 
and (2) th(* static pressure' on the suction surface* to luM-oustant and eepml 
to the exit static pressure. In eepiation form, 


where 



(4-5) 


xpz Zweif(*l loading coefficient 

]>i inlet total pressure, X/m-; Ib/ft- 

pz exit static pressure, X/in-; Ib/ft- 

The S(?cond coefficient is similarly ch'fined (*xc('pt that the' assuiiK'd 
constant static pre*ssure on the* suctie>n surface* is e*epial te> the* minimum 
value of static pressure (see fig. 4-1) on that surfae^e*. 4'hi.'^ h>ading ce>effi- 
cient can never exceed a value' of 1, and for all prae‘tie‘al pur])oses, it must 
always be less than 1. The Zwe'ifel co(4hcie*nt, e)ii the* e)the*r hand, can e'x- 
ceed a value of 1 . In eeiuation fe)rm, this se*(‘e)nd le>ading coe'fficie'iit }p is de*- 
fined as 


where p,,„nn is 
X'/m‘ e>r Ib/ft-. 

The velocity 
pressed as 



, — - (4-G) 

pi 

the minimum static jire'ssure on the suction surfae*e* in 
components in terms of ve'Iocity and flow angle are ex- 


and 


Vu ~ F sin a 

(4-7) 

Vx = V cos a 

(4-8) 
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when* 

r fluid v('l()cit y, in/s('c; ft/soc 

a fluid flow an^lc', d(‘g 

Substituting (‘([imtions (4-r)) or (4-()), (4-7), and (4-8) into (equation 
(4 4) and usin^ the tri^ononu'tric n'lation sin 2a = 2 sin a cos a yields 





(ifnvA 

1 ^ I (7v — 1 ) sin 2«2 


(4-9) 


K is dll' ratio of taiig('iitial velocity component ( V„,i) at the blade 
inl(‘t to that ( I n, 2) at th(^ blad(' (^xit. 

I)(rmit}o}i of i}ico)}iprei4Hihle~jlow Ib'Iations involving solidity, 

v(‘locity diagrams, and loading arc usually (‘volvc'd by assuming incom- 
pr('ssibl(‘ flow with no loss. With this assumption, dtaisity p is constant, 
and B(‘rnouIIi’s (‘(juation 

pr- (4-10) 

can b(‘ us(‘d. Substituting (‘(juation {4-10) into (Hpiation (4-9) yields 


( K—\ ) sin 20-2 ( A' — 1 ) sin 2^2 



(4-11) 


wlu‘r(‘ Vn.oT is the v(‘U)city on the* suction surface where p = ps,min- 
Let us now d(4ine a suction-surface diffusion param(‘ter as 





(4-12) 


Afany ])aram(4(Ts of this typ(' hav(‘ IxMai us(*d to represent a measure of 
th(‘ d(‘(*elerat ion of thc‘ flow on th(‘ suction surfac(\ This deceleration is an 
indication of tiu' sus(*('])tibility of th(' flow on the blade* to se[)arate. Using 
this d(‘finition ((*(j. (4-12)) in eeiuation (4-11) yields 


( A' — 1 ) sin 2 q (2 ( A' — 1 ) sin 2«2 


(4-13) 


hkjuation (4-18) shows that the* solidity param(4(*r or is 

constant for ('ar‘h particular v(‘locity-diagram reHiulreanent. Since loading 
coefh(*ient xp, whicli cannot (‘xc(*<‘d a valiK* of 1 , do(‘S not vary greatly, it 
can be seeai that d('cn*asing solidity n^sults primarily in increased suction- 
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surface diffusion (higher £)«), the conse(iu(‘nc(‘ of which will be discussed 
later in this chapt(T. The solidity param(‘t(T is |)lotted against th(' tan- 
gential velocity ratio K for several values of (^xit flow angles in figure 
4-2(a). A valiK' of /v^O repn^scaits a redaction blad(^ with axial inl(‘t, a 
value of K= — \ represents an inipuls(‘ bladc^, and a value of /v< — 1 
repr('S(‘uts a n(‘gativ(‘ reaction blade. Positive valiu's of K r(‘pr<‘S(‘iit inlet 
and exit tangcaitial velocities in th(‘ sarn(‘ dire^ction and arc‘ (‘iicountered 
primarily in the tip se'ctions of rotor blad(‘S. As seaai from (‘(juation (4-13) , 
solidity parame'teT is (‘qua! to zero for all (‘xit angle's fen* 1. This n']ire- 
sents the' case whe^re' theTe is ne> turning of the' flow. The' sejlidity ])aram(‘ter 
increase's with de'cn'asing K values. Thus, if e'xcf'ssive' suctie)n-surface 
diffusie)ii is to hr avoide'd, solidity must increase' as the' ve'locity diagrams 
move' from re'ae*tie)n toward impulse. It can be' se'e'ii that fe>r any give'ii 
value e)f A", a maximum value' of solidity parame'ter is e)btained with an 
e'xit angle' e>f 45"^. 

Ee juation (4- 13) can be modifie'd te) a function of the' ink't and C'xit 
angle's to yie'ld the- {*e[uatie)n de'rive'd in re'fe're'iice' 1. 

2 (‘OS «2 , 

cTj. = — sm (ai — a^) (4-14) 

\l/s cos 

For bre'vity, this is ('xpressed e>nly in terms of the coe^fficie'iit lujuation 
(4-14) shows that the' solidity parame'te'r ajr\p^ can be' e'xpre'sse'd in te'rms 
of the flow' angle's only. Solidity parame'te'r is })le)tte‘d against e'xit flej\v 
angle' for se've'ral value's of the' ird(*t flo\v angle' in figure' 4-2 (b). For a 
give'll e'xit angle', solidity parame'te'r increase's with incre'asing inle't angle*. 
In the* re'gion of most inte're'st («i>(P, « 2 < —45°), solidity paramete'r for 
each inlet angle de'cre'ase's w ith de'cre'asing exit angle. 

A third re'lation can be* e'volve'd, this one in te'rms of blade re'action A, 
which was de'fined in chapte'r 2 as 

13 ' 

= (4-15) 

1 2 “ 


Substituting ('([nation (4-8) inte> e'ejuation (4-15) yie'lds 


(4-16) 


for the* two-dime'iisional, incomjire'ssible'-flow case*, where lT,i=Fx, 2 . 
Substitution of ('ejuation (4-16) back into e'ejuation (4-14) then yie'lds 


9 

= - \Zl-Rsm Aa (4-17) 

iz 

whe're* Aa is ai — az. 
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Equation (4-17) oxprosses the solidity parameter in terms of blade 
reaction and turning angle. The solidity parameter is plotted against 
reactir)!! for s('V(ual values of turning angh^ in figure 4-2 (c) . It can be seen 
that, as indicated previously, the solidity parameter decreases with in- 
CH'asing redaction. The solidity parameter is a maximum for a turning 
angle of 90° and varies little with turning unk^ss very high or very low 
turnings are used. 

Radial variation . — Chapter 3 discussed the radial variations in velocity 
diagrams that must occur in order to satisfy both the varying blade speed 
and radial (Hiuilibrium, Sinc(^ axial solidity was shown to vary with 
var^dng velocity diagrams, there will be a radial variation in the desired 
value of axial solidity. The nature of this radial variation will be illustrated 
by an C'xample. Consider a single-stage turbine having axial inlet and exit 
flows (zero inlet and (*xit swirls), constant axial v(4o(*ities, a constant 
hub-to-tip-radius ratio of 0.7, an impulse rotor hub with a stator-hub exit 
flow angle of —70°, and free-vortex swirl distribution. For this case, the 
flow angles at the hub and tip and the corresponding solidity-parameter 
values comput('d from CHjuation (4-14) arc^ shown in the following table: 



Stator 

Rotor 


Inlet 

angle, 

(leg 

Exit 

angle, 

deg 

Solidity 

parameter, 

T nlet 
angle, 
deg 

Exit 

angle, 

4eg 

Solidity 

parameter, 

Hub 

0 

-70 

0.64 

54 

-54 

1.90 

Tip 

0 

-62 

.83 

-2 

-63 

.79 


Note again that the angle convention being used in this chapter is some- 
what diflferent from that of previous chapters. Herein, stator exit angles 
are negative. Assume that the loading coefficient is to be maintained 
constant radially. This is a reasonably desirable condition, and the 
assumption enables us to proportion solidity directly to the solidity 
parameter. 

Let us now determine how the hub and tip values of solidity parameter 
shown in the preceding table can be made physically consistent. The axial 
solidity variation in any blade row must be inversely proportional to 
radius (because blade spacing is directly proportional to radius) and 
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directly proportional to axial chord. For the stator, the axial solidity 
parameter at the hub is 0.64. If axial chord were held constant, then the 
corresponding tip value of the axial solidity parameter would be 
0.64X0.7 = 0.45, which is almost half of the desired value of 0.83. There- 
fore, a considerable axial taper from tip to hub is often used so that the 
axial chord can increase with radius and yield the higher solidities desired 
at the tip. 

In the case of the rotor, the axial solidity parameter at the hub is 1.90. 
If axial chord were held constant, then the corresponding tip value of the 
axial solidity parameter would be 1.90X0.7 = 1.33, which is still larger 
than the desired value of 0.79. Therefore, axial taper from hub to tip is 
often used in rotor blades so that axial chord can decrease with increasing 
radius and yield the lower solidities desired at the tip. Taper from hub to 
tip in the rotor is not only aerodynamically desirable, but is also mechani- 
cally desirable from the standpoint of reducing blade stress. To simplify 
fabrication in many cases, especially for smaller turbines, axial taper is not 
used, and there results a radial variation in loading coefficient. With the 
axial solidity selected on the basis of the mean-section velocity diagrams, 
this radial variation in loading coefficient in many cases, especially those 
where the blading is not highly loaded, will not have a severe effect on 
turbine performance. 

Effect of compressibility , — The term 




in equation (4-9) reduces to !//)« for incompressible flow conditions, as 
shown by equation (4-13). For a compressible flow case having the same 
loading coefficient ip as for incompressible flow, division of equation (4-9) 
by equation (4-13) yields 


Zq 

Ox ^ 


O'z.inc Vl Pi, min 


(4-18) 


where is the incompressible flow value as determined from an 

equation such as (4-13), (4-14), or (4-17). By introducing the relations 
between critical velocity ratio, density, and pressure (eqs. (1-3), (1-52), 
(1-61), (1-63), and (1-64)) and using the definition of D, (eq. (4-12)), 
equation (4-18) is modified to 
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CTx 

I ,inc. 


7 + 1 


.-1 /_F {I\ 

iVFcrAJ \Vj,^‘ 


y+ 


r 7-1 /rv 

1-' — ) As 

7+1 \l cr/2 


(4-19) 


whc're 

7 ratio of sp{‘cific heat at constant pressure to specific heat at 
constant volume 

Ver critical velocity, m/s('c; ft /see 

Then, by using binomial expansion and by neglecting the secondary terms, 
equation (4-19) can be approximated as 


= 1 - 


(fJ <fj 

' 7+1 2(7+iy 


(4-20) 


The approximation represented by c([uation (4-20) is quite good for 
( l//T/-r )2 valu(‘s up to about 1. The solidity ratio is plotted against 

suet ion-surf ac(‘ diffusion jiaranu'tcT for several values of critical velocity 
ratio in figur(‘ 4-3. The compressibility effect becomes more pronounced 
as D, (‘ith(T incr(\as(*s or decreases from a value of 2. At = 2, there is no 
compressibility (‘ff('ct for any value of ( r/ITr ) 2 . Tor values of less than 
2, the recpiired solidity decrc'ases with increasing values of (F/Fcr) 2 - 



Finciu-: 4-S. - Effect of comjiressifiility on axial solidity. 
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For D. values of more than 2, a region that is only of academic interest 
because it is beyond the limits of good design practice, the solidity ratio 
increases with increasing (V/Vcr) 2 . Experience has shown that D, values 
should be maintained below about 2 to avoid excessive losses. 


Relation of Loss to Solidity 

It is well recognized that the loading of a turbine blade or of a com- 
pressor blade is an important function of both solidity and reaction. 
Correlation of blade loss with a compressor diffusion parameter was 
described in reference 2 and is used widely within the compressor field. 
This parameter includes two terms, one reflecting reaction and the second 
reflecting turning and solidity. An analogous diffusion parameter was 
evolved for the case of the turbine in reference 3, where an overall diffusion 
parameter is defined as the ratio of the sum of the decelerations in kinetic 
energy on the suction and pressure surfaces to the exit kinetic energy. If it 
is assumed that the pressure surface minimum velocity is low enough to 
neglect (Fp.„,in = 0), then the overall diffusion parameter is defined as 


D = 


V 


2 

mnfx 


~V2^+Vi^ 


V2^ 


(4-21) 


With the use of the definitions of D, (eq. (4-12)) and R (eq. (4-15)), 
equation (4-21) reduces to 


D = D,-R (4-22) 

As seen from equation (4-13), 

= (4-23) 

Substitution of equations (4-23) and (4-14) into equation (4-22) then 
yields 


2 cos «2 . 

D = — sin Aa — R (4-24) 

(Tx4' cos ai 

This relation is like that for compressors, with the two terms involving 
reaction and solidity. 

Attempts have been made to correlate turbine blade loss with both 
overall (ref. 4) and suction-surface (ref. 5) diffusion parameters. A 
definite trend of increasing loss with increasing diffusion was established, 
but complete correlation could not be obtained. Such a correlation of 
blade loss with diffusion parameter alone would not be expected, since 
different values of reaction and solidity giving the same value of D do 
not give the same loss. 

Consider first the effect of reaction on loss, as shown qualitatively in 
figure 4-4 (a). As reaction is reduced from a relatively high value near 
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unity, there occurs a gradual increase in blade loss. Further reductions in 
reaction to negative values cause the loss to increase rapidly. This varia- 
tion in loss with reaction is caused by the change in boundary-layer 
characteristics (which are discussed in chapter 6) as the nature of the 
flow varies from highly accelerating to diffusing. The negative reaction 
regime, although desired in many applications, is usually avoided because 
of the high loss encountered when conventional blading is used. 

The effect of solidity on loss is indicated in figure 4-4 (b). A minimum 
loss occurs at some optimum solidity. As solidity increases, the amount of 
frictional surface area per unit flow is increasing. As solidity is reduced, 
on the other hand, the loss per unit surface area is increasing because of 
the increased surface diffusion required. A minimum loss occurs as a result 
of these opposing factors. The value of the suction-surface diffusion 



Axial solidity, 


(a) Reaction. 

(b) Solidity. 

Figurk 4-4. — Loss trend with reaction and solidity. 
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parameter corresponding to the optimum solidity is a function of many 
factors such as Reynolds number, shape of suction surface velocity dis- 
tribution, and rate of turning. In general, as mentioned previously, values 
not exceeding about 2.0 arc used. 

Selection of Optimum Solidity 

Both analytical and experimental attempts have been made to identify 
optimum solidity. According to reference 1, minimum loss occurs when the 
Zweifel loading coefficient ypz is equal to 0.8. By using this value in equa- 
tion (4-14), optimum axial solidity can be determined as a function of 
the blade-row inlet and exit flow angles, and this is plotted in figure 4-5 (a) 
for a wide range of angles. The dashed (long-short) curve represents the 
locus of points for impulse blading. 

In order to determine the optimum values in terms of actual solidity, 
it is necessary to determine the stagger angle because 

(T ^ — - — (4-25) 

cos a. 

An analytical blade model was used in reference 6 to relate stagger angle 
to the flow angles and the axial solidity. Thus, optimum values of actual 
solidity were obtained as a function of inlet and exit angles, as shown in 
figure 4-5(b). The authors of reference 6 compared an optimum solidity 
determined in this way with the data of reference 7, where efficiency was 
measured with four different rotor solidities, as shown here in figure 4-6. 
The solidity determined as optimum in reference 6 from a figure such as 
figure 4-5 (b) is seen to be quite close to that yielding maximum efficiency 
for this case. 

Loss coefficients based on cascade data are presented in reference 8 as a 
function of pitch/chord ratio (inverse of solidity) and exit angle for 
reaction blades (ai = 0) and impulse blades (ai= — aa). These coefficients, 
in relative terms, are replotted here in figure 4-7 against solidity for 
various exit angles. These curves indicate the importance of selecting 
optimum solidity. For the larger (more negative) values of exit angle, the 
curves arc rather flat in the region of minimum loss, and some deviation in 
solidity from optimum does not cause any significant increase in loss. As 
the exit angle gets smaller, the minimum loss region becomes more pro- 
nounced and the loss penalties become more severe as solidity departs 
from the optimum value. It must be recognized that curves such as those 
of figure 4—7 are usually obtained by using a given blade shape and varying 
the spacing. Thus, the blade shape and resultant velocity distribution 
cannot be optimized for each solidity, and the significance of such a 
correlation is somewhat clouded. 


13 







m.ADK DKSICN 



Fkjukk 4-0. -Viiriation of officiiMicy witli solidity for four turbines ol reien*nre 7. 



Kkuikk 1-7. Kffeet of solidity iuid exit nngle on hbule-loss eoeffieUMit . 


The optimum solidities obtained from the easeade results shown in 
figun^ 4-7 are ])lott(Hl against (‘xit angh‘ in figim^ 4 -S and an' (‘omjian'd 
with those ohtaiiK'd analytieally and shown in figure 4-5(1)). It is obvious 
that agreement betwec'ii the (‘xperimental and the analytical n'sults is 
not good for most exit-angle values. Although th(^ ('Xj)('rimental and the 
analytical curves do cross each otluT for both tlu' reaction («i = ()) blading 
and the impulse (ai= -ao) blading, the indicattal variations in o])tinuim 
solidity with c'xit angle an^ just not similar. All that can be said at this 
tim(‘ is that the analytical results involve many assumptions, th(' ex]ieri- 
mental results j)('rtain to ori(* particular l)lad(' profile', and there' are' many 
factors tliat act te) de'termine eeptimum se)lidity in a manne'r that we do 
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not yet fully understand. Analytical n‘.sults, suc'h as those figure' 4 5, 
ar(‘ moi'(* frenjiaaitly us('d te> d('t(‘rniine o|)1imuin solidity ifian arc* (*xp(*ri- 
mental I’esults, such as those* shown in figure* 4 7. < ’urreiit d(*sign ])ractic<! 
is to use* valu(‘s {)f O.U to 1.0, which is slightly high(*r than tlie* O.S 
r(*connne*nd(Ml iri lefere'iice* 1. 

Ultralow-Solidity Blading 

Ifi the* past, the limitation to re'ductions in solidity has h(*e*n se*paration 
oce'urriiig on tla* suction surface of the* hlade*. 4'o a<-hiev<* low(*r soliditie*s, 
some* modifi(*at i(U' in i)lad(* conce*pt must Im* utili/cd siudi that se*paration 
is su})[)re*sse‘d and the* assoe‘iat(*d high Iossi*s do not oce‘ur. 

44ie* tre*atnient of the* l)e)undary lay(*r in the* re'gion of s(*])aration is e)iie* 
a})])roae*h to r(‘du(‘<*d solidity. Siiedi t i e*a1 me*nts euaihl include* re*moving tin* 
i)oundary la\e*r ley siu'tion, energizing the* Ixuuiidmy la\er by blowing, or 
incre*asitig the turl.eule'ne'e* e>f the* l><»undary laye*r by use* of tiirbulators e>n 
the* blade*. (\*rtain of the*se e‘e)m‘(‘pts have* be*en (‘xplore-d with marginal 
succe'ss. dhve) alternate* l)lade* cone*e‘[)ts 1 hat have*, i)(*rha])s, ])e*ttc*r ])ot(*ntial 
an* tfie* lande*m and j(*t-fla[) blaele*s, w liiedi are* lllust r at e el in figure* 4 0. 
Studie's ap|)lying tlie* boundary-la\ e*r tr(‘atm(*nt e*ein(*e>pls as w(*ii as the* 
alt(*rnat(* hlade* e*e>nce*])ts to stator blael(*s anel rotor hlade*s are* summariz(*d 
in r(‘tVre'nce*s t) aTid 10, n*spe*e*t ively. ( 'ascade* tests eef leew-seelidity plain, 
tandi'in, and je*t-flap blade's are* j>re‘se*nte*el in r(*fe*re*neM*s 1 1 to 14. ddirbitie* 
test rc'sults with low-solielity tand(*m and j<‘t-flap rotors are* pn'se*nt(*(.l in 
refe*n*ne*e*s lb and Ki, r(‘sp(‘ctivc‘ly. 

4di(* tande*ni lelaele* op(*rat(*.s on the* prine*i[)Ic* t hat , all he sigh a high value 
e)f suet ion-surfaea* dirfusie)n is utilize'd (pei'haps 2). the* front foil is ter- 
minati'el at altout the* point of se*parat iejin The* I’eunaining diffusion th(*ri 
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Tandem blades 


Jet-flap blades 



r'Kinu-; 1 9. I.ow-snlidit \' runrc'pfs. 


tiikc's |)lac(' on tlu' r(‘ar foil with a t‘l(‘aii }>f)n]HlMry layor and witli pcrlinps 
20 to do of t]i(‘ inaiuHtroam air i!;oin^ tliroiigli tin* slot . 

TIh' j(‘t-flap l)lad(^ opci-at ('s with a S(‘con(l:n‘v air strc^nm jHliipu; out tlu^ 
trailing (*d»;(‘ p(“rp(ai(li(‘u]ar to tlu' main stronm. ddiis jet mov{‘s tin* rear 
sta^nat ion ])oint around t lu‘ I r ailing <'d<i'(‘, I }i(a (‘1)y suhstaiit ially iii(*r(‘asin^ 
tli(‘ lift. In addition, tln^ j(‘t d{'liv(‘rs souhv for(*(‘ to th(‘ blade throuuli its 
own ua)nH‘ntum. kd»ur(' 1 10 shows (‘XperinnadMl v(‘lof‘ity distributions 



i‘’i<a'Ki: t-10. ,I(d-flap (‘xperinu'ntal v(‘lo(‘ify (list l il mt Inns. 
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around on(‘ such blade witli the j(d ofT and on. Witli tlu^ j(d on, tluTO is no 
lon^('r arcciuir(^ni(‘n1 for th(’ velocities on tlu' suet ion and prc'ssure surface's 
t o lx* (‘(pial at t h(* blade trailing- (‘d^e*. Idu* loading diagram now approaches 
a r(‘(*tangular shape, with tlx* load coe'iiicie'iit \p inor(* close'ly a])))roachinjr 
unity. Also, th(* diffusion on tlu* suction surface* is sul)st antially r(*duc(*d, 
tlius suppressing' tin* t(*nd(‘n(*y to s(‘parat(*. 

Both tlu* j(*t-flap-blade and the tand(‘in-])lad(* conc(*pts otter tin* ])oten- 
tial for solidity r(*du(*tions. dlu* j(*t fla]), how(‘V(*r. will ])rol)abiy lx* con- 
sld(‘i{*d only for ap])licat ions wh(*n' a s(*condary air flow is napiired for 
oth(*r pur[)oses, sueh as blade* coolinj;'. 


BLADE-PROFILE DESIGN 

After tlx* blade* (‘liord l<*n«;th has lx*(*n s(*l(‘ct(*d and the* blade spacing 
d(*t(‘rniin(*d from solidity (‘onside'rations, tlx* blade itsell must lx* d(‘sign(*d. 
'This involves (h*t (*rmination of tlx* inl(*t and (*xit g(*e>m(‘tri(*.s and the 
comx'cting surface* profiles. Tlx* inl<*t and exit parts of the blade must 1)0 
d(‘sigrx'd to provide a snxK>th, (*fh(‘i(*nt transit ie)n b(‘tw(*(*n tlx* blade 
chaniH*! arxl tlx* fr(*(* str(*am. d’lx* surface* ])rofile*s conne*cting the* inlet nm\ 
(*xit must ])rovid(* t lx* re*(iuire*el Wow t urning wit h minimum loss. 

Exit 

ronsid(*ration of the blaele* e*xit s(*e*tion ineluel(*s tlx* trailing e*dge, the 
thn>at, and tlx* sue‘tion surfae*e* b(*twe*(*n the* throat [uxl tlx* trailing e*dge*. 

Trained — In tlx* de*sign of turbirx's, it is wise* te> utilize* the* smalk'st 
trailing e*dge* (‘onsiste*nt with nxx'hanical ce)nsid(*rations. As shown in 
re*fe*n*nc(* 17. an in(*n*ase* in t railing-e*dge thie*kix*ss <’aiis{*s an incre*ase in 
the* blade* loss, 'rhis e*ff(*ct is elise*usse*d furthe*r as part of tlx* turbitx*-loss 
eliscussie)n in chapte«r 7. In addition, t railing-(‘dge* thie-km*ss alse> has a 
sigrufi(*ant (*(Te*e*t on tlx* fle>w blex'kiige* in tlx* blaele* e*xit re*gie)ii. 

( \>nsid(*ration of the* blockage* e*lTe*ct will l)e* maele* with the* use* e)f figure* 
4 -1 1 , whi(‘h shows example* lelade* s(*ctions with the* ne)nx*nclature* use*d. A 
ne*w e*xit-va*le)cit v diagram is ceenst rue*t(*d at station 2ep which is hx*ate‘d 
just within the* blaele* t railing-e*dge* re*gion. 44x* re'duc(*el are*a d\ie* to tlx* 
t railing-e*dge ' bloe'kage* re*sults in a higlx'r ve'loe'ity at station 2a than at 
station 2, whie*h is lex*at(*d ju.st Ix'vond tlx* blade* trailing-e*dge* re*gion. 
44x* e*eiuations that have* lx*e*n use*d te> obtain this ‘Avithin-tlx'-blade*'’ 
diagram at 2(t irx*hxle* e*onse*rvatie)n of tange*ntial inom(*ntum: 

and continuity : 

{pVs}2.. ( I — ) = (pe;)2 (4-27) 

\ .s (*e)S a 2a/ 
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4 11. Hljide section and nonuaiclatiii'c. 


^vho^e t is tho trailirig-odo;o thickness, in ni(‘t('rs or f(‘(d. The flow an«;I(‘ ao^ 
is cl(d('rminecl from (‘([nations (4“2()) and (4 27) by assuming; th(‘ flow 
h(‘tw(‘(‘n stations 2a and 2 to be (‘ith(‘r in(*ompr(‘ssi})](‘ (sinc(‘ tlu‘ (dianj>;es 
are usually small) or is('ntro[)ic. Th(‘ l>lad(‘ must b(‘ (hvsi^ned to have* an 
exit an^l(‘ of a^a in ordca* to [)roduc(‘ a va‘locity-dia<i;ram angle at station 
2 outside th(‘ })lad(‘ row. 

1 he Ma(‘h iuimb(‘r at station 2a (*an also b(‘ d(‘t(‘rmiii(‘d from tlu‘ 
pr(‘C(‘ding (‘([nations and assumptions. B(‘caus(‘ th(‘ angl(‘ ao is oft(‘n large* 
((>5° or gr(‘at(‘r) and th(* flow Mach miml)(‘r at tiu' blade* (*xit (station 2) 
is ott(‘n sp(*cifi(‘d to be* in the* liigh subsonic re'gion, the* trailing-(‘dge 
})lockag(‘ can c‘ause station 2u to Ix'conu* edioked. It is, lh(*r(‘fore*, imjiorlant 
to de‘t(*rmin(‘ wh(*th(*r edioking inside* tlu* blade* row will occur su(‘h that 
the* d(*sigri flow rate* cannot lx* obtain(‘d. 

Throat.~^\\Wi\ in ge*n(‘ral, a turbine* blade row o[x*rat(‘s as a nozzle*, wit h 
the* flow acce*l(*rating ii[) to the* throat, or minimum are*a, the* (l(*t('rmination 
of tile* throat ope*ning o (see* fig. 4-11) become's a nithe*r critical aspe*e*t of 
the* d(*sign pro(‘e'dure*. One* te*chni(|ue* use*d succe‘ssfull\’ to give* this dime*n- 
sion make*s use* of the* “inside*-the*-f railing-(‘elge*” v(*loe*ity eliagram. If one* 
assumes ne> cliange* in flow conditions anel a straight suction surface* 
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lH‘t\v(M'u ili(‘ thrniit and stali<»n 2a, tluai tli(‘ thi’oal dimcaisiori can 1)C 
obtained from tin* volo(*ity diagram at station 2a l)v using tlu' following 
Of [nation : 

1— ■ ) cos {T 2 „ (4-28) 

N cos 

wh(‘r(‘ a is th(' tliroat opcaiing, in m('tcrs or f(M^t. 

If it is assunu'd lhat the velocity and loss do ikjI ehang(‘ l>etw(‘(‘ii tli(‘ 
throat and the “fiMK^-strcam” station 2, tlaai 

c<)s a-j (1 2!)) 

.V 

W'hen thi< nu'thod is used, the <4tect of 1 railing-(‘dge lhickii(‘s.s (dmnges 
th(‘ angle of the throat position but not its haigth. Hoth imhhods (('([S. 

( 1 2S ; and ; 1 2h ' ) give similar 1 hmat dinuaisions. I b f(‘r(‘nc(‘ S compar(‘S 
iu(‘asur(‘d ('xit~fl<A\ angles with lhos{‘ pr(Mli('t<“d by (Mjuation (4 29). 'This 
comparison i!idieat('s clos(* agr(Maiu'nt at (‘xit angh‘s gix^ater than ()0° and 
d(‘viaiions of up to o'^ Pu’ (‘xit angles down to M-v'. This (hwaation could be 
du(“ to low('r solidities as w('ll as larger gradi(‘iits that would occur across 
t h(‘ t hroat , 

4'lir throat'Op('ning dinaaisioii as (ha (aniim'd fi“om (M[uati()u (4 2S) 
or i t 29 I applies to the cas(> w here the bhuh‘-row (‘xit flow is subsonic. 
If th(' flow within the l>lade row (‘xpands to a sup(M’soni(‘ v(T)citVj tluai 
t his com])ute<l t hroat dini(Misi< t!i must bi' moditi(‘d t o account for (‘xpansion 
from th(‘ sonic f'Miidition at tli<‘ throat to tla^ super‘sonic condition at th(‘ 
('xit. k'or (‘xit Maeh numb(‘rs great er t han aljout bib t h(‘ ch(4\ing sectioti 
( throat ) must i)c lo(‘at{‘d back witliin th(' channel smTi t hat a cf)nv(Tg(ait- 
(liv(‘rg(‘nt passag(‘ is obtaiiK’d. l or low suj)ersonic Mach numl)crs (uj) to, 
p('rhaps, 1.0 ). it has Ixaai found that satisfa<eory p(‘rfoniianc(' can b<" 
a(4ii('ved if tlu‘ tlirojit is still loeat(‘d at ilx' exit of ihe (*haniieh and tlu' 
additional flow (“xpaiision occurs dtw^ list I’eam liom tin* thixmt. In this 
cas(\ th(‘ rcfpiired f liannel (‘xit dimension o would In- computc(l by tlu‘ 
f< illow ing juat ion ; 

whi'va 

(r., throat f^p^‘'dng (‘om])uted from (‘(juation (1-28) or {4-29) for 

siiper<oi]ic v(Tx-it y. m ; ft 
A,r flew an 'a f<a’ soni(‘ flow, m-: ft- 

.1,,- How ai’('a for sup(‘rsorii(‘ flow. m~: ft- 
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I'Mjrui, } 12. \':u'i:it inti iti flow willi <ni)(‘rsnnir flow Marh iiimu1h>i'. 

This ‘diva (*()iT(a‘t w itli assunuHl is(‘ntropi(* flow In'twcaai tliroat and 
('xit, is shown in figure \ 12. 

Surf /Off sin fdcr doirfisfrram fn^fn (liroat . — Th(‘ S(T‘ctinn oi“ (h<‘ typ(' of 
surlacc- l)ctw(H‘n tlu‘ throal an<l trailing on th(‘ su<‘tion surface' must 
})(‘ ina(l(‘ Ironi such (‘onsidc-rat ions as stnudural intc^i-ity in tlu^ trailinj>;- 
cdj>;(‘ rc^^ion, Ma(‘h nuiniKa’ and asso(‘iat(‘d loss(>s, dcsirt^d of 

suet ion-snrla(*(‘ diffusion ( />.,), and hlad(‘ surface' are a re'sultin^ fr*e>in tlie' 
design. 

A “straight hack'’ ek'si^ti is use'd wIk'Ii le>w value's e>f ( apf)roximatcly 
unity ) are' spe'cilie'd anel le)ii”; t railing e'dge's are' pe'rinissihle'. sul)Soni(‘ 

or transonie* hladinj^, as wejuld he' indicate'd by the' dis(‘ussie)n in llu' next 
para^ra[)h, use's this type' of surface' in oreh'r te> ])re've'rit flow acce‘I('ratie)n 
on the' tail of the' blade- anel ke'e-p tlie' assoe*iate'd le)sse's Ie>w. Principal 
pre)bl(‘ins with a slrai^^ht surface' are that the le)W J), value-s j)re'clude' le)w- 
se>lielity de-si^iis anel the- le>n»; trailing e'dge- e*an be-conu' structurally flimsy. 

.\le)st e'on ve-nt ie)nal gas-turbine' blaeling utilize's some amount of curva- 
ture' be-twe-ere the- throat anel t railing-e-elge- re-gion. This ])e'i'mits some- dif- 
iusie>n and aelelit ional leeaeling e)H t he- tail e>f t he- blaele', anel it aelels conside-r- 
ably to the- strue-tural inte-grity e>f the- blaele- by irit retdue'ing a we'dge' angle* 
at tlie- ('xit. If e-e)n ve-nt ie>nally leeade-el blaeling is use-el, the* e-ffe-ct of this 
curve-el surlae-e- on loss is not gre-at. As inelicate-el by figure- 4 Id (\N’hich is 
fmm I'e'l. S), it the- e-xit-fle)\\' Mach numf)e-r is le-ss than 0.8, tlie- e-urvature- 
e-fte-e't is small. At highe-r e-xit Mae*h numbe-rs (gre-ate-r than 0.8), the' e-ffe-ct 
on kiss e*aie be-e-ome- se-ve-re-. The‘re*fore‘, design e'lirvature-s she)ukl be- lowe-r in 
the- highe-r Mae-h numbe-r re-gieins. The* type- e>f curvature' sc-le-e'te-d fe>r the- 
suctie>n surface- be-twe-e-n the- tliroat and trailing e-elge- has an e-ffe-e-t on the* 
suet ie)n-surfae*e' ve-leicity (list ributiem. In gene-ral, the- velocity distribution 
is imjireive'd if the' ('urvature' de-cre-ases from threiat te) trailing enlge- inste-ad 
of re-maining cemstaiil . 
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Fnanu: \ V:tri:ili«m of profile loss with Much nuinlMT and surfin-e curvature 

l>et\veeu throat and exit (from ref. S). 


Inlet 

1 h(* l(‘julin^-<*d^(‘ ^('onK'try ol ti turl)iii(‘ hhi(i(‘ row is iisutilly It'ss cTitictil 
than tlu- cxit-K'Kion geometry. At the l)I:ule itilet, a relatively largi; 
leading-edge radius can usually he used, heeause tin- .Maeh mimher is gen- 
erally low at the inlet and then increases through the hlade row. '1 he lead- 
ing edge l)('(*oin('s a serious eoiua'rn for low-r(aiet ion hlading and high .Maeh 
number lihuling. In the eas(‘ of low-ree-tioii hlading, excessively high 
vi'loeities in the inlet region can lead to high valui's of suet ion-surf aee 
dilTusion and a t('iideney toward inerea.sed losses. With high inlet .Mach 
numhers, care must he taken tied tlie area contraction is not so severe as 
to choke the hlade at the inlet. Kciuations (4-‘2()) and (4 27), which were 
u.sed for the hlad(‘ exit, can also he us('d to determine a hlade-inlet opening 
and “wilhin-the-hlade” flow angh- and .Maeh mimher to cheek for hlade- 
inl(d (‘hokiniz;. 

Altliou^li circular l('a(Iin»; an^ \isually this is arl)itrary 

and could limit llu' frc(‘dom of vchxaty-distrihution s(d(‘ction in th(‘ 
h‘adins-(‘d^(‘ r('»;ion. TIk^ lar^(^ (uirvatun'S associatfsl with circular hmlinj^ 
(‘dj*:t‘S can n^sult in undcsirahh' v('locity pcaiks on both th(^ suction- and 
pr(*ssur(‘-s\irfac(‘ ])ortions oi th(* leading (‘dgc. Otlu'r gconi(dri(‘s, su('h as 
cllips(‘s, which pfuanit variat ions in curvat un* around t h(‘ hauling (‘dgc, can 
be us<‘d to minimiz(^ or (‘liminatc^ tlu* v('locit\' p(^aks. 

Blade-Surface Profile 

Onc(‘ th(‘ li'ading- and trailing-edg(‘ g(M)iiu‘tri('s have Ixsai s(‘l(‘ct(‘d, the 
task rc'inaining is to join them witli a prolil(‘ tluit yi(‘lds tlu' r(‘(|uir(‘d flow 
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turniiig and a satisfactory velocity distribution around the blade. The 
design procc^dure must describe the flow conditions through th(? blade rows’ 
to an accuracy sufficic^nt to impose design controls (e.g., diffusion limits). 
Two of th(' major flow considerations are illustrated in figure 4-14. 
Velocity gradients occur across the channel from the suction to the 
pr(\ssur{‘ surfac<‘ as a result of th(» static-prc'ssure diffcTence n'ciuired to 
turn the flow. Radial variations in streamliru' position and, therefore, 
velocity occur as a result of radial-equilibrium considerations. Since both 
of th(‘se factors influence th(^ blade-surface velocity distribution, the design 
procedures us(‘d should l)e at least of a (piasi-threc'-dimcmsional nature. 
The channel flow analysis theory that s(tv(^s as the basis for these design 
procedures and the computer programs available to perform the com- 
putations arc discussed in the next chapter. 


Pressure r Suction 
surface 7 \surface 



(a) 



(a) Cross-channel variation. 

(b) Radial variation. 

Figukk 4-14. — Turbine blade-row velocity variations. 
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SYMBOLS 

A flow area, m^; 

c chord, m; ft 

D diffusion parameter 

F force, N ; lb 

g conversion constant, 1 ; 32.17 (Ibm) (ft)/(lbf) (sec^) 

K ratio of inlet to exit tangential components of velocity ( F«,i/F„, 2 ) 

0 throat opening, m; ft 

p absolute pressure, N/m^; Ib/ft^ 

R reaction 

s blade spacing, m; ft 

t trailing-edge thickness, m; ft 

V absolute velocity, m/sec; ft/sec 

X axial distance, m; ft 

a fluid absolute angle from axial direction, deg 

a, blade stagger angle from axial direction, deg 

y ratio of specific heat at constant pressure to specific heat at 

constant volume 
p density, kg/m^; Ib/ft* 

(T solidity 

loading coefficient defined by equation (4-6) 
loading coefficient defined by equation (4-5) 

Subscripts : 

cr critical 

inc incompressible 

max maximum value 

min minimum value 

opt optimum 

p pressure surface 

s suction surface 

55 supersonic 

u tangential component 

X axial component 

1 blade row inlet 

2 blade row exit 

2a within trailing edge of blade row 

Superscript : 

absolute total state 
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Channel Flow Analysis 


By Theodore Katsanis 


The design of a proper blade profile, as indicated in the last section of 
chapter 4, requires calculation of the blade-row flow field in order to 
determine the velocities on the blade surfaces. This chapter presents the 
analysis theory for several methods used for this calculation and also 
discusses associated computer programs that were developed at NASA 
Lewis Research Center. 

The actual velocity distribution throughout a blade-row flow field 
cannot be calculated at this time because of the extreme complexity of 
nonsteady, viscous, three-dimensional flow through geometHcally complex 
passages. To calculate a theoretical velocity distribution, therefore, certain 
simplifying assumptions must be made. The three-dimensional flow is 
simplified to flow on or through various two-dimensional surfaces. Such 
surfaces are illustrated in figure 5-1 for the case of a radial-inflow turbine. 
Similar surfaces are used for an axial-flow turbine. A flow solution on the 
mean hub-to-shroud stream surface (commonly called the meridional 
surface), shown in figure 5-1 (a), does not yield blade-surface velocities 
directly, but provides information required for the blade-to-blade surface 
(fig. 5-1 (b)) and orthogonal surface (fig. 5-1 (c)) solutions, which yield 
the desired blade-surface velocities. 

There are two parts to a method of analysis to obtain a velocity dis- 
tribution over one of these surfaces. The first part is the mathematical 
formulation of the problem, and the second part is the numerical solution 
of the mathematical problem. For the mathematical formulation of the 
problem, we will discuss stream- and potential-function methods and 
velocity-gradient (stream-filament) methods. The stream- and potential- 
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Hub-to-shroud Blade-to-blade 



(a) (b) 



(c) 


(a) Hub-to-shroud stream surface. (b) Blade- to-blade surface, 

(c) Orthogonal surface acros.s flow passage. 

Figure 5-1. — Surfaces used for velocity-distribution calculations. 
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function methods will be described relative to the blade-to-blade surface 
solution. A similar type of analysis can be made for the meridional surface. 
The velocity-gradient equation to be presented is general and can be used 
for solutions on any of the surfaces. 

The following assumptions are made in deriving the various methods 
of analysis discussed herein: 

(1) The flow is steady relative to the blade. This means that the 
surface velocity at any given point on the blade does not vary with time. 
Thus, if the blade is rotating, the flow would not be steady relative to a 
fixed coordinate system. 

(2) The fluid obeys the ideal-gas law 

V = pRT (5_1) 

where 

V absolute pressure, N/m“; Ib/ft* 

P density, kg/m®; Ib/ft® 

R gas constant, J/(kg) (K) ; (ft) (lbf)/(lbm) (°R) 

T absolute temperature, K; °R 

or is incompressible (p = constant) . 

(3) The fluid is nonviscous. A nonviscous fluid has no boundary layer. 
The blade-surface velocity is calculated, therefore, as if the free stream 
extends to the blade surface. 

(4) The fluid has a constant heat capacity. 

(5) The flow is isentropic. 

(6) The total temperature and total pressure are uniform across the 
inlet. 

(7) For the stream- and potential-function analyses, the additional 
assumption is made that the flow is absolutely irrotational. Therefore, 

curlV = VXV =0 (5-2) 

where V is the absolute velocity vector. Intuitively, this means that 
particles do not change their absolute orientation with time, although 
their shape may change. For example, figure 5-2 shows a hypothetical 
particle at times t and f-f Af. In the absolute frame of reference, the 
particle changes its location and shape at a later instant of time, but the 
net rotation is zero. Of course, in a frame of reference relative to the 
blade, the particle has rotated, because the frame of reference has rotated. 
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Direction of rotation 




(Absolute frame of reference) 



(Relative frame of reference) 


Figure 5-2.— Absolutely irrotational flow. 


Some numerical techniques for solving the mathematical equations will 
also be discussed. However, it must be emphasized that there are many 
techniques for solving these equations, and we will discuss only a few. An 
excellent theoretical discussion of flow in two-dimensional cascades is 
given in Chapter IV of reference 1. 

STREAM- AND POTENTIAL-FUNCTION ANALYSES 
Stream-Function Method 

The stream function can be defined several ways, but perhaps the 
simplest is in terms of streamlines. Suppose we consider two blades of a 
cascade as shown in figure 5-3. It is assumed that there is two-dimensional 
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Figure 5— 3. Streamlines for a stator cascade. 


axial flow here, so that the radius r from the centerline is constant and 
there^ is no variation of the flow in the radial direction. There may be 
rotation about the centerline. 

Shown in figure 5-3 are a number of streamlines. The mass flow between 
the blades is w. The number by each streamline indicates the fraction of w 
passing between the upper surface of the lower blade and the given stream- 
line. Thus, the upper surface (which is a streamline) has the value 0, and 
the lower surface of the upper blade has the value of 1, while the remaining 
streamlines have values between 0 and 1. Note that a value can be asso- 
ciated with any point in the passage. This value is called the stream- 
function value and can be used to define the stream function. 

It will be recalled that mass flow can be calculated for a one-dimensional 
(or uniform) flow by 


w = pVA 


( 5 - 3 ) 


31 



TURBINE DESIGN AND APPLICATION 


where 

w rate of mass flow, kg/ sec ; lb/ sec 

V fluid absolute velocity, m/ sec ; ft/ sec 

A flow area normal to the direction of the velocity F, m^; ft^ 

This can be extended to a varying flow by using an integral expression: 


w = 



dA 


( 5 - 4 ) 


Since this stream-function anal3^sis applies to both stationary and rotating 
cascades (blade rows), the fluid velocity will be expressed in terms of 
relative velocity Wj which for a stationary blade row reduces to absolute 
velocity F. We will assume that our cascade has a uniform height b. Then, 
the mass flow Wi,2 between any two points Qi and Q2 in the passage (see 
fig. 5-4) can be calculated by 


Wl,2 = 



(^ 5 ) 



Figure 5-4. — Arbitrary curve joining two points in flow passage. 
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where W„ is the relative velocity component in the direction of the rieht- 
hand normal of the line going from Q. to Q,. This sign convention means 
^at wi,i will be negative if Qj is below a streamline passing through Q,. 
The integral is a line integral between the points Qi and Q* and is in- 
dependent of path for steady flow relative to the cascade. 

With the use of equation (5-5), an analytical expression can be written 
for the stream function w at a point (x, y) : 


w(x, y) = 



(5-6) 


where Qo is any point on the upper surface of the lower blade, and the 
integral is taken along any curve between Q« and (x, y). This is indicated 
in figure 5-5. 

Since the integral in equation (5-6) is independent of path, it is rela- 
tively easy to calculate the partial derivatives of u. For example, we will 
calculate du/dx at the point (x, y). Let x«<x such that the point (x», y) 
18 still m the flow passage, as shown in figure 5-6. Then 


«(x, y) =- 


f pWnbdq+f pW„bdq 
•'ci Jct 


w 


(5-7) 



Figure 5-5.— Curve joining (x,y) with a point on the upper surface of the lower blade. 
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Figure 5-6. — Curve joining horizontal line through (x,y) with a point on the upper 

surface of the lower blade. 


where Ci is an arbitrary curve between Qo and (xo, y), and C2 is a hori- 
zontal line between (xo, y) and (x, y). The integral along C, does not 
depend on x. Along Cj, we have W n = — Wy and dq = dx. Hence, 


du 


d pWyb dx 

dx w 


(5-8) 


du pW yh 
dx w 


(5-9) 


In a similar manner, we can calculate 


du pWJ) 
dy w 


(5-10) 


Now we will make use of the fact that the flow is absolutely irrotational. 
From the definition of the curl operator and the above assumption, 


/dV, dVy\. /dV. dVA {dVy 


where i, j , and k are the unit vectors in the x, y, and z directions, respec- 
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lively, and F;c, Fj,, and F* are the absolute velocity components (in 
m/sec or ft/sec) in the x, y, and z directions, respectively. Since we are 
considering two-dimensional flow only, 

F, = 0 (5-12) 

and 


dV, dVy 

dz dz 


(5-13) 


Hence, equation (5-11) requires only that 


dVy dV, 

dx dy 


(5-14) 


Since 

Fx = TFx (5-15) 

and 

Vy = Wy+icr (5-16) 

where w is the angular speed (in rad/sec) and the radius r is constant, 
equation (5-14) can be expressed in terms of relative velocities as 


dWy dW. 
dx dy 


(5-17) 


Actually, the flow is irrotational with respect to the moving coordinates 
in this particular case. Now, from equations (5-10) and (5-9), 


w 

" pbdy 

(^18) 

W du 
~pb dx 

(5-19) 

Substituting equations (5-18) and (5-19) into equation (5-17) yields 

d /I du\ d /I du\ 
dx \p dx) dy \p dy) 

(5-20) 

since w and b are both constant. 

For incompressible flow, p is constant, and 


^ d^u 

V"w = — +— = 0 
dx^ dy^ 

(5-21) 


which is Laplace’s equation. Any function satisfying Laplace’s equation is 
called a harmonic function. There is a great deal of theory concerning 
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harmonic functions that is related to the theory of analytic functions of 
complex variables. 

The important thing to know here is that there are a tremendous 
number of functions that satisfy equation (5-21), and we must find a 
solution that satisfies certain boundary conditions. The solution to either 
Laplace^s equation (5-21) or equation (5-20) will be determined by 
specifying two things: (1) a finite region, and (2) a boundary condition 
along the entire boundary of the region. 

The first thing that must be specified is the solution region. A typical 
two-dimensional cascade is shown in figure 5-7. Since the flow is the same 
in every passage, we can consider a finite solution region as shown in 
figure 5-8. It is assumed that AH is sufficiently far upstream so that the 
flow is uniform along this part of the boundary and that the flow angle 
fiin is known. Similarly, it is assumed that the flow is uniform along 
DE, and that the flow angle ^out is known. From the way the stream 
function was defined, we can specify boundary conditions on the entire 
boundary ABCDEFGHA. Along BC, u = 0; and along FG, u=l. Along 
AB, HG, CD, and FE, a periodic condition exists; that is, the value of u 
along HG and FE is exactly 1 greater than it is along AB and CD. Along 
AH and DE, du/dr) is known, where 17 is the distance in the direction of 
the outer normal. 



Figure 5-7. — Two-dimensional infinite cascade. 
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Consider the differential of u in the direction of the velocity W : 

, du du , 

= = 0 (&- 22 ) 

The differential is 0 because the stream function is constant along a 
streamline, and the velocity vector must be tangent to a streamline 
Along AH, 


du 


du 


drf dx 

and substitution from equation (5-22) yields 


However, 


du du dy 
drj dy dx 


(5^23) 


(5-24) 


^ = tan^ (5-25) 

Further, du/dy is constant along AH, since it is assumed that the flow is 
uniform there. Therefore, 


du Cu(H)-u(A)] 1 
dy s ~s 


(5-26) 


where s is the blade spacing in the y direction. Substituting equations 
(5-25) and (5-26) inequation (5-24) gives along AH 


/ tan /3,„ 


(5-27) 
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Similarly, along DE, one can calculate 



We now have a boundary condition along the entire boundary of the region 
shown in figure 5-8. These boundary conditions will always determine a 
unique solution to Laplace’s equation (5-21). For compressible flow (eq. 
(5-20)), a unique solution is always determined if the flow is strictly 
subsonic throughout the region. 

There are numerous techniques for solving equation (5-20) or (5-21). 
After the stream function is obtained, blade-surface velocities and velocities 
throughout the passage can be obtained by differentiation of the stream 
function. This is what is known as the direct problem. A method of solving 
this problem will be discussed later. The indirect, or inverse, problem is to 
specify a desired velocity distribution on the blade surface and from this 
determine a blade shape that will give this velocity distribution. This will 
not be discussed here. 

Potential-Function Method 

For two-dimensional irrotational flow, a potential function can be 
defined. If lines of equal potential are drawn, they will be orthogonal to 
streamlines. The potential function will not be defined in the same detail 
as the stream function, but the main properties and relations will be given. 
If the potential function ^ exists (i.e., the flow is irrotational) , then it can 
be defined so that 


and 



(5-29) 


— =F, (5-30) 


We will refer to absolute velocities here, since we must have flow irrota- 
tional relative to the coordinate system used. This, coupled with the 
assumption of absolute irrotational flow, implies that the coordinate 
system does not rotate. This does not exclude use of the potential function 
for pure axial flow, since the rotation has no effect if there is no change in 
radius; that is, the flow is actuall}" irrotational with respect to the blades, 
as w^e saw in the discussion of the stream function. 

From the continuity relationship for steady flow, 


a(piT) ^ djpVy) 
dx dy 


(^31) 
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Substituting equations (5—29) and (5—30) in equation (5—31) yields 
d / d ( d^)\ 

If the flow is incompressible, p is constant, and 


524. 

V2$=—_ -I =0 

dx^ 


(5-33) 


So, the potential function satisfies Laplace^s equation. Thus, for incom- 
pressible, irrotational flow, both the stream function and the potential 
function satisfy the same differential equation (Laplace's equation). The 
difference lies in the boundary conditions. 

We can consider the same solution region shown in figure 5-8. We can 
specify boundary conditions over the entire boundary as follows* Alone 
BC and FG, 


dn) 

where Fn is the velocity normal to the blade surface. Along AH, 

(d.= 


and along DE, 


( 


/, 




The inlet and outlet axial velocities are given by the equations 


(F.),n = 


Pinbs 


(5-34) 


(5-35) 


(5-36) 


(5-37) 


and 


(F^)ouf = — ^ (5-38) 

Poutbs 

Along AB, GH, CD, and EF, a periodic condition exists. Since the flow is 
uniform along AH, 


Substituting 


/^\ _ [4>(H)-<I>(A)] 

WAn S 


Vy= Vx tan 


(5-39) 

(5-40) 
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into equation (5-39) yields 

^(H) = ^(A) +s{ V,) vn tan (5-41 ) 

Because of the periodicity, ^ is exactly s{Vx)<n tan /Sm greater along HG 
than along AB. Similarly, at the outlet, 

4>(E) = 4»(D) +s(yx)o«t tan fiout (5-42) 

Equation (5-42) gives the difference in 4> along the lines FE and CD. 

This completes the boundary conditions for equation (5-32) or (5-33). 
The boundary conditions, however, do not determine a unique solution, 
but only a solution within an arbitrary additive constant. If the value of 

is specified at one point, these boundary conditions will determine a 
unique solution to equation (5-33) for incompressible flow, or to equation 
(5-32), for strictly subsonic compressible flow throughout the region. 

As for the stream function, there are numerous methods for solving 
equation (5-32) or (5-33) subject to the preceding or equivalent bound- 
ary conditions. A method for solving the inverse problem of specifying 
the velocity distribution to determine the blade shape is described in 
references 2 and 3, 

Choice of Stream- or Potential-Function Method 

If the flow is steady, irrotational, and incompressible, there is little to 
choose between the stream function and the potential function. In this 
case, the choice is made on the basis of ease of solution for the boundary 
conditions (the differential equation is the same: Laplace's equation). 
However, if any of the three assumptions (steady, irrotational, or incom- 
pressible flow) is not applicable, then we may be restricted as to the 
choice of stream function or potential function. 

The existence of the stream function is proven from the continuity 
equation. For the stream function to be defined, the mass flow crossing 
a line between two points must be independent of path. This requires 
that the flow be either incompressible or steady. Some additional assump- 
tion is necessary for the flow to be unique. We used the assumption that 
the flow was absolutely irrotational, which turned out to be irrotational 
relative to the blade for the axial-flow case considered. However, other 
assumptions could be made for other problems. Another restriction on 
the stream function is that it can be defined only for two-dimensional flow. 
This can easily be seen since the stream function is defined as a percentage 
of mass flow between two points, and this is meaningless in three 
dimensions. 

The existence of the potential function can be shown if the flow is 
irrotational relative to the given coordinate system. This is necessary 
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because we must have equality of mixed second partial derivatives • that 
is, if ’ 


then 


dxdy dydx 


(5-43) 



(5-44) 


and the flow must be irrotational. A similar situation exists in three- 
dimensional flow; that is, the potential function exists only if the flow is 
irrotational with respect to the coordinate system being used. Finally, 
an assumption must be made to assure a unique solution. This can be 
done by using the continuity equation. 


Finite-DiflFcrence Solution for Stream-Function Method 

As stated before, there are many ways of solving various problems 
posed by stream-function or potential-function theory. We will consider 
in further detail the finite-difference solution of the direct problem for the 
stream function for the simplest case of steady, incompressible, irrotational 
flow. In this case, we must solve Laplace’s equation subject to the bound- 
ary conditions discussed in the section on the stream function. The method 
of solution for the potential function is quite similar, but with a lower 
rate of convergence for the finite difference solution. 

The first step is to establish a rectangular grid of mesh points in the 
region shown in figure 5-8. A typical grid is shown in figure 5-9. Then a 
finite-difference approximation to Laplace’s equation (eq. (5-21)) can be 
written at each mesh point where the stream function is unknown. A 
typical mesh point with four neighboring mesh points is showm in figure 
5-10. The point in consideration is labeled 0, and the four neighboring 
points are labeled 1 to 4, as showm. The distance between points 1 and 0 is 
denoted h,, and similarly, the other distances are h^, h, and h, as indicated 
in figure 5-10. The value of u at points 0 to 4 are labeled uo to u,, respec- 
tively. With the use of a Taylor series expansion for u in the x- and y-direc- 
tions, equation (5-21) can be approximated by using only values of u at 
mesh points. (Further explanation of this is given in ch. 6 of ref. 4.) 
When this is done, the follow ing expression is obtained; 


r 2wi 2ui 

2uq 

1 2uz 

2m4 

2mo1 



L^3(^3 + ^4) 


A3A4J 


(5-45) 
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Solving equation (5-45) for Uo yields the expression 


where 


4 


Uo= X) 

(5-46) 

^3 + ^4 

a»hi 

(5^7) 

02 = , 
O 0 O 2 

(5^8) 

03 — 

O 0 /I 3 

(5-49) 

h\-{-h2 

ai- , 

O0/I4 

(6-50) 

(hi) 

(5-51) 


Equation (5-46) holds at every interior mesh point. If one of the neigh- 
boring points is on a blade surface, then the value of u at that point can be 
used. At other points along the boundary, equation (5-46) cannot be used, 
but the boundary conditions can be used to obtain alternate equations at 
these points. For example, along the upstream boundary AH in figure 5-9, 
duldr) is given by equation (5—27). If point 0 is on line AH, then, a finite 
difference approximation gives 

u.-u.+lu(^) (^ 52 ) 

Similarly, if point 0 is on line DE, 

= (^53) 


For the points along AB and CD, equations can be derived by using 
the periodic boundary condition. If the point 0 (fig. 5—11) is on the 
boundary between A and B, the point 1 is outside the boundary. However, 
it is known that ui = mi,b— 1, where the point l,s is a distance s above 
point 1 in the y-direction, as shown in figure 5-11. Substituting this 
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Figure 5-11. — Mesh point on line AB. 
condition in equation (5-46) gives 




• GiUi^5 + 2. d\Ui—a\ 
i-2 


(5-54) 

This equation holds along CD (fig. 5-8) also. 

points along HG need not be considered, since they are just 1 
mesh li corresponding point along AB. The equation for the first 

HG In^^th'^ modified, since point 2 is on line 

beL no nr2T’tlT^“''7'^^’ ^ s 

SnhsTit t .1 the negative y-direction, as indicated in figure 5-12 
Substituting this condition in equation (5-46) gives 

Ko = aiUi+ajM2.-.d-a3U3-|-a4U4-|-a2 (5-55) 

This equation also applies to the first mesh line below FE (fig 5-8) 

One of equations (5-46) or (5-52) to (5-55) can be appHed to each 

hlteL^t°r unknown in the region of 

knoZs Th**'''’ number of linear equations as there are un- 

rTrrTcU r '^'^known will be 

referred to simply as unknown mesh points. 

f there are n unknown mesh points. We then have n equa- 

^Th^v^r if consecutively from 1 to 

and so forlTim T P»“t, 

annlv The the last point. At each point, one equation will 

pp y. The equation at a typical point, i, could be written 
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^aaUi = ki (5-56) 

j-i 

The values of the are determined by one of equations (5-47) through 
(5-55). All but five, at most, of the aij are zero, and the aii= — 1. The 
value of ki is always zero, except for the outermost unknown points 
around the boundary. It can be shown that the a a matrix is always non- 
singular; hence there is always a unique solution for the Uj. 

A numerical solution to equation (5-56) can be obtained by iterative 
techniques. These techniques are particularly valuable in solving systems 
of linear equations of this type ; that is, where there are a large number of 
unknowns, but few terms in each equation. Storage requirements are 
small, and roundoff error is minimized with iterative methods. To start 
the iteration, an initial estimate of u at every unknown mesh point is 
required. The simplest iterative procedure is relaxation. This consists of 
changing the estimated value of u at each point in succession so as to 
satisfy the equation for that point. After this is done at every point, the 
procedure is repeated until there is negligible change in the values of u. 
The procedure is simple and it always converges for this problem. How- 
ever, the convergence rate is extremely slow, so that excessive computer 
time is required. The convergence can be accelerated greatly by increasing 
the change in u at each iteration by a factor w, called the overrelaxation 
factor. When co= 1, the procedure is straight relaxation, and when cj>l, 
it is overrelaxation. It is proven in reference 4 that overrelaxation (or 
underrelaxation) is convergent if 0<o?<2. However, the greatest rate of 
convergence occurs when Kco<2. In fact, there is an optimum value of o) 
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between 1 and 2 which gives the most rapid convergence. This optimum 
overrelaxation factor can be calculated as explained in reference 4. 

To give an explicit expression for the overrelaxation procedure, we will 
use a superscript on the uu That is, Wi"> is the m^h iterate of m. The initial 
estimates are denoted Ui® and may be any value. For example, an initial 
estimate of Ui» = 0 is satisfactory. Then, if Ui™ is known for all i, we can 
calculate for i= 1, 2, . . . , n in succession by 




Ui"' + CO 


— 

j-i 


aijUj^+fci 

j-i+I 



(5-57) 


After a solution for it is obtained by overrelaxation (or any other 
method), it is necessary to calculate the velocities with the use of equa- 
tions (5-9) and (5-10) as 



and 

Wy= (5-59) 

po 

The partial derivatives du/dx and du/dy must be estimated from the 
calculated discrete values of Ui. This can be readily done, either by finite 
differences, or by fitting a smooth curve, such as a spline curve, through 
the points. The resultant velocity is calculated from the two components 
at unknown mesh points. On the blade surface, the velocity is calculated 
from one component and the blade tangent angle. 

Computer Programs for Stream-Function Analyses 

As can be seen, the solution of Laplace^s equation and the calculation 
of velocities is a lengthy calculation procedure which is best done by com- 
puter. Several computer programs have been written at the NASA Lewds 
Research Center for the analysis of flow through turbomachine blading 
by stream-function methods. Most of these programs are for blade-to- 
blade analysis (region shown in fig. 5-9). The program called TURBLE, 
which is described in reference 5, can be used to analyze axial, radial, or 
mixed flow. In accordance with the constraints associated with the 
stream-function method, the flow must be subsonic throughout the entire 
solution region. The TSONIC program, described in reference 6, super- 
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sedes TURBLE in that it performs all the same calculations and, in 
addition, extends the solution to transonic (local supersonic velocities) 
flow problems. Transonic solutions are obtained by using a velocity- 
gradient equation of the type described in the next section to extend a 
preliminary (lower mass flow rate) subsonic stream-function solution. A 
program called TANDEM, which is described in reference 7, can be used 
to analyze flow in tandem or slotted blade rows or blade rows with 
splitters. Another program, called MAGNFY and described in reference 
8, obtains a detailed solution in the leading- or trailing-edge regions of any 
blade or in the slot region of tandem or slotted blades. The TANDE]\I 
and MAGNFY programs are restricted to subsonic flow. 

Flow in the meridional plane (mean hub-to-shroud flow surface, as 
indicated by fig. 4-14 (b) or fig. 5-1 (a)), of any axial- or mixed-flow 
turbomachine can be analyzed by a program called MERIDL, which is 
described in references 9 and 10. Transonic solutions can be obtained by 
the use of a velocity-gradient equation to extend a preliminary subsonic 
stream-function solution. 


VELOCITY-GRADIENT ANALYSIS 

As indicated previously, the stream-function and potential-function 
methods of analysis are limited to solutions that are entirely subsonic 
within the computation region. By use of a velocity-gradient equation and 
additional assumptions, however, the subsonic solution can be extended to 
give an approximate solution in the transonic flow regime. It is also 
possible to use a velocity-gradient method of analysis alone to obtain sub- 
sonic, transonic, or supersonic solutions without assumptions other than 
the basic ones indicated earlier. The velocity-gradient analysis is often 
called a stream-filament analysis because the velocity-gradient equation 
involves the streamline, or stream-filament, curvature and position. 

A velocity-gradient method of analysis can only give solutions within a 
guided passage; that is, a passage where both ends of all streamline 
orthogonal s intersect a solid boundary. Therefore, the usefulness of this 
method depends on the degree of flow guidance provided by the turbine 
blades. For a well-guided passage (high solidity and/or small angles), 
such as shown in figure 4-11, most of the suction surface is within the 
guided region, and the associated surface velocity distribution can be well 
defined. On the other hand, for a low-solidity blade row, such as that 
shown in figure 5-9, less than half of the suction surface is within the 
guided region, and surface velocities can be computed only on the front 
half of the suction surface. In this latter case, the stream-function analysis 
must be used if better definition of the suction-surface velocity distribution 
is required. 
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Method 

The idea of a velocity-gradient method can be demonstrated by con- 
sidering a simple case. Suppose we have two-dimensional flow through a 
narrow passage as shown in figure 5-13. We assume the height of the 
passage to be 5, and the width d. If the mass flow is known, the average 
velocity can be calculated approximately from continuity by 


W 


avQ 


W 

pbd 


(5-60) 


However, there is a variation in velocity across the width of the passage, 
and in turbomachinery it is this velocity difference we are interested in. 
With a force-equilibrium equation, by balancing centrifugal force against 
the pressure gradient as was done in chapter 3 for consideration of radial 
equilibrium, it can be shown that 


^_W 

dq Tc 


(6-61) 


where q is the distance from the suction (convex) surface, and Tc is the 
radius of curvature for the streamline. The sign convention for is 
important; is positive if it is concave upward, and negative if it is con- 
cave downward. For the simple case shown in figure 5-13, equation (5-61 ) 
can be integrated along a radial line by assuming the streamline radius of 
curvature to be equal in magnitude to the passage radius. There results, 
for integration from the inner radius to any point in the passage, 


W~r 


(5-62) 



Figure 5-13. — Flow through a curved passage. 
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where 


Ws relative velocity on inner, or suction, surface, m/sec; ft/sec 
r, radius of inner, or suction, surface, m; ft 
r radius of passage, m ; ft 

The mass flow through the passage is expressed as 


w = 



(5-63) 


and substitution of equation (5-62) into (5—63) and integration, with 
constant density assumed, yields 




w 


In 



(5-64) 


In a similar manner, the outer, or pressure, surface velocity can be com- 
puted as 


w 

= (5-65) 

A.(l+0l„(l+0 

Thus, an estimate of the blade-surface velocities can be obtained simply 
by using equation (5-62), which is a velocity-gradient equation. We are 
not necessarily restricted to two-dimensional flow. If there were some 
variation of velocity in the height of the passage, a velocity gradient could 
be calculated in that direction also. 

We will now consider a very general velocity-gradient equation. Since 
we are interested in turbomachinery, we will use a rotating cylindrical 
coordinate system with radius r, angle By and axis Xy as shown in figure 
5-14. Also indicated are the velocity components, Wr, Wzy and Wj. The 
meridional component is the resultant of Wr and W*. The meridional 
plane is a plane containing the x axis. Also shown in figure 5-14 are a, the 
angle between Wm and the x axis, and /S, the angle between W and the 
meridional plane. The following relations hold for the components: 


= sin 

(5-66) 

Wm-W cosfi 

(5-67) 

Wr—Wm sin a 

(5-68) 

Wx = Wm cos a 

(5-69) 
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Figuke 5-14.~Cylindrical coordinate system and velocity components. 


In addition to the r-, 0-, and x-coordinate, it is convenient to use an 
m-coordinate. The m-coordinate is the distance along a meridional stream- 
line, as shown in figure 5-15. The m-distance is less than the true stream- 
line distance if the angle The meridional streamline is the projection 
of a streamline in the meridional plane; that is, the ^-coordinate is 
neglected. The curvature of the meridional streamline is l/r^, where Vc is 
the radius of curvature of the meridional streamline. The sign of Tc is 
positive if the streamline is concave upward. 

We want the velocity gradient along an arbitrary curve. Let q be the 
distance along this curve. For the case of constant total temperature and 
constant angular momentum (rVu) at the inlet, 


where 


dq 


dr dx dd 


(5-70) 


W cos a cos^ ^ W sin^ p . dWm 

a — hsin a cos — 2oi sm /3 (5-71) 

Tc r dm 
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Figure 5-15. — The w-coordinate. 


W sin a cos^ /S , ^ dW„ 

h= — -+COS a cos p — 

dm 


(5-72) 


c=W sin a sin /3 cos 0+r cos 0 


/dW, „ . \ 

( l-2w sin a ) 

\dm / 


(5-73) 


These equations are derived as equations (B13) and (B14) of reference 11. 
In using any velocity-gradient equation, it is necessary to solve a differ- 
ential equation involving streamline-geometry parameters, such as cur- 
vature, a, and These are not known precisely in advance. However, 
for a guided channel, these parameters can be estimated reasonably well. 

A great number of special cases can be obtained from equations (5-70) 
to (5-73). For example, suppose we have an annular passage with no 
blades, as shown in figure 5-16, and no velocity component in the tangen- 
tial (fl) direction (into page). We can calculate dW/dn, where n is the 
distance normal to the streamline. Let g = n in equation (&-70). Since 
Wt = 0, then dd/dn = 0 and (3 = 0. Further, from figure 5-16, it can be 



Figure 5 - 16 .— Annular passage with no blades. 
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seen^at dr/dn = cos a and dx/dn^ -sin a. Then, from equations (5-70) 


dW W 
dn ~ re 


Thus, for this case, equation (5-70) reduces to 
tion (5-61). 


the simple form of equa- 


v^mputer Jrrograms 

Several computer programs for the analysis of flow through turbo- 
machine blading by velocity-gradient methods have been written at the 
Research Center. One that was used for many years is the 
C i TD program which is described in reference 12 and is limited to axial- 
ow turbines. This program has now been superseded by the more general 
and to use CHANEL program, which is described in reference 13. 
Ihe CHANEL program can be used to analyze axial-, radial-, or mixed- 
flow turbines or compressors. Velocity-gradient equations are used to 
e ermine velocity variations both from hub to tip along meridional- 
streamline orthogonals and from blade to blade along hub-, mean-, and 
ij^streamhne orthogonals. This results in a flow solution for an orthogonal 
surface, as illustrated in figure 5-17, which satisfies a specified mass flow 
rate. Computations are made for a number of these surfaces along the 
blade passage. This program can also be used to compute the maximum 
(choking) mass flow rate for the channel. The program gives good results 
for medium- to high-solidity blading. As indicated previously, more 
definition than can be provided by this program may be needed for low- 
solidity blading, because solutions can only be obtained for fully guided 
sections of the passage. 

Velocity-gradient methods have also been used to obtain meridional- 
plane and blade-to-blade plane solutions. The basic method for a meridi- 
onal-plane analysis for mixed-flow centrifugal impellers is presented in 
reference 14, which uses the velocity-gradient equation along streamline 
orthogonals. Since the orthogonal lengths are not known in advance it 
was more convenient to base a computer program on the use of the 
ve ocity-gradient equation along fixed straight lines, which were called 
quasi-orthogonals. Such programs for meridional-plane analysis are 
presented in reference 11 for a radial-inflow turbine impeller and in 
reference 15 for backward-swept or radial impellers and vaned diffusers of 
centrifugal compressors. A program for a blade-to-blade plane analysis 

that uses quasi-orthogonals for a radial-inflow turbine impeller is described 
in reference 16. 

A further use of the velocity-gradient equation, as mentioned pre- 
viously in this chapter, is to extend a subsonic stream-function solution to 
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obtain local supersonic velocities. The subsonic solution is used to obtain 
the flow angles and streamline curvatures required for the velocity- 
gradient equation. Programs for transonic-flow solutions based on this 
method are presented in references 9 and 10 for a meridional solution and 
in reference 6 for a bladc-to-blade solution. 



Figure 5-17. — Turbine blades with three-dimensional orthogonal surface across 

flow passage. 
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SYMBOLS 

A flow area, m^; 

ai coefficients for equation (5-46) 

b cascade height, m; ft 

d passage width, m ; f t 

h distance between mesh points, m; ft 

ki constant in equation (5-56) 

m distance along meridional streamline, m ; ft 

n distance normal to streamline, m ; ft 

p absolute pressure, N/m^; \b/W 

q distance along arbitrary curve, m; ft 

R gas constant, J/(kg) (K) ; (ft) (lbf)/(lbm) (°R) 

r radius, m; ft 

s blade spacing, m; ft 

T absolute temperature, K; °R 

t time, sec 

u stream function 

V absolute velocity, m/sec; ft/sec 

W relative velocity, m/sec; ft/sec 

w mass flow rate, kg/sec ; Ib/sec 

a fluid absolute angle of inclination from axial direction in the 

meridional plane, deg 

/3 fluid flow angle, relative to blades, out of the meridional plane 

(in the tangential direction), deg 

7) distance in direction of outer normal to cascade boundary, m; ft 

B angular distance in direction of rotation, rad 

P density, kg/m^; Ib/ft^ 

potential function 
f angular velocity, rad/sec 
\overrelaxation factor 

Subscripts: 

c curvature 

in inlet 

m meridional component 

n component normal to streamline 

out outlet 

p pressure surface 

r radial component 

8 suction surface 

X axial component 

y component in y-direction 

z component in z-direction 
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B tangential component 

0 I 

1,2,> mesh-point designations 

3,4j 
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CHAmR 6 


Introduction to Boundary- 
layer Theory 

By William D. McNally 


As shown in chapter 2 , the pressure ratio across a turbine provides a 
certain amount of ideal energy that is available to the turbine for pro- 
ducing work. The portion of the ideal energy that is not converted to 
work is considered to be a loss. One of the more important and difficult 
aspects of turbine design is the prediction of the losses. 

Before losses can be predicted, it is necessary to understand their 
causes. The primary cause of losses is the boundary layer that develops 
on the blade and end-wall surfaces. Other losses occur because of shocks, 
tip-clearance flows, disk friction (windage), flow incidence, and partial- 
admission operation. This chapter gives an introduction to boundary- 
layer theory, which is used to calculate the parameters needed to estimate 
viscous (friction) losses. ^Methods for determining the basic viscous loss 
and the associated trailing-edge and mixing losses are presented in the 
next chapter. 


NATURE OF BOUNDARY LAYER 

When a real fluid (such as air) flows past a turbine blade at normal 
velocities, the influence of viscosity on the flow is confined to a relatively 
thin layer in the immediate neighborhood of the blade. This layer is called 
the boundary layer. At the outer edge of this layer the flow is frictionless, 
and conditions there agree with those calculated with the use of ideal 
(frictionloss, nonviscous) flow assumptions. At the wall, on the other 
hand, the velocity of the fluid is zero in all directions (no-slip condition). 
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It is the frictional, or viscous, forces in this thin layer that reduce the fluid 
velocity from its free-stream, frictionless value to zero at the wall. 

A boundary layer on a turbine blade is illustrated in figure 6-1. The 
boundary layer develops from a small finite thickness at the stagnation 
point at the leading edge of the blade and grows along both the suction 
and pressure surfaces. The initial portion of the boundary layer is always 
laminar. In a laminar boundary layer, fluid layers parallel to the blade 
surface slide over each other. Any minute local fluctuations in velocity are 
sufficiently damped so that they have negligible influence on the smooth- 
ness of the overall flow. The velocity at a point is either steady with time 
or changes in some smooth way, as figure 6~2(a) indicates. 

]\Iost flows being ducted to a turbine, or entering it from a combustor, 
are turbulent in nature. The fluctuating components of velocity have a 
significant influence in this type of flow. With this overall flow, the 
boundary layer on the blades cannot remain laminar for any great dis- 
tance. It usually passes through a transition region and becomes a tur- 
bulent boundary layer. In the transition region, weak disturbances in the 
flow are amplified, and this leads to the random fluctuations in velocity 
that are characteristic of turbulent flow. In the turbulent boundary layer, 
as in turbulent flow, the velocity at any point oscillates in a random 
fashion about a mean value, as figure 6-2 (b) indicates. 

Figure 6-1 also shows a separated region in the turbulent boundary 
layer. Separation can likewise occur in the laminar boundary layer. When 
a boundary layer separates, the fluid moves away from the blade surface. 
The manner in which this happens is illustrated in figure 6-3. As the free- 
stream velocity decreases along the rear portion of the suction surface of 
a turbine blade, static pressure correspondingly increases. This positive 



Figure 6-1. — Boundary layer on blade- 
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Steady 


Steady 


Unsteady 



(a) 


t 


Unsteady 




(a) Laminar flow. 

Figure 6-2. — Variation of 


(b) Turbulent flow, 
velocity with time at a point. 



Figure 6-3.— Boundary-layer separation. 


pressure gradient (adverse pressure gradient) retards the flow in the 
boundary layer and causes it to lose energy. The flow in the boundary 
layi'r can be retarded to such a degree that very close to the wall it moves 
in a direction opposite to that of the mean flow passing the blade. This is 
separation. The point at which the flow reverses itself is the separation 
point. The laminar boundary layer at the leading edge of a turbine blade 
can also separate and immediately reattach itself to the surface as a 
turbulent boundary layer. This is illustrated in figure 6—4. 

Finally, it should bo noted that both laminar and turbulent boundary 
layers can be either incompressible or compressible, depending on the 
level of the Mach number. Just as there are diflferent equations to repre- 
sent laminar and turbulent boundary-layer flow, there are different 
equations for the incompressible and compressible variations of each. 
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Turbulent 

boundary 



Figure 6-4. — Laminar separation and reattachment. 


Boundary layers should be considered compressible if the free-stream 
relative Mach number exceeds values of 0.3 to 0.4. The boundary-layer 
equations for these various cases are derived and solution methods are 
discussed in this chapter. 

DERIVATION OF BOUNDARY-LAYER EQUATIONS 

The general equations of motion of viscous fluids are called the Navier- 
Stokes equations. In normal coordinate systems, there are three such 
equations, one for each of the coordinate directions. The boundary-layer 
equations can be derived from the Navier-Stokes equations. The Navier- 
Stokes equations themselves can be derived by applying the law of con- 
servation of momentum to a fluid element. This exercise is length}^, and 
will not be repeated here. References 1 and 2 both have the complete 
derivation, in two somewhat different forms. 

There are various forms of the Navier-Stokes equations, depending on 
what assumptions are made during their derivation. The following equa- 
tion represents the Navier-Stokes equations combined into vector form 
for a compressible fluid with constant viscosity: 

a a 1 u 

-T = g^-~ Vp+- V^uT-- V(V-u) (6-1) 

at p p 6 p 

where 

u general velocity vector, m/sec; ft/sec 
t time, sec 

g conversion constant, 1; 32.17 (Ibm) (ft)/(lbf) (sec^) 

f general body force acting on a unit mass of fluid, N/kg; Ibf/lbm 
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p density, kg/m*; Ibm/ft* 

p static pressure, N/m^; Ibf/ft^ 

/i dynamic viscosity, (N) (sec)/m*; lbm/(ft) (sec) 

In this equation, u represents a general velocity vector with components 
u, V, and ie in the three coordinate directions x, y, and 2 , respectively. 

u = wi+yj+iuk (b-2) 

where i, j, and k are the unit vectors in the three coordinate directions. 
The total, or substantial, derivative of u is du/dL In any of the coordinate 
directions, 


d 

dt 


d 

dt 


d d 


+ U— + V— + W 

dx dy 


dz 


(6-3) 


In equation (6-1), the Laplacian operator is applied to the vector u 
rather than to a scalar function. If the term V^u is expanded into simple 
vector quantities, equation (6-1) becomes 

^ = vp + - V(V-u)--CVX(VXu)]+i- V(V-u) (6^) 

dt p P P dp 

Expressing the V operator in terms of gradients, curls, and divergences, 
which may be more familiar to the reader, equation (6-4) becomes 


— = ^f— - grad p+- grad(div u) — - curl(curl u) - grad(div u) 
dt p p P dp 

(6-5) 


In order to derive the boundary-layer equations, equation (6-1) has 
to be expanded into three scalar equations, one for each of the coordinate 
directions. The three resulting equations are 


du du du du g dp u / d^U d^U d^u\ 

dt dx dy dz p dx p \dx^ dy^ dz^ / 


a. M 

3 p dx \dx dydzj 

dv dv dv dv . 9 f 

dt~^^ dx~^^ dy^^ dz p dy p di/^ dzy 

(fr-7) 

3pdy\dx dy dz) 
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dw dw 

\-u (- 

dt dx 


dw dw ^ ^ M /' 

dy dz pdz p\ 


d’^w dhv d^w\ 


Ifld /du 

+o “ I “+ 


3 p dz \dx 

where /r, fj,^ and ft are the components of the body force f. 


dv die\ 
dy dz ) 


( 6 - 8 ) 


Laminar Incompressible Boundary Layer 

In order to derive Prandtl's boundary-layer equations for laminar 
incompressible flow, the following assumptions will be made: 

(1) Viscosity is a constant. This has already bec^n assumed in the 
writing of the preceding equations, 

(2) Flow is incompressible. Since for incompressible flow the con- 
tinuity equation is 


V* u = div 


/ du dv div\ 

u = {— +— H 1 = 0 

\dx dy dz / 


(6-9) 


the final terms in equations (G-6) to (6-8) can be eliminated. 

(3) Flow is two-dimensional. This eliminates equation (6-8) from 
consideration, as well as all terms involving w or d/dz in equations (6-6) 
and (6—7) . 

(4) Flow is steady. This eliminates d/dt terms, 

(5) Body forces are negligible in relation to inertia and viscous forces. 
Thus,/x and/i/ can be discarded from equations (6-6) and (6-7). 

With these assumptions, the Navier-Stokes equations reduce to the 
following two equations for the x- and ^/-directions ; 


du du gdp u 1 

dx dy p dx p * 

^ d^a d^u\ 
\dx^^dy'‘) 

(6-10) 

dv dv gdp Ilf 

dx dy pdy p" 

^dhr dh\ 
\dx^^dy‘‘) 

(fr-H) 

Likewise, the continuity equation becomes 



du dv 

— H = 0 

dx dy 


(6-12) 


In order to make equations (6-10) to (6-12) suitable for the analysis of 
boundary-layer flow, the equations arc traditionally made dimensionless, 
and an ordcr-of-magnitude check is performed on the various terms to 
show that some are negligible with respect to others. Figure 6-5 shows the 
velocities and coordinate directions pertinent to the boundary layer. 
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The following dimensionless parameters are defined : 


II 

(6-13a) 

II 

(6-13b) 

II 

(6-13c) 

II 

(6-13d) 

p_ gp 

pW 

(&-13e) 

^ pLU, 
Re — 

M 

(6-13f) 


where 

X dimensionless x-coordinate 

L characteristic length (in this case, the blade chord) , m; ft 

Y dimensionless ^-coordinate 

U dimensionless velocity in x-direction 

Uq free-stream velocity upstream of blade, m/ sec ; ft/ see 

V dimensionless velocity in y-direction 
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P dimensionless pressure 

Re Reynolds number 

From figure 6-5, we see that since x is proportional to L, X is of order 1. 
And since y is proportional to the boundary-layer thickness bfuii, Y is of 
order 5/un/^ = e, a quantity much less than 1. likewise, since u is of order 
Uoy U = u/Uo is of order 1. And V — v/Uo is of order e, since velocities in 
the y-dircction in the boundary layer are much smaller than those in the 
x-direction. 

In order to put equations (6-10) to (6-12) in terms of dimensionless 
quantities, equations (6-10) and (6-11) are multiplied by L/W, and 
equation (6—12) is multiplied by L/f/o. The resulting dimensionless 
equations are 


U 


u 


dU 

dX~ 

dV 

dx' 


dU dP 1 

/d^U dW\ 


dY^~dX'^7{e 


(6-14) 

dV dP 1 , 





(6-15) 

dU dV 

= 0 

dX dY 


(6-16) 

the various terms in these equations 

can now 


be compared with each other. Since X and U are of order 1, and Y and V 
are of order c, 


at/ 1 
ax~i~^ 

(6-17a) 

dU_l 

dY~t 

(6-17b) 

dV_e 

ax~i~* 

(6-17c) 

dV t 

(6-17d) 

1 

aX2“i-i~^ 

(6-17e) 

a*f/_ 1 1 

dY^ € • e 

(6-17f) 


64 



INTRODUCTION TO BOUNDARY-LAYER THEORY 


dW 


(&-17g) 


dW € 1 

dY^~f€~( 


(6-17h) 


Furthermore, the change in P with respect to X is of the same order of 
magnitude as the change of U with respect to X, so that dP/dX is of 
order 1. 

Relating these orders of magnitude to the terms in equations (6-14) 
to (6-16) yields 


^dU dP 1 /d^U dW\ 

U f-F — = 1 ( 1 I 

ex dY dX Re\dX^ dYy 

(l)(l)+(t) 0=-l+(t^) 

^ dX~^^ dY~ ~ dY~^ Re 

(l)(*) + (*)(l) = -«+(*^) (t+J 

dU dV 

I = 0 

dX dY 


(6-18) 


(6-19) 


( 6 - 20 ) 


1 + 1 

By examining equations (6-18) to (6-20), the following conclusions 
can be reached : 


(1) In boundary-layer theory, it is assumed that the viscous terms 
l/Re\i(d^U/dX^) + (d^U/dY^)'] are of the same order of magnitude as the 
inertia terms U{dU/dX) -\-V{dU/dY). For this to be true in equation 
(6-18), \/Re must be of order c^, since d‘^UfdY^ is much larger than 
d^UjdX^ and dominates the two terms in parentheses. Therefore, the 
Reynolds number must be relatively large. 

(2) In equation (6-19), with 1/Re of order and with dW/dY^ 
dominating d'^V/dX^j the terms are of order e. Therefore, unless dP/dY 
is to dominate, it too must be of order e or less. Therefore, dP/dY is much 
smaller than dP/dX^ and P can be considered a function of X alone. There- 
fore, P = P{X) or p = p(x), and dP/dX — dP/dX or dp/dx = dp/dx. This 
allows us to assume that the pressure across the boundary layer in the 
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^-direction is essentially constant. It can be assumed equal to the potential 
flow pressun* existing at the outside of the boundary layer. 

(3) Since the first equation is of order 1, and the second equation is of 
order e, the second equation can be neglected completely. 

(4) In equation (6-18), dW/dX^ can be neglected because it is so 
small in comparison with dWIdY'^. This leaves the following dimension- 
less equations : 


^ dY'^ dx'^RedY'^ 


(^ 21 ) 


dX'^dY 


( 6 - 22 ) 


These are PrandtFs boundary-layer equations in dimensionless form. 

The boundary-layer equations in this form are useful in determining 
the influence of the Reynolds number on the size of the boundary layer 
for different fluids. From e(juation (6-21) we see that as Re increases in 
magnitude, the viscous-force terms {I /Re) (d^U/dY'^) will get smaller and 
smaller. The boundary-layer thickness will correspondingly decrease. So, 
as Re increases, d/uii decreases. Furthermore, increasing Re corresponds to 
decreasing viscosity if pLUo is constant. So, as a general rule, the thickness 
of the boundary layer decreases as the viscosity decreases. 

The boundary-layer equations can be put in terms of dimensional 
variables by multiplying equation (6-21) by U(//L and equation (6-22) 
by Uq/L. The resulting equations are 


du du g dp d^u 

dx dy p dx p dy"^ 


(6-23) 


du dv ^ 

I = 0 

dx dy 


(6-24) 


These are Prandtl’s boundary-layer equations for two-dimensional, 
laminar, incompressible flow. Density and viscosity are assumed constant 
and known. The pressure gradient along the blade surface, dp/dx, is also 
known from an ideal-flow solution. The remaining unknowns are u and r, 
and equations (6-23) and (6-24) arc sufficient for their calculation. 

It should be noted that the boundary-layer equations are not valid in 
the presence of shock waves (i.c., where instantaneous adverse pressure 
gradients of large magnitude occur). Just as flow phenomena in the 
boundary layer depend on mainly the Reynolds number, conditions in a 
shock wave depend on primarily the Mach number. Since the influence of 
Mach number is not included in the boundary-layer et^uations, they tell 
us nothing about the interaction of shock waves and boundary layers. 
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The boundary-layer equations are not completely reliable as separation 
is approached. One of the assumptions used in their derivation is that the 
velocity v is much smaller than w. Very close to the separation point, the 
boundary layer grows rapidly, and v begins to be of the same order as u. 
Nonetheless, the boundary-layer equations are gen(*rally used in cal- 
culations right up to the separation point, since the region where V is 
significant is very small, and little error in the location of the point of 
separation is incurred. However, these equations should not be used for 
detailed calculations in the neighborhood of a separated flow region. 

The Navier-Stokes equations (6-6) to (6-8) used in the development 
of the boundary-layer equations were derived for an orthogonal system 
of coordinates in which the radius of curvature of each of the coordinate 
axes is quite large (i.e., where curvature effects are negligible). The 
question arises as to how the boundary-layer equations would change for 
flow over a curved wall. If a curvilinear orthogonal coordinates system 
(fig. 6-6) is introduced wherein the x-axis is in the direction of the curved 
wall and the ?/-axis is normal to it, a new set of Navier-Stokes equations 
can be derived for flow in such a system. These eejuations are given in 
reference 1. The terms in the equations arc very dep(‘ndent on the radius 
of curvature r at a position x along the blade surface. The relative orders 
of magnitude of the individual terms can be estimated in the same manner 
as was done previously. With the assumption that the boundary-layer 
thickness is small compared with the radius of curvature of the wall, and 
for the case where no largo variations in curvature occur, so that dr/dx ^ 1 , 
the same boundary-layer equations result as were obtained for flat walls. 
Therefore, the flat-plate boundary-layer equations may be applied to 
curved walls as well, provided there are no large variations in curvature, 
such as would occur near sharp edges. 



Figure 6-6. — Curvilinear coordinate system on a blade. 
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Laminar Compressible Boundary Layer 


An order-of-magnitude analysis can also be performed to derive the 
equations for a compressible boundary layer. In the incompressible case, 
viscosity and density were assumed constant, temperature variations were 
neglected, and the energy equation was not used. For the compressible 
case, density is no longer constant, viscosity is considered a function of 
temperature, the equation of state is used to relate pressure and density to 
temperature, and, if the process is not isothermal, some form of the energy 
equation is required. The boundary-layer equations for compressible, 
nonisothermal, variable-viscosity flow will involve three parameters which 
can be related to temperature. These are viscosity, specific heat, and 
thermal conductivit 3 \ 

There are sev(‘ral relations for viscosity as a function of temperature. 
The most common is probablj' Sutherland’s relation (ref. 1) 


M _/TY^To+S 
ho~\tJ T+S 


(6-25) 


where 

tiQ dynamic viscosity at the reference temperature To, (N)(scc)/m2; 
Ibm/ (ft) (sec) 

T absolute static temperature, K; 

To reference temperature, K; 

S a constant, K ; °R (for air, = 1 10 K or 198° R) 

A less complicated, but also less accurate, temperature-viscosity relation 
is the power la^^' 


~ = (^) 0.5<co<1.0 (6-26) 

MO \To/ 

where oj is a constant. For air, w is approximately 0,65. 

Specific heat and thermal conductivity can be relatc^d to temperature 
by least-squares poljuiomial-curve fits for the particular gas and tem- 
perature range involved. With these variables related to temperature, 
the unknowns in the compressible-boundary-layer problem reduce to u, 
r, p, and T. The four equations relating th('S(^ variables will be the con- 
tinuity equation, one component of the momentum equation, the state 
equation, and the energy equation. 

The order-of-magnitude analysis of the contiriuit}^ and momentum 
(Navier-Stokes) equations for comf)ressible flow is almost identical to 
that for incompressible flow. For compressible flow with nonconstant 
viscosity, the equations analogous to (6-10) to (6-12) are the following: 
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du du 

pu \-pV T“ 

dx dy 


dp a r„ au 2 /au aw\] a r 


(6-27) 


dv dv 

pu - — \-pV 

dx dy 


dp a 


V dv 2 /a« ,a«\] a r /dv duV 


(6-28) 


djpu) a(pti) ^^ (6-29) 

dx dy 


If an order-of-magnitude analysis is performed on these equations similar 
to that for the incompressible-flow equations, the following boundary- 

layer equations result : 


du du dp d / du\ 

(6-30) 

d{pu) d(piQ_^ 
dx dy 

(6-31) 


The equation of state is also required for the solution of compressible 
boundary-layer flow. The state equation is 

p = pRT (^^ 2 ) 

where R is the gas constant, in J/(kg) (K) or (ft) (lbf)/(lbm) (°R). 

The final equation required besides the continuity equation, the 
momentum equation, and the equation of state, is the energy equation. 
The energy equation for a compressible boundary layer is derived from 
the energy equation for a perfect gas by means of another order-of- 
magnitude check. The following is the energy equation for compressible, 
two-dimensional steady flow of a perfect gas, written in full . 


/ dT V-v> (6-33) 

^^Vdx'^^dy) J dx^J dy^dx\ dx^dyK dy ) gJ 


where 

Cp specific heat at constant pressure, J/(kg) (K) ; Btu/ (Ibm) (°R) 

J conversion constant, 1; 778 (ft) (lbf)/Btu ^ 

k thermal conductivity, W/ (m) (K) ; Btu/ (sec) (ft) ( R) 


and 


(6-34) 


r/auV /dvy] , duV 2/du dvV 

w +y 
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If an order-of-magnitude check is performed on the above equations the 
following boundary-layer energy equation results : 


pOp 


^ dT dT\ 
w T — f-r — 1 
^ dx dy / 


1 * (i 

Jdx^dy V SJVW 


(6-35) 


Equations (6-30), (6-31), (6-32), and (6-35) are the laminar boundary- 
layer equations for nonisothermal, two-dimensional, compressible flow of 
a gas obeying the ideal gas law. 


Turbulent Boundary-Layer Solution Methods 

It is desirable to have a turbulent boundary layer over the major portion 
of turbine blades. If the boundary layer is not turbulent, separation will 
probably occur on the blades, with a resulting decrease in their per- 
formance. Turbulent flow has irregular fluctuations (mixing or eddy 
motions) superimposed on the main fluid motion (see fig. 6-2). These 
fluctuations are so complex that closed-form solutions are not feasible at 
present Yet, the mixing motion is very important, since the stresses in 
the fluid due to fluctuating components of velocity arc often of greater 
magnitude than those due to the mean motion. 

There are two approaches to the solution of turbulent boundary-layer 
flow . The first is the exact solution of the time-dependent, three-dimen- 
sional, Xavier-Stokes equations. The three-dimensional equations are 
required, since two-dimen.sional calculations could never represent the 
stretching of eddies, which is a prime mechanism of turbulent flow'. How- 
evc'r, even the largi'st computers available at the present time cannot 
handle three-dimensional solutions of these equations on a small enough 

mesh to represent the fluctuating components of velocity of turbulent 
flow. 

The second approach is to write the equations of continuity, momentum, 
and energy in terms of mean and fluctuating components of pressure] 
density, temperature, and velocity. In this approach, the time average 
of the u component of velocity, for example, is denoted by u and the 
velocity of fluctuation by u'. So the velocities, den.sity, pressure, and 

temperature are written as follows: 

u = u+u' (6-36a) 

v = v+v' (6-36b) 

P = P+P (6-36c) 

P = P+P' (6-36d) 

T^f^-r ( 6 - 36 c ) 

The fluctuations in viscosity, thermal conductivity, and specific heat are 
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negligible and are not considered. So these three parameters are calculated 
as functions of the time-averaged value of temperature. 

If the flow properties listed in equations (6-36) are substituted into the 
continuity, momentum, and energy equations for incompressible and 
compressible flow, a new sot of stress terms arises in thees equations. 
These are called the “apparent” turbulent stresses, or Reynolds stresses. 
They are of the form pu^ and pi??, where uV is the average over tirne of 
the product of u’ and v'. These new terms in the equations add additional 
unknowns to the boundary-layer problem for which additional equations 
arc not presently available. For this reason, empirical expressions or 
approximations are substituted for the Reynolds stress terms before the 
turbulent boundary-layer equations can be solved. 


Turbulent Incompressible Boundary Layer 


Substituting the relations of equations (6-36) into equations (6—10) , 
(6-11), and (6-12), and then performing an order-of-magnitude analysis 
yields the following equations for turbulent, incompressible, boundary- 
layer flow : 


du dU dp d / du 

pu f-p5— =— + I 

dx dy dx dy\dy / 


(6-37) 


dx dy 


(6-38) 


These equations are analogous to equations (6-23) and (6-24) for laminar 
flow. Notice, however, the presence of the Reynolds stress term in the 
momentum equation. This adds a new unknown {u'v') to the original two 
(m and v), thereby making three unknowns with only two equations. 


Turbulent Compressible Boundary Layer 

Substituting the relations of equations (6-36) into equations (6-27), 
(6-28), (6-29), (6-32), and (6-33) and then performing an order-of- 
magnitude analysis yields the following equations for turbulent, com- 
pressible, boundary-layer flow : 

dipu) djjiv) d(p V ) p (6-39) 

dx dy dy 


du du , ^du dp d / du 
pu — \-pv—+pv — = -g-j-+—\p - — ) 
dx dy dy dx dy\ dy / 


p=prT 


(6-40) 

(6-41) 
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. f.df, df\ , —dTt Udp ^ d /, dT,\ 

pCp I u f-p — ]-\-Cpp V — = — I k — 1 

dx dyr dy Jdx^dyX dy) 


+ — — 
gj dy 


where T i is the absolute total temperature, in K or °li, and is defined as 


T, = T 


-I 

2gJcp 


(6-43) 


Wc have now derived the basic boundary-layer equations for two- 
dimensionalj laminar and turbulent, incompressible and compressible 
boundary-layer flow. Wc should note at this time that this is really only 
the starting point as far as boundary-layer solutions are concerned. 
These equations are only the basis for the many, many methods which 
presently exist for obtaining boundary-layer solutions under various 
circumstances. 


SOLUTION OF BOUNDARY-LAYER EQUATIONS 

After velocity profiles are discussed and the important boundary-layer 
parameters defined, some of the solution methods will bo discussed. 
Included will be the flat-plate, incompressible solution, as well as com- 
pressible methods. 

Velocity Profiles 

One of the principal results obtained from most boundary-layer solu- 
tions is a description of the velocity profile in the boundary layer along 
the blade surface (fig. 6—7). The velocity profile describes mathematically 



Figure 6-7. — Boundary -layer velocity profiles. 
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Outer edge of 
boundary layer 

y 





(b) 

^777777777777 


(a) Laminar 
profile. 

Figuke 6-8.— Laminar and 


(b) Turl)ulent 
profile. 

turbulent velocity profiles. 


the dimonsionless velocity w/ue as a function of the dimensionless distance 
yl^juii from the blade. The velocity u is the velocity in the boundary layer 
at a distance y from the surface, and the velocity \it is the external free- 
stream velocity at a distance equal to the boundary-layer thickness, 
5/uu, from the surface. Alternately, bfuii is often defined as that distance 
from the blade where the velocity differs by 1 percent from the external 
velocity, w,. 

Velocity profiles for laminar flow (fig. 6-8(a)) tend to be parabolic in 
shape, while those for turbulent flow are blunted (fig. 6-8(b)). A com- 
monly used mathematical expression for u/u^ in laminar flow is that 
originated by Pohlhausen (sec ref. 1) : 


— +6 Y+c Y+^ Y 

Ue bfull \bfuU/ \bfuU/ 


(6--i4) 


The constants a, 6, c, and d are defined in terms of a dimensionless shape 


parameter 

II 

(6-45) 

where 




a — 2+- 


(6— 46a) 


6 - 


X 

2 


(6-46b) 


c = 



(6-46c) 
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Figuue 6-9, — Laminar velocity profiles. 


rf=l-- (6^6d) 

6 

Velocity profiles for various values of X are shown in figure 6-9. 

Velocity profiles for turbulent flow are often represented by the power 
law 



Pipe-flow experiments show that the exponent n is a mild function of the 
Reynolds numb(T and varies from 4 up to about 10. The value of n = 7 
is most appropriate for boundary-layer flow on a flat plate. The exponent 
n can be related to other boundary-layer parameters, namely the dis- 
placement thickness 5 and the momentum thickiu'ss 6, which are described 
in the next section. 

Definitions of Important Boundary-Layer Parameters 

Solutions of the two-dimensional boundary-layer equations are most 
often obtained in terms of three important param(4ers. These are the dis- 
placement thickness 5, the momentum thickness 6, and the form factor H . 
In order to define these parameters, it is necessary to first define the 
thickness of the boundary lay^T, d/uii- The definition of boundary-layer 
thickness is rather arbitrary, since transition from the velocity inside the 
boundary layer to that outside it takes place asymptotically. This is of 
little importance, however, because the velocity in the boundary layer 
attains a value which is very close to the external v(4ocity at a small 
distance from the wall. It is possible to define the boundary-layer thick- 
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(a) Actual velocity profile. (b) Equivalent profile for 

equal mass flow. 

Figure 6-10.— Displacement thickness of a boundary layer. 


ness as that distance from the blade where the velocity differs by 1 percent 
from the external velocity Ue. 

The displacement thickness 5, for compressible boundary-layer flow, 
can be defined with the help of figure 6-10. As seen from figure 6-10 (a), 
the decrease in mass flow within the boundary layer due to the influence of 
friction is given by 


/ V = bfult 

{peUe—pu) dy (6-48) 

where p« is the density, in kg/m* or Ibm/ft*, in the free stream outside of 
the boundary layer. This integrated mass defect can be represented by a 
distance 5, the displacement thickness, as shown in figure 6-10(b). It is 
the distance by which the external potential field of flow is displaced 
outward as a consequence of the decrease in velocity in the boundary 
layer. 

As figure 6-10 shows, the distance b can be defined by the equation 

/ V = ^full 

ipeiie-pu) dy (&-49) 

V-0 

Solving for b gives 


'— /■ 
P.M. •'v- 


S= / (peU, — pu)dy=l 

V-»0 \ Pe^c/ 


dy 


The displacement thickness for incompressible flow reduces to 


(6-50) 


-a 


y=^Svn 


{Ue~u) dy 


rU-S,ul, / X 

= / (^-r)dy 


e *' 


(6-51) 


The loss of momentum in the boundary layer due to the presence of 
friction is given by 
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/ y=^fuii 

pu{ue—u) dy (6-52) 

j/=0 


This momentum defect from the momc^iitum of purt'ly potential flow 
can be represented by a distance 0, defined by the ecpiation 

pu{ue—u) dy (6-53) 

T/-0 


Solving for d in this equation giv(‘s the dc’finitiori of the momontum thick- 
ness for compr(\ssiblo boundary layers as 


Y rV=^full 


pu{jle—u) dy = j 


y = ifuii 



(G-54) 


The momentum thickness for incompressible flow reducc-s to 


0 = — [ u{u,-u) dy-= I — ( 1 ) dy (6-55) 

■'y_0 •^y-O V 


The form factor H for both compressible and incompressible flow is 
defined as the ratio of displacement thickness to momentum thickness: 



(6-56) 


There are many other boundary-layer parameters besides 5, d, and H for 
two-dimensional, and especially for three-dimensional, boundary layers. 
These three, however, are the principal parameters used in general 
boundary-layer studies. 

Physical Interpretation of Separation 

When separation of flow from a blade or a casing occurs, some of the 
retarded fluid in the boundary layer is transported away from the surface 
toward the main stream. When a region with an adverse pressure gradient 
exists along a surface, the retarded fluid particles cannot, in general, 
penetrate too far into the region of increased pressure because of their 
small kinetic energy. Thus, the boundary layer is deflected away from the 
surface and moves into the main stream. In general, the fluid particles 
behind the point of separation follow the pressure gradient and move in a 
direction opposite to the external stream. The point of separation is 
defined as the limit between forward and reverse flow in the layer in the 
immediate neighborhood of the wall. At separation, 

(-) =0 (6-57) 

\9y/y-a 

Figure 6-1 1 illustrates separation occurring along a surface. 
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By oxamining Prandtrs boundary-layer equations and considering the 
relation between pressure gradient dp/dx and velocity distribution u(y)^ 
it is possible to infer that separation in a steady flow will occur only in 
the presence of an adverse prc'ssure gradient (i.e., decelerated flow), 
dp/dx>Q. From equation (C-23), with the boundary conditions at the 
surface being ?/ = = 0, we have 



We can now relate velocity profiles to du/di/j d'^ujdy^y and finally to 
dp/dx through equation (0-58). The eejuation indicates that in the 
immediate neighborhood of the wall, the curvature of the velocity profile, 
d^u/dy’^j deptmds only on the pressure gradient, dp/dxj and the curvature 
of th(! velocity profile at th(^ wall chang(*s its sign with the pressure 
gradient. 

figure 6-1 2 (a) shows a velocity profile' that would exist in a boundary 
layer subjected to a decreasing pressure. For such a profile, figure 6-12 (b) 



(a) Velocity (b) Velocity (c) Velocity- 

profile. gradient. profile 

curvature. 

Figuke 6-12.— Velocity distribution in a boundary layer with pressure decrease. 
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indicates that bujdy is positive for all y and decreases as y increases. 
Furthermore, figure 6-12(c) indicates that d^^ufdy^, which is the slope of 
dufdy, is negative for all y. From equation (6-58) , we know that negative 
d^ufdy^ corresponds to negative dpidx. Consequently, a boundary layer 
subjected to a decreasing pressure (negative dpfdx) will have velocity 
profiles which are not indicative of impending separation (the form of 
fig. 6-12(a)). 

Figure 6-13 (a) shows a profile which would exist in a boundary layer 
with decelerated flow due to an increasing pressure (adverse pressure 
gradient). Here, figure 6-13(b) indicates that dujdy has a positive slope 
near the blade surface; that is, d^u/dy^ is positive (fig. 6-13(c)). This 
corresponds to positive dp/dx. However, since in all cases d^u/dy^ must 
be less than zero at some distance from the surface, there must exist a 
point for which d^u/dy^ = 0. This is a point of inflection of the boundary- 
layer velocity profile. It follows that in a region of retarded potential flow 



(a) Velocity (b) Velocity (c) Velocity- 

profile. gradient. profile curvature. 

Figure 6-13. — Velocity distribution in a boundary layer with pressure increase. 



Figure 6-14. — Pressure distribution on a turbine blade. 
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(positive dp/dx), the velocity profile in the boundary layer will have a 
point of inflection. Since the velocity profile at the point of separation 
(with du/dy = Q at the surface) must have a point of inflection, it follows 
that, with the assumptions used in deriving equation (6-23), separation 
can occur only when the potential flow is retarded (i.e., in regions of 
adverse pressure gradient) . 

Figure 6-14 indicates a typical pressure distribution on the surface of a 
turbine blade. The danger zone, as far as separation is concerned, is 
readily seen to be the rear portion of the suction surface, where the major 
part of the blade diffusion is taking place. 


Laminar Incompressible Boundary Layer on a Flat Plate 

Prandtrs boundary-layer theory was first reported in 1904 in Germany. 
It was later translated and published in 1928 as an NACA Technical 
Memorandum (ref. 3). The first mathematical solution of Prandtl's 
equations to be published was the flat-plate solution of Blasius in 1908. 
This German work was also later translated by NACA (ref. 4). 

On a flat plate with steady flow at zero incidence, the velocity from the 
potential solution is constant. Therefore, p{x) is constant and dp/dx = 0. 
The boundary-layer equations, therefore, reduce to 


du du d^u 

u [-V — = v — 

dx dy dy^ 


(6^59) 


where p is the kinematic viscosity ft/p, in mVsec or ftVsec, and 


du dv 
dx dy 

The following are the boundary conditions: 
u = i^ = 0 at v = 

u = Ue at y= 


(6-60) 


(6-61) 


With the use of a stream function Blasius transformed the partial 
differential equation (6-59) into the following ordinary differential 
equation : 


. dj d^f 

where /is a normalized stream function 


fiv) = 


Vi 


PXUe 


(6-62) 


(6-63) 
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which depends on the dimensionless y-coordinate, where 


V 



(6-64) 


This equation has the following boundary conditions: 


= 0 at 77 = 0 
dtj 


— = 1 at 77 = 00 

dy 


(6-65) 


Equation (6-62) cannot be solved exactly. Blasius obtained an approxi- 
mate solution in the form of a power series expansion about 77 = 0 and an 
asymptotic expansion for 77 = qo , the two solutions being joined at a suit- 
able point. More recently, Howarth (ref. 5) solved the Blasius equation 
(6-62) with a high degree of accuracy, and provid(*d tabular values for/, 
df/dy, and d‘^f/di/ as functions of 77. Since df/dy — ufue^ the solution gives 
the velocity profile of figure 6-15. This profile possesses a very small 
curvature at the wall and turns rather abruptly further from it in order 



Figure 6-15. — Blasius-Howarth velocity profile for flow on a flat plate. 
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to reach the asymptotic value. At the wall itself, the curve has a point of 
inflection, since for y = 0, d^u/dy^=0. 

From the order-of-magnitude analysis performed to obtain Prandtl’s 
boundary-layer equations, we had the relation 

Re 

(6-66) 

For a semi-infinite flat plate, the Reynolds number can be expressed as 

Rex = — 

V 

(6-67) 

In order to make equation (6-66) dimensionally correct, 

we can say 

X* Rcx 

(6-68) 

or 


^ Ue 

(6-69) 

The constant of proportionality can be obtained from Howarth's numerical 
solution and is equal to 5. So, for a semi-infinite flat plate at zero incidence 
in laminar flow, we obtain the useful relation for the boundary-layer 
thickness 

3/«ii = 5.0 V- 

^ Ue 

(6-70) 

With the use of Howarth^s solution to the Blasius equations, the follow- 
ing relations for other important boundary-layer parameters for laminar 
flow on a flat plate can also be obtained: 

5 = 1.72 

u. 

(6-71) 

e=oMiJ^ 

^ U. 

(6-72) 

^ = 0.332^=®-/.^ 

PM« ’ UeX V RCz 

(6-73) 

r. 1-328, , 

D “ h\/ fiplUe^ 

g 

(6-74) 
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C/= 1.328 


f7 1.328 


(6-75) 


where 

Tw shear stress on the surface, N/m-; Ibf/ft- 
D total drag on both sides of flat plate, N ; Ibf 
h width of flat plate, m; ft 

I length of flat plate, m; ft 

Cf dimensionless drag or skin friction coefficient for flat plate 
Rei Reynolds number based on plate length I 

It should be noted that all of these relations arc valid only for laminar 
flow; that is, they arc valid only where 10®, a value that is indicative 
of laminar flow over the entire length of the plate. For Rei'> 10®, transition 
to turbulent boundary layer will probably occur, and the expressions in 
equations (6~71) to (6-75) will be valid only from the leading edge of the 
plate to the transition point. If transition to turbulent boundary layer does 
occur, then the drag will be larger than that calculated by equation 
(^74). 

Integral Methods for Solving the 
Laminar-Boundary-Layer Equations 

The two principal means of solving the laminar-boundary-layer equa- 
tions are by integral methods and by fin it e-difference methods. Both 
means provide approximate solutions, since exact solutions are extremely 
cumbersome. 

Integral methods are based on von Karmen’s momentum integral 
formula. Von Karman^s original work was published in 1912 in Germany 
and was later translated by NACA (ref. 6). Von Kdrmdn realized that it 
was not necessary to satisfy the boundary-laj^er equations for every fluid 
particle. Instead, he satisfied the boundary-layer equations close to the 
wall and in the region where external flow is approached by satisfying the 
boundary conditions. In the remaining region of fluid in the boundary 
layer, only a mean over the differential equation is satisfied. Such a mean 
is obtained from the momentum equation (eq. (6-23) or (6-30)) by 
integration over the boundary-layer thickness. If equations (6-23) and 
(6-30) are integrated from ?/ — 0 to y = ^juih ^*^d if the definitions of dis- 
placement thickness (eq. (6-50)) and momentum thickness (eq. (6-54)) 
are introduced, the following equations result. For laminar, incom- 
pressible flow. 


= — ( 6 - 76 ) 

dx dx p 
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For laminar, compressible flow, 

(6-77) 

dx dx Pe 

where the subscript e denotes conditions at the outer edge of the boundary 
layer. 

Equation (6-76) or (6-77) leads to an ordinary differential equation 
for the boundary-layer thickness, provided that a suitable form is as- 
sumed for the velocity profile, u/ue. This allows us to calculate the dis- 
placement thickness, 5, the momentum thickness, 0, and the shearing 
stress at the wall, tu,. Pohlhausen was the first to use eciuation (6-76) to 
obtain a solution for incompressible boundary layers. His work was 
published in 1921 (refs. 7 and 1). The velocity profile assumed by 
PohIhaus(Ti was discussed earli(T in this chapter, under ‘‘V(^locity Pro- 
files.’" Although Pohlhausen’s solution is probably the simplest, it is known 
to give poor results in regions of rising pressure. As a rc'sult, various authors 
have tried to improve and extend his method by assuming different 
families of velocity distributions. 

A famous work among those which followed Pohlhausen’s was that of 
Thwaites (ref. 8). Thwaites collected and compared all known velocity 
distributions from exact and approximate solutions for laminar incom- 
pressible flow. Thwaites" method does not require the solution of ordinary 
differential equations. He relates the wall sh(‘ar, its derivative at the wall, 
and the form factor to one another without specifying a type of velocity 
profile. To do this, nondimensional forms of these quantities were defined 
and evaluated with the use of exact solutions for the laminar boundary 
layer. It developed that a nearly universal relation existed among these 
quantities for favorable pressure gradients. For adverse gradients, 
Thwaites selected a single representative relation, A unique correlation 
was chosen that reduced the solution of an incompressible problem to the 
evaluation of a single integral. Thwaites" method was extended to com- 
pressible fluids by Pott and Crabtree (ref. 9) . They recognized that when 
heat transfer is negligible, and the Prandtl number is equal to 1, a trans- 
formation proposed by Stewartson (ref. 10) could be used to relate 
compressible to incompressible boundary-layer solutions. 

One of the best integral methods to appear to date for the solution of 
laminar boundary layers is that of Cohen and Reshotko (refs. 11 and 12). 
Their method applies to compressible or incompressible flow over two- 
dimensional or axially symmetric surfaces. It handles arbitrary free- 
stream pressure distribution and performs well in areas of adverse pres- 
sure gradient. A surface temperature level may be specified, and heat 
transfer is calculated. Cohen and Reshotko"s method is based on Thwaites" 
correlation concept, Stewartson"s transformation (ref. 10) is first applied 
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to Prandtrs equations. The resulting nonlinear, first-order differential 
equations are expressed in terms of dimensionless parameters related to 
the wall shear, the surface heat transfer, and the transformed free-stream 
velocity. Then Thwaites' concept of a unique interdependence of these 
parameters is assumed. The evaluation of these quantities is then carried 
out by utilizing the exact solutions of reference 11. With the resulting 
relations, methods are derived for the calculation of all the important 
boundary-layer parameters. In 1960, Luxton and Young published a 
method (ref. 13) which is as general as Cohen and Reshotko’s, but which 
allows the Prandtl number to have values slightly different from L 


Finite-Difference Methods for Solving the 
Laminar-Boundary-Layer Equations 

Finite-difference methods for solving the boundary-layer equations 
have recently come into prominence because of the development of 
digital computers. Smith and Clutter have done a considerable amount of 
work in developing this technique (refs. 14 and 15). Another recent 
reference of interest is that of Krause (ref. 16). These methods give very 
good results with relatively short running times on the computer. 


Eddy-Viscosity and Mixing-Length Concepts in 
Turbulent Boundary-Layer Flow 

Before referencing any of the current methods for solving turbulent 
boundary-layer flow, the concepts of “eddy viscosity’’ and ‘^mixing 
length” should be discussed. These approximation concepts have been 
used in many of the methods developed to date to relate the Reynolds 
stresses produced by the mixing motion to the mean values of velocity 
components. By this means, the Reynolds stresses arc given a mathe- 
matical form which, upon substitution into the governing equations, 
leads to differential equations containing only mean values of density, 
velocity, and pressure. These transformed differential equations con- 
stitute the starting point for the calculation of the mean boundary-layer 
flow. 

Boussinesq first worked on this problem in 1877. In analogy with the 
coefficient of viscosity in Stokes’ law for laminar flow 


^ du 

g dy 


(6-78) 


where n is the laminar shear stress, in N/m^ or Ibf/ft^, he introduced a 
mixing coefficient, Ar, for the Reynolds stress in turbulent flow by putting 


Ar du 

g dy 


(6-79) 


84 



INTRODUCTION TO BOUNDARY-LAYER THEORY 


where t, is the turbulent shear stress, in N/m* or Ibf/ft^ In 1880, Reynolds 
introduced the concept of eddy, or virtual, viscosity, t, where 


Ar 

< = — (6-80) 

p 

Thus, the eddy viscosity is analogous to the kinematic viscosity v = fi/p. 
Turbulent stress can then be expressed as 


p du — p 

r< = -e— = — mV ( 6 - 81 ) 

9 9y g 

With the use of this concept, terms in equations (6-37) and (6-40) such as 


d 

dy 




du 

dy 



can be written as 


— (p-fA,) 

dy L 




A similar concept can be applied to the energy equation where an eddy, 
or a virtual, conductivity can bo defined. The difficulty with applying the 
eddy-viscosity method is that A, and hence t depend on velocity. It is, 
therefore, necessary to find empirical relations between these coefficients 
and the mean velocity. 

In 1925, Prandtl introduced a completely different approximation for 
the Reynolds stresses. His argument is called PrandtFs mixing-length 
hypothesis, since the mixing length is somewhat analogous to the mean 
free path in the kinetic theory of gases. The main difference is that kinetic 
theory concerns itself with the microscopic motion of particles, whereas 
Prandtrs concept deals with the macroscopic motion of large clusters of 
fluid particles. Deriving PrandtPs expression for shear stress requires a 
good deal of discussion of his physical model of turbulent flow, all of which 
is contained in reference 1. His final expression is 


Ti 




du 

dy 


du 

dy 



(6-82) 


where I is the mixing length, in m or ft. 

On comparing PrandtPs expression (eq. (6-82)) with that of Boussinesq 
(eq. (6-81)), it appears that little has been gained. The unknown eddy 
viscosity e of the first expression has merely been replaced by the unknown 
mixing length I of the second expression. However, PrandtPs equation 
(6-82) is generally more suitable for the calculation of turbulent motion 
than is equation (6-81). Turbulent drag is roughly proportional to the 
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square of velocity , and the same result is obtained from (6-82) if the 
mixing length is assumed to be independent of the magnitude of velocity. 
So, mixing length is a purely local function, although we cannot say it is a 
property of the fluid. It is far simpler to make assumptions about the 
mixing length I than about the eddy viscosity t, and this constitutes the 
superiority of Prandtl’s expression over that of Boussinesq. 

Integral Methods for Solving the 
Turbulent Boundary-Layer Equations 

Just as with the laminar-boundary-layer equations, there are both 
integral methods and finite-difference methods for solving the turbulent 
boundary-layer equations. Both of these provide approximate solutions, 
since exact solutions for turbulent flow arc now impossible. 

Gruschwitz was the first to propose a method for solving the equations 
for an incompressible turbulent boundary layer. His work was published 
in Germany in 1931. A rash of works followed, most of them making 
improvements to the calculational technique and empirical data used by 
Gruschwitz. Ludwieg and Tillmann, whose work was published in Ger- 
many in 1949 and was translated by NACA in 1950 (ref. 17) , proposed an 
empirical relation for the skin-friction term in the momentum integral 
equation. This relation is still ustd in many current methods. Stewartson’s 
transformations (ref. 10) arc likcwusc used in many methods for solving 
the turbulent-boundary-layer equations. 

IMaskell, in 1951 (ref. 18), proposed an improved method for incom- 
pressible turbulent boundary layers. He replaced the momentum equation 
by an empirically determined approximation which is directly integrable 
and thus determines the momentum thickness. A profile parameter is 
obtained from an empirical auxiliary differential equation. The Ludwieg- 
Tillmann skin-friction formula is used to calculate the skin-friction dis- 
tribution and to determine a separation point for flows with adverse 
pressure gradient. 

Truckenbrodt, whose work was published in Germany in 1952 and was 
translated by NACA in 1955 (ref. 19) , proposed solutions for both laminar 
and turbulent incompressible boundary-layer flows. 1 he method is simple 
and, like Haskell’s method, does not use the momentum integral equa- 
tion. It applies to both two-dimensional and rotationally symmetrical 
flows. Because of its simplicity and relatively accurate results, Trucken- 
brodt’s method is still used for incompressible turbulent boundary layers. 

Compressible turbulent boundary layers were first treated adequately 
with the use of integral methods by Reshotko and Tucker in 1957 (ref. 
20) . Prior to their work, the Karman momentum integral equation had 
been utilized with an assumed boundary-layer velocity profile, usually the 
power law', and one of several empirical skin-friction relations. When 
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pressure gradient was present, an auxiliary equation, usually the moment- 
of-momentum equation, was used. (This equation is obtained by mul- 
tiplying the integrand of the momentum integral equation by a distance 
normal to the surface and then integrating with respect to that distance.) 
The momentum integral equation and the auxiliary equation were then 
solved simultaneously. 

Reshotko and Tucker^s method, applicable to compressible flow with 
heat transfer and pressure gradient, also uses the momentum and moment- 
of-momentum integral equations. These arc expressed in incompressible 
form and are uncoupled with the use of Stewartson's transformation 
(ref. 10) and the results of Maskcll (ref. 18). The Ludwieg-Tillmann 
skin-friction relation is used in a form suitable for compressible flow with 
heat transfer through application of Eckert's reference-enthalpy concept 
(ref. 21). An approximation for the shear-stress distribution through the 
boundary layer and the power-law velocity profile are used to simplify the 
moment-of-momentum equation. Separation is located as the point where 
the skin friction, when extrapolated, becomes zero. This method, until 
several years ago, ^vas the best available for compressible turbulent 
boundary layers. It is still widely used in many computer programs today. 

One of the best integral methods available today for compressible 
turbulent boundary layers is that of Sasman and Cresci (ref. 22). It is 
simply an extension of the Reshotko-Tucker method. It uses somewhat 
the same analysis, but no attempt is made to uncouple the momentum 
and moment-of-momentum integral equations. These equations are solved 
simultaneously after introduction of boundary-layer shear-stress dis- 
tributions obtained from recent numerical results of equilibrium turbulent 
boundary-layer analysis. The Sasman-Cresci analysis is better than that of 
Reshotko-Tucker at predicting separation in regions of adverse pressure 
gradient. McNally (ref. 23) has developed a computer program based on 
the Cohen-Reshotko (refs. 11 and 12) and Sasman-Cresci (ref. 22) tech- 
niques. An additional source of information on compressible turbulent 
boundary-layer analysis is the work of Herring and Mellor (ref. 24). 

Finite-DifFerence Methods for Solving the 
Turbulent Boundary-^Layer Equations 

Finite-difference methods for solving the turbulent boundary-layer 
equations have recently begun to appear. Cebeci and Smith have done a 
large portion of this work to date (refs. 25, 26, and 27) . Bradshaw, Ferriss, 
and Atwell have also developed methods for the turbulent boundary layer 
(refs. 28 and 29) based on the use of the turbulent energy equation. 
Patankar and Spalding have developed still another method for handling 
the turbulent boundary-layer equations (refs. 30 and 31). A great deal of 
work is going on in this field at the present time, and no method is yet 
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clearly superior to any of the others. Two relatively recent publications 
(refs. 32 and 33) compare many of the most prominent methods, both 
integral and finite difference, for solving the turbulent boundary layer. 


CONCLUDING REMARKS 

The selection of a method of solution suitable to a particular boundary- 
layer problem requires some familiarity with the various methods avail- 
able. This can be achieved by studying some of the more recent references 
that have been mentioned herein. The present discussion of the methods 
of solution has been intended to show the historical development of 
solution techniques, the variety of methods available, and the complexity 
of the whole boundary-layer problem, especially where turbulent flows 
are involved. 
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SYMBOLS 

Ar turbulent flow mixing coefficient, (N) (sec)/m**; lbm/(ft) (sec) 
a constant in eq. (6-44) 

^ /width of flat plate, m; ft 
(constant in eq. (6-44) 

C f skin-friction coefficient for a flat plate 

c constant in eq. (6-44) 

Cp specific heat at constant pressure, J/(kg)(K); Btu/(lbm) (°R) 
D total drag on flat plate, N ; Ibf 
d constant in eq. (6-44) 

/ Blasius dimensionless stream function defined by eq. (6-63) 

f general body force vector, N/kg; Ibf/lbm 

fx component of body force f in x-direction, N/kg; Ibf/lbm 

fy component of body force f in ^-direction, N/kg; Ibf/lbm 

fz component of body force f in ^-direction, N/kg; Ibf/lbm 

g conversion constant, 1; 32.17 (Ibm) (ft) /(Ibf) (sec^) 

H form factor, defined by eq. (6-56) 
i unit vector in the x-direction 

J conversion constant, 1 ; 778 (ft) (Ibf) /Btu 
j unit vector in the ?/-direction 

k thermal conductivity, W/ (m) (K) ; Btu/(sec) (ft) (°R) 

k unit vector in the ^-direction 

L characteristic length (e.g., the blade chord) , m ; ft 
^ fPrandtl mixing length, m; ft 

/length of flat plate, m; ft 

Mg Mach number external to the boundary layer 
n exponent on the turbulent velocity profile, eq. (6-47) 

P dimensionless pressure, defined by eq. (6-13e) 
p static pressure, N/m^; Ibf/ft^ 

R gas constant, J/( kg) (K) ; (ft) (lbf)/(lbm) (°R) 

Re Reynolds number based on L and C/o, as defined by eq. (6-13f) 
Rei Reynolds number based on Z, as defined in eq. (6-75) 

Rex Reynolds number based on x, as defined by eq. (6-67) 

r radius of curvature of blade surface, m; ft 

S constant in eq. (6-25), K; °R 
T absolute static temperature, K; °R 

T t absolute total temperature, K; °R 

7^0 reference temperature used in eq. (6-25) , K; °R 
t time, sec 

U dimensionless velocity in x-direction, defined by eq. (6-13c) 

C/o free-stream velocity upstream of blade, m/sec; ft/sec 

u component of general velocity vector u in the .x-direction, m/sec; 

ft/sec 
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general velocity vector, m/sec; ft/sec 

free-stream velocity at the outer edge of the boundary layer, 
m/sec; ft/sec 

dimensionless velocity in i/-direction, defined by eq. (6-1 3d) 
component of general velocity vector u in the y-direction, m/sec; 
ft/sec 

component of general velocity vector u in the ^-direction, m/sec; 
ft/sec 

dimensionless x-coordinate, defined by eq. (6-1 3a) 
fx-coordinate, m ; ft 

^coordinate parallel to boundary surface, m; ft 
dimensionless y-coordinate, defined by eq. (6-13b) 
f ^-coordinate, m; ft 

\coordinate perpendicular to boundary surface, m ; ft 
z-coordinate, m ; ft 
displacement thickness, m ; ft 
boundary-layer thickness, m; ft 
( eddy viscosity defined by eq. (6-80) , m^/ sec ; ftV sec 
(a dimensionless quantity much less than 1 
Blasius transformed y-coordinate defined by eq. (6—64) 
momentum thickness, m; ft 

dimensionless shape parameter defined by eq. (6-45) 
dynamic viscosity, (N) (sec)/m’^; Ibm/ (ft) (sec) 
dynamic viscosity at reference temperature To, (N) (sec)/m®; 
lbm/(ft) (sec) 

kinematic viscosity, mVsec; ftVsec 
density, kg/m*; Ibm/ft* 

free-stream density external to the boundary layer, kg/m*; 
Ibm/ft* 

laminar shear stress, N/m*; Ibf/ft^ 
turbulent shear stress, N/m*; Ibf/ft^ 
shear stress at the wall, N/m*; Ibf/ft^ 
function defined by eq. (6-34) 

Blasius stream function, mVsec; ftVsec 
constant in eq. (6-26) 


Superscripts: 

time average 
fluctuating component 
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CHAPTER 7 


Boundary-layer losses 


By Herman W. Pnist, Jr. 


The primary cause of losses in a turbine is the boundary layer that 
builds up on the blade and end-wall surfaces. In particular, these losses 
are the friction loss resulting from the flow of the viscous fluid over the 
surfaces, the pressure-drag loss resulting from the flow of fluid past the 
blade trailing edge, and the loss downstream of the blades resulting from 
the mixing of the low-velocity boundary-layer fluid with the high-velocity 
free-stream fluid. Chapter 6 presented an introduction to boundary-layer 
theory, by means of which the surface boundary-layer buildup can be 
analytically described. This chapter covers analytical and experimental 
methods for determining the friction, trailing-edge, and mixing losses 
associated with the boundary layer. The theory presented herein refers 
primarily to two-dimensional blade-section boundary layers. Methods for 
obtaining three-dimensional blade plus end-wall losses from the two- 
dimensional results are also discussed. 

A fundamental objective in blade-row design is to minimize the energy 
loss resulting from the flow of fluid through the blade row. Therefore, the 
final expressions for loss developed in this chapter are in terms of kinetic- 
energy loss coefficients. These coefficients express the loss in fluid kinetic 
energy as a fractional part of the ideal kinetic energy of the actual flow 
through the blade row. Efficiency based on kinetic energy can be obtained 
by subtracting these coefficients from unity, and this is consistent with 
the blade-row efficiency definition used in chapter 2. 

Before proceeding with the discussion of boundary-layer parameters 
and loss coefficients, the blade-row station locations and associated 
pressure and velocity distributions will be introduced with the aid of figure 
7-1. These pressure and velocity distributions and the associated dis- 
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cussion refer to an attached boundary layer only. A separated boundary 
layer, with its associated reversal of flow at the surface, is thicker, yields a 
higher loss, and cannot be analyzed in the same manner, if at all. 

Figure 7-1 (a) indicates the four station locations that will be referred 
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Figure 7-1. — Station locations and associated pressure and velocity distributions. 
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to in this chapter. Station 0 represents the inlet to the blade row. At this 
station, a uniform total pressure po' is assumed, as indicated in figure 
7-1 (b). Station la is just within the trailing edge of the blade. The 
boundary layers developed on the blade surfaces result in velocity and 
pressure profiles as shown in figure 7-1 (b). Velocity varies from the free- 
stream value F/,,ia to zero at the blade surfaces. There is, of course, no 
flow through the region of the solid trailing edge. Total pressure varies 
from the frcc-stream value = Po' to the static pressure vu at the blade 
surfaces. This static pressure is assumed constant across station la, as is 
the flow angle «ia. At station la, only the surface friction loss has occurred. 

Station 1 is just beyond the blade trailing edge, where the boundary- 
layer fluid has filled the void, but where little mixing with the free stream 
has occurred. This is indicated in figure 1(b) by the station-1 profiles 
showing flow throughout the entire wake region. Here too, static pressure 
and flow angle are assumed constant across the station. Between stations 
la and 1, the traiiing-edge loss occurs. Station 2 is located at a distance 
sufficiently downstream of the blade row that complete mixing, with the 
associated mixing loss, has taken place. The velocity and total-pressure 
profiles are again uniform. 

In order to simplify analysis and discussion, a number of variables 
have been assumed constant across the various stations. Uniformity of 
inlet conditions is a universal convenience that usually can be approached 
in component tests but seldom exists in actual applications. Experiments 
have shown that static pressure and flow angle do vary somewhat across 
both free stream and boundary layer at stations la and 1. In some in- 
stances, which will be later identified, this variation can be accounted for. 
Although some downstream mixing of the flow does take place, a com- 
pletely uniform downstream state is merely a hypothetical convenience. 

BOUNDARY-LAYER PARAMETERS 

When a real fluid flows over a surface, a loss results due to both friction 
between the fluid and the surface and friction between the layers of fluid 
in the region adjacent to the surface. As shown by figure 7—2, the fluid 
velocity in the boundary-layer region varies from zero velocity on the 
surface to free-stream velocity F/, at the full boundary-layer height 
bfuw To describe the losses in flow, momentum, and energy resulting from 
the presence of the boundary layer, certain parameters are used. Some of 
these (displacement thickness, momentum thickness, form factor) were 
introduced in the last chapter and will be reviewed here; in addition, others 
specifically used for obtaining the desired kinetic-energy coefficients will 
be introduced and defined. 
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Free-stream velocity. 



Figure 7-2. — Typical boundary-layer velocity profile. 


The displacement thickness 6, which is indicative of the loss in mass 
flow, is defined by 


KpV)f.= / 

•' n 


(pV)„dY 


•'a 


(pF) dY 


(7-1) 


where 


5 displacement thickness, m ; ft 

dfuii boundary-layer thickness, m ; ft 

V fluid velocity, m/sec; ft/sec 

p fluid density, kg/m®; Ib/ft® 

Y distance in direction normal to boundary layer, m ; ft 

( ) /, free-stream (ideal) conditions 


Equation (7-1) states that the loss in mass flow of the fluid in the bound- 
ary layer is equal to the ideal flow which would pass through a length 
(or an area) equal to the displacement thickness. Solving for 5 yields 


5 = 




pV 

(pF)/. 


dY 


(7-2) 


The momentum thickness which is indicative of the momentum 
loss, is defined by 


e{pV^)/.= {pVV^.)dV- pV^dY (7-3) 

•'o •'o 


where 6 is the momentum thickness, in m or ft. Equation (7-3) states that 
the loss in momentum of the fluid in the boundary layer is equal to the 


96 



BOUNDARY-LAYER LOSSES 


ideal momentum of the ideal flow which would pass through a length (or 
an area) equal to the momentum thickness. Solving for 6 yields 


6 = 



pV 

(pV)f. 



pV^ 

(pV^)f. 


dY 


(7-4) 


The loss in kinetic energy can be similarly expressed in terms of an 
energy thickness defined by 

1 1 rhull 1 rhull 


where is the energy thickness, in m or ft. Equation (7-5) states that 
the loss in kinetic energy of the fluid in the boundary layer is equal to the 
ideal kinetic energy of the ideal flow which would pass through a length 
(or an area) equal to the energy thickness. Solving for yields 



pV 

ipV)„ 





dY 


(7-6) 


Ratios of the aforementioned thickness terms are also used as basic 
boundary-layer parameters. The form factor H is defined as 



(7-7) 


Substituting equations (7-2) and (7-4) into equation (7-7) and defining a 
dimensionless distance y as 


yields 



(7-8) 


(7-9) 


An energy factor E is defined as 



(7-10) 


Substituting equations (7-6), (7-4), and (7-8) into equation (7-10) 
yields 
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/■' pV , (' tV' , 

i 

E=~ (7-11) 

C pV , /•* pV^ . 


Velocity profiles for turbulent flow are often represented by a power 
profile of the type 


( 7 - 12 ) 


where the value of the exponent n is most often between 0.1 and 0.25. 

Note that this power profile is here expressed as while the same 
profile in chapter 6 (eq. (6-^7)) is expressed as The exponent ex- 
pressed as 1/n is consistent with general boundary-layer theory usage. 
The exponent expressed as n, however, is consistent with reference 1, 
wherein the equations that follow are derived. Therefore, the specific 
numerical value to be used for n will depend on the form being used for 
the exponent. 

Withthis velocity profile, equations (7-9) and (7-11) can be integrated 
in series form, and the form and energy factors for turbulent flow can be 
expressed in terms of the exponent ?i and the free-stream critical velocity 
ratio V/Ver- The resulting equations derived in reference 1 are 


1 I 3A,. 5A^. 

n+1 3n+l 5n+l 

H = (7-13) 

I ^ I ^ l_ 

(n+l)(2n+l) (3n + l)(4n+l) (5n+l)(6n4-l) 


and 


4^ 


L(n+l)(3n+l) (3n+l)(5n-|-l) (5n+l)(7n+l) 


+ 


E = 


4^ 


(n+l)(2n+l) (3n+l)(4a+l) (5n+l)(6n+l) 


+ 


where 



(7-14) 

(7-15) 
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and y is thf ratio of specific heat at constant pressure to specific heat at 
constant vjlume, and Vcr is the fluid velocity, in m/sec or ft/sec, at the 
critical (Mach 1) flow condition. For incompressible flow, \\^here V/Vcr 
approaches zero, equations (7-13) and (7-14) reduce to 

^inc = 2n+l (7-16) 

and 


Einc=^ 


2(2n+l) 

3n+l 


(7M7) 


Values of the form and energy factors for turbulent compressible flow 
are shown in figure 7-3 for V /Vcr varying from 0 to 1.4 and n varying 
from 0 to 1.5. It can be seen that the form factor varies much more than 
does the energy factor. For any constant exponent n, the energy factor 
is almost independent of F/ F^r. 

The boundary-layer parameters just presented are general and can 
refer to a boundary layer on any type of body. They are directly useful in 
certain aerodynamic work. For instance, the drag of a body can be ob- 
tained directly from the momentum thickness. In turbine work, however, 
where the flow is confined to the physical boundaries of the blade row, it is 
simpler and more meaningful to express the losses as a fractional part of 
the ideal quantities that could pass through the blade row. The thickness 
parameters so expressed are herein termed “dimensionless thickness 
parameters’^ and are defined on the basis of zero trailing-edge thickness. 


Free -St ream 



Figure 7-3. — Effect of compressibility on variation of energy factor with form factor. 

(Data from ref. 1.) 
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These dimensionless thickness parameters must represent the sum of the 
suction- and pressure-surface thicknesses. 

With the assumption that flow conditions in all channels are the same, 
the dimensionless thickness parameters are obtained by dividing the 
losses in flow, momentum, and energy for a single blade-row channel by 
the corresponding ideal quantities that could pass through one blade- 
row channel. The total losses for one channel, as indicated in figure 7-4, 
are composed of the suction-surface loss plus the pressure-surface loss, or 



(7-18) 

= + 

(7-19) 

^tot = 4^M+ypp 

(7-20) 

where the subscripts toty s, and p denote total value, suction-surface value, 
and pressure-surface value, respectively. Thus, in terms of the previously 
defined boundary-layer thicknesses, the dimensionless boundary-layer 
thicknesses are expressed as 

^iot 

scosa(pF)/, scosa 

(7-21) 

Otoi(pV^) fg Slot 

s cos a (pF *) /g s cos a 

(7-22) 

✓jv s cos a 

scosa f-j (pF*)/. 

(7-23) 
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where 

6* dimensionless displacement thickness 
6* dimensionless momentum thickness 

4^* dimensionless energy thickness 

8 blade spacing, m; ft 

a fluid flow angle from axial direction, deg 

Equations (7-21), (7-22), and (7-23) express the losses in flow, momen- 
tum, and energy, respectively, as fractions of their respective ideal quanti- 
ties for the blade row if the trailing-edge thickness is assumed to be zero. 
These equations can be subscripted to apply at either station la, within 
the trailing edge, or station 1 , beyond the trailing edge. 


BLADE-ROW LOSS COEFFICIENTS 

As mentioned previously, the losses are to be expressed in terms of 
kinetic-energy loss coefficients. In this section, methods for evaluating the 
friction, trailing-edge, and mixing losses and expressing them in terms of 
the kinetic-energy loss coefficients will be presented. 

Surface -Friction Losses 

The kinetic-energy loss coefficient eu, defined as the loss in kinetic 
energy as a fraction of the ideal kinetic energy of the blade-row actual 
flow, can be expressed in terms of the boundary-layer dimensionless 
thicknesses as 


cos au{pV^)f,.la 

(s cos aia-5u« cos au—t) (pT») /.,i„ 


(7-24) 


where t is the blade-row trailing-edge thickness, in m or ft. (Refer to fig. 
7-4 for the nomenclature in the region of the trailing edge of the blade.) 
Since this coefficient is referenced to station la, just within the blade-row 
trailing edge, it represents only the surface-friction loss. If a trailing-edge 
dimensionless thickness is expressed as 



S cos aia 


equation (7-24) reduces to 


(7-25) 


Sla = 


it 


(7-26) 


In order to evaluate the loss coefficient ci„ from equation (7-26), it is 
necessary to know the values of the dimensionless energy thickness 
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and the dimensionless displacement thickness These can be evaluated 
either experimentally or analytically, as will be discussed herein. 
Experimental determination . — In determining experimental loss values, 
it is impractical to measure the density and velocity directly. Instead, 
pressure data are taken, and the density and velocity are related to 
pressure functions. The pressure data required for computing the friction 
loss consist of (see fig. 7-1) the upstream total pressure po^ the blade-exit 
static pressure puy and the total-pressure loss survey data V^'^Via 
blade space. Since the dimensionless boundary-layer thicknesses express 
the losses of the blade row as a fractional part of the ideal quantities w^hich 
could pass through the blade row, the dimensionless displacement thick- 
ness can be expressed in terms of the flow across one blade pitch as 

5 cos au{pV)f,,u-i*s cos aia(pF)/,.u-cos aia ( (pF)ia du 

* •'o 

^ia= 7-^7^ (7-27) 

SCOSaia(pF)/.,i„ 

where u is the distance in the tangential direction, in m or ft. Equation 
(7-27) simplifies to 


In a similar manner, the dimensionless momentum and energy thick- 
nesses can be expressed as 


and 








(Px.V}.)u 

■/.'['-(aiteX'C) 


(7-30) 


Assuming that the total temperature T' and the static pressure pu in 
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the boundary layer are the same as in the free stream, the density ratio 
(p/p/*)ia can be related to the pressure ratio p[Jpo as follows: From the 
isentropic relation, 



(7-31) 

(7-32) 


Since P/,,u = Po^ and p[J pf,,ia^p[a/Po (from the ideal gas law, with 
= division of equation (7-31) by equation (7-32) yields 



(7-33) 


The velocity ratio (F/F/,)u can be related to the pressure ratios 
(p/p')ia and pia/pQ as follows: From the total-temperature definition and 
isentropic relation, equations (1-51) and (1-52) of chapter 1, we can 
write 


where 


F2 ^ T ^ {p 


(7-34) 


g conversion constant, 1; 32,17 (Ibm) (ft)/(ibf) (sec^) 

J conversion constant, 1; 778 (ft) (lb)/Btu 

Cp specific heat at constant pressure, J/(kg) (K) ; Btu/(lb) (®R) 

Subscripting equation (7-34) once for station la and again for free-stream 
values at la, dividing the first of these equations by the second, and 
recalling that P/,,ia = Po' and yields 



(7-35) 


With the density and velocity ratios expressed in terms of the measured 
pressures by equations (7-33) and (7-35), it is now possible to integrate 
equations (7-28), (7-29), and (7-30) and evaluate the dimensionless 
boundary-layer thicknesses. Then, the kinetic-energy loss coefficient eia 
can be computed from equation (7-26) . 

The kinetic-energy loss coefficient thus determined is a two-dimensional 


103 



TURBINE DESIGN AND APPLICATION 


coefficient; that is, it is based on data either from a two-dimensional 
cascade or from a constant radius of an annular cascade. (The annular 
cascade can be, and often is, the full stator or rotor from a turbine.) In 
order to obtain a three-dimensional loss coefficient for a blade row, data 
are taken at a number of radii sufficient to adequately cover the annulus, 
and the two-dimensional dimensionless boundary-layer thicknesses are 
calculated as shown previously for each radius. Three-dimensional bound- 
ary-layer thicknesses are then obtained by radial integration from hub 
to tip : 


^lo.SO ~ ■ 


/ »■< 

C(P^)/-.ia COS auT dr 
rh 

f ipV)f.. la COS aia r dr 


(7-36) 


I /..la COS au r dr 


— ■ 




COS aia r dr 


(7-37) 


^i«(p^’) /..1« cos «u r dr 

(7-38) 

(p^*) /..I. cos aur dr 


/ 


In terms of the measured pressures, these integrals are expressed as 

,(7-l)/-ylV2 




fU r \(7-1)/7-|1/2 

<p-)‘'’[i-(-:) J 


(7-39) 
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4'ta.lD — ■ 


,r, r /„ \(t-1)/t-13/* 

,r, r /„ \ (t-1)/t-]3/* 

/ (pJ»'[l-(-:) J »os».r* 


(7-40) 


(7-41) 


The three-dimensional kinetic-energy loss coefficient is then obtained 
in a manner similar to equation (7-26) : 


^la.SD 


*ta. 




(7^2) 


where is the trailing-edge dimensionless thickness at the mean radius 
and is used to represent the average value for the blade row. 

Analytical determination , — The kinetic-energy loss coefficient C\a can 
also be evaluated with the use of analytically determined boundary-layer 
thickness parameters. While not as reliable as experimental values, ana- 
lytical values are much less costly and time consuming to obtain. Ana- 
lytical methods for calculating the basic boundary-layer parameters are 
discussed and referenced in chapter 6. The boundary-layer solutions are 
not simple, and the better methods require computer solution. Boundary- 
layer computer programs currently in use at the NASA Lewis Research 
Center include one (ref. 2) based on an integral method solution and 
another based on the finite difference method of reference 3. 

An equation used in the study of reference 4 to compute turbulent 
boundary-layer momentum thickness was 

0.231 



( pV\ 


X 


r/ oF \ / F / y.JM 


(7-43) 
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where 

A/, parameter defined by equation (7-15) 

X distance along blade surface from forward stagnation point, m; ft 

I blade-surface distance from forward to rear stagnation point, 

m; ft 

M viscosity, (N) (sec)/m2; lb/(ft) (sec) 

The development of this equation is presented in reference 4. It is as- 
sumed that the boundary layer has a power-law velocity profile. In 
reference 4, the exponent n was obtained from the referenced equation 



Equation (7-43) must be evaluated for both the suction and pressure 
surfaces of the blade. The free-stream velocities and densities required for 
equations (7-^3) and (7-44) are those free-stream values adjacent to the 
blade-surface boundary layers. These can be obtained by any of the 
channel flow analysis techniques discussed in chapter 5. 

Values of the form factor H as required in equation (7^3) and of the 
form factor H i® and energy factor E\a at station la for each surface can be 
obtained from equations (7-13) and (7-14). With and Exa known 

for both the suction and pressure surfaces, the various boundary-layer 
thickness parameters and the kinetic-energy loss coefficient eia can be 
evaluated from the equations presented earlier in this chapter. For the 
turbine stator blade studied in reference 4, the analytical values, as 
calculated from equation (7-43), of the boundary-layer momentum 
thickness for the blade and for the two surfaces individually were reason- 
ably close to the experimental values. In general, however, results ob- 
tained from equation (7-^3) will not be as accurate as those obtained 
from the computer programs of references 2 and 3. 

Three-dimensional boundary-layer parameters could be calculated 
directly from equations (7-36) to (7-38). The two-dimensional thickness 
parameters would have to be analytically determined at a number of radii 
sufficient to establish the variation over the blade length and would also 
have to be determined, somehow, over the end-wall surfaces. Such a 
procedure would require considerable effort, so the simplified method of 
reference 5 for predicting three-dimensional losses from two-dimensional 
mean-section losses is commonly used. Results obtained by this method 
have shown good agreement with experimental results. 

In the method of reference 5, the following assumptions are made: 
(1) The average momentum loss for the blade surface can be represented 
by the dimensionless momentum thickness at the blade mean section; 
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Approximate area 



Figuke 7-5— Schematic diagram of equivalent two-dimensional blade used to calcu- 
late the effect of end-wall area on blade loss. 


(2) the momentum loss per unit area on the inner and outer end walls is 
the same as the average momentum loss per unit area on the blade surface; 
and (3) the blade configuration can be satisfactorily approximated by an 
equivalent two-dimensional blade, as shown in figure 7—5, having a con- 
stant cross section, spacing, and stagger angle equal to those at the mean 
section of the given blade. The surface area of one equivalent blade 
(see fig. 7-5) is 

Ab = 2ck (7-45) 

where 

Ah total surface area (sum of suction-surface and pressure-surface 
areas) of one blade, m^; ft^ 
c blade chord, m; ft 

h blade height, m; ft 

The inner and outer end-wall area for one passage is 

= 2sc cos ( 7-46 ) 

where 

Aw total surface area of passage end walls (sum of inner and outer 
end-wall areas) , m^; ft^ 
blade stagger angle, deg 

Now^, taking the average momentum loss 6*^^^ over the blade radial length 
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and modifying it to include the end-wall losses yields 

* ^ /Ab+Ay\ ^ / s cos aA 


The three-dimensional energy and displacement thickness parameters 


are then calculated as 




’/'la, 3D ~ 

(7-48) 

and 




^la,3D ~ -m^la.3D 

(7-49) 


Mean-section values are used for the energy and form factors. Although 
the energy and form factors were originally defined in terms of individual 
boundary-layer thicknesses, it is indicated in reference 4 that they can be 
satisfactorily used as is done in equations (7-48) and (7-49). The three- 
dimensional kinetic-energy loss coefficient is then obtained from equa- 
tion (7-42). 


Trailing-Edge Loss 

The kinetic-energy loss coefficient €te that represents the loss associated 
with flow past the blade trailing edge can be determined either experi- 
mentally or analytically. 

Experimental determination . — Experimental values of blade trailing- 
edge loss coefficient are obtained from differences between experimental 
two-dimensional loss coefficients 6i, which include both surface-friction 
loss and trailing-edge loss, and loss coefficients which include only the 
blade surface-friction loss. Thus, 


Cfe e\ e\(i (7 50) 

Loss coefficients e\ay which include only surface-friction loss, are ob- 
tained as described previously. Loss coefficients Ci, which include both 
surface-friction loss and trailing-edge loss, are determined in exactly the 
same manner except that the total-pressure loss and static pressure are 
measured at different locations. The surface-friction loss coefficients were 
based on data obtained just within the blade trailing edge at station la, 
where the trailing-edge loss has not yet occurred. To determine the loss 
coefficients w^hich include both surface-friction loss and trailing-edge loss, 
the measurements must be made at a location just downstream of the 
blade row, corresponding to station 1 in figure 7-1, w here the trailing-edge 
loss, but little mixing, has occurred. 

Analytical determination . — In reference 6, experimental drag coefficients 
are presented for a large number of surface discontinuities. Included in 
the reference are experimental data for sheet-metal joints of different 
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geometry, bolt and rivet heads of different geometry, and airfoil trailing 
edges. It is indicated that the pressure-drag loss due to the discontinuity 
behaves similarly regardless of the flow direction over the discontinuity. 
Therefore, the loss due to flow past a trailing edge will be treated ana- 
lytically as if the loss were due to a body placed in the path of a boundary 
layer. 

As indicated in reference 6, the drag of a small body of height hj equal 
to or less than the full boundary-layer height d/uih placed in a turbulent 
boundary layer, as shown in figure 7-6, corresponds approximately to the 
effective dynamic pressure of the part of the boundary layer equal to the 
height of the body. Thus, 

D^q^ffhCn (7-51) 


where 


D drag on body, N/m; Ib/ft 

h height of body, m ; ft 

Cd drag coefficient 

and the effective dynamic pressure Qe// is expressed as 


Qeff 


hJ. 


pV^ 

2g 


dY 


Drag is related to momentum thickness as 


(7-52) 


(7-53) 


Therefore, a dimensionless momentum thickness 6*^ representing the 
trailing-edge loss is obtained by combining equations (7-51) and (7-53) 
with a jiroperly subscripted form of equation (7-22) : 
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^?e = - 


Qe/fhCD 


S COS ai 


(pn> 


(7-54) 


The flow angle ai is related to the angle au as discussed in chapter 4 (eqs 
(4-26) and (4-27)). 

Before equation (7-54) can be evaluated, the effective dynamic pres- 
sure must be determined. The ratio of the effective dynamic pressure to 
the free-stream dynamic pressure is equal to 





(7-55) 


For turbulent flow, the variation of velocity in the boundary layer can be 
expressed with the use of the simple power profile presented previously. 
Combining equations (7-8) and (7-12) yields 


V 

Vfs 



(7-56) 


Assuming that the total temperature and static pressure in the boundary 
layer and free stream are the same and using the ideal gas law and equa- 
tion (1-64) of chapter 1 gives 


p 


Til 

To' 


P/» 


T 

TV 



(7-57) 


Substituting equations (7-56) and (7-57) in equation (7-55) and using 
the parameter A/, defined by equation (7-15) yields 


Qeff 


1 

h 


a-Af.) 



Performing a binomial expansion and integrating then gives 


Qeff 

Qfs 


(1-^/.) 



A/s 

(4n+ 1) 


(7-58) 


+ 



4^ 




(6n+l) 


(7-59) 


no 
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Substituting equation (7-59) in equation (7-54) and using trailing-edge 
thickness t in place of body height h finally yields 


e 


* 

te 


iCp 

2s cos «i 


( 1 -^/.) 


( t Y 1 / 1 Y A„ 

\^fuii/ 2n+l \5fuii/ (4n+l) 

+ -f- 

WiiJ (6n+l) 



(7-60) 


The boundary-layer thickness 5fun to be used in equation (7-60) should 
be the sum of the suction- and pressure-surface values. 

Equation (7-60) is for compressible flow. In many cases, at least when 
n is not well known, the following simplified equation, which assumes 
incompressible flow and n = 3^ (commonly used for turbulent flow), 
is adequate : 


=0.375 


4 . 


tCi) 
dfuii S cos ai 


(7-61) 


The information in reference 6 indicates that the drag coefficient Cd can 
be set equal to 0.16 for a rounded trailing edge and 0.22 for a square 
trailing edge. The corresponding values reported in reference 7 and con- 
verted to the same basis as equation (7-61 ) arc 0. 14 for a rounded trailing 
edge and 0.22 for a square trailing edge. Frequently, 8 tot instead of 8/uu 
will be available. In such a case, for incompressible flow. 


/I \ 

full = ^tol ( 1 ) 

\n ) 


and for compressible flow, 

5/u n = 

1-(1-A/,) 


^tot 


\n-f- 


Af. Al, 

H — I ~ + 

+ 1 3n+l 5n+l 


•) 


(7-62) 


(7-63) 


Equations (7-60) and (7-61) give the fractional loss in blade-row 
momentum due to the blade trailing edge. To find the corresponding 
kinetic-energy loss coefficient, it is necessary to find the fractional losses 
in flow and kinetic energy. As a simple approximation, the form and energy 
factors, evaluated from equations (7-13) and (7-14) for compressible 
flow and from equations (7-16) and (7-17) for incompressible flow, are 
used to obtain 

8l = H0l (7-64) 

and 

K = (7-65) 
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At station 1, which is just downstream of the blade trailing edge, fluid has 
flowed into the area behind the trailing edge and there is no longer a void 
due to trailing~edge blockage. Therefore, a kinetic-energy loss coefficient 
is obtained as 

«<.=7%r (7-66) 

This loss coefficient expresses the loss in kinetic energy as a fraction of the 
ideal kinetic energy of the flow that would exist if the trailing-edge loss 
were the only loss. The trailing-edge kinetic-energy loss coefficient for 
incompressible flow is plotted against trailing-edge thickness in figure 7-7 
for several values of the ratio of trailing-edge thickness to boundary-layer 
thickness. This figure is based on the momentum loss as expressed by 
equation (7-61). The flow loss associated with blade-surface friction is 
not included in equation (7-66). Therefore, this trailing-edge kinetic- 
energy loss coefficient is approximately, but not rigorously, additive with 
the surface-friction loss coefficient. Expression of the combined friction 
and trailing-edge loss in terms of a kinetic-energy loss coefficient is dis- 
cussed in the next section. 



Figure 7-7. — Effect of trailing-edge blockage on kinetic-energy loss coefficient. Form 
factor = 1.3; energy factor E — 1.8; drag coefficient Ci> = 0.16. 
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Combined Friction and Trailing-Edge Loss 

As stated in the discussion of the trailing-edge loss, the boundary-layer 
thickness parameters and a kinetic-energy loss coefficient expressing the 
combined friction and trailing-edge loss can be obtained by making the 
experimental measurements at a location corresponding to station 1, 
which is just downstream of the trailing edge. In this way, we obtain 
experimental values of ^i*, and from appropriately subscripted 
versions of equations (7-28), (7-29), and (7-30). The value of ei is then 
obtained as 

= ^ (7-67) 


Analytically, the boundary-layer thickness parameters at station 1 are 
obtained by adding the surface-friction loss to the trailing-edge loss. 
Before the friction and trailing-edge boundary-layer thickness parameters 
can be added, they must be expressed on the basis of the same ideal flow. 
The friction-loss dimensionless thicknesses at station la (6f^, are 

expressed in terms of an ideal flow without trailing-edge blockage. How- 
ever, there is a trailing-edge blockage at station la, where the ideal flow 
with blockage must be comparable to the ideal flow at station 1 , where 
there is no blockage. Therefore, the friction-loss boundary-layer thickness 
parameters are adjusted to account for the true (with blockage) ideal 
flow as follows : 


and 


♦ / scosau \ 
cos 

_ / scosaig \ 

~ As cos ara-t) 


\5 cos aia^t/ 


(7-68) 

(7-69) 


(7-70) 


where the subscript/ refers to the loss due to surface friction. Adding the 
friction and trailing-edge loss parameters then yields the combined loss 
parameters at station 1 ; 


and 


«l*=c+c 

(7-71) 


(7-72) 


(7-73) 


And the value of ei is then obtained from equation (7—67). 
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After-Mix Loss 

The after-mix loss is the total loss that includes the surface-friction loss, 
the trailing-edge loss, and the mixing loss. The after-mix loss coefficient 
62 is determined as described in this section, and the mixing loss, if desired, 
is obtained by subtracting the previously determined ei from 62. 

To determine the after-mix loss experimentally would require that the 
pressure measurements be made downstream of the blading where com- 
plete mixing has occurred. This is impractical for several reasons: (1) 
The length for complete mixing, while quite long, is unknown; (2) the 
after-mix loss would have to be corrected for side-wall friction for the 
mixing length, thus leading to possible error; and (3) after the flow had 
mixed, values of after-mix po' — p 2 ^ would be constant and small enough 
that the possibility of measurement error would be relatively large. For 
these reasons, values of after-mix loss are obtained analytically with the 
use of either experimentally or analytically determined before-mix 
(station 1) loss parameters. 

The basic equations for determining the after-mix conditions are 
those for conservation of mass, momentum in the tangential direction, 
and momentum in the axial direction during mixing. Equating the mass 
flow rate before mixing (station 1) and after mixing (station 2) yields 


j (pF) 1 cos aid = cos a2(pF)2 (7-74) 

From conservation of momentum in the tangential direction we get 

j (pF^) 1 sin ai cos ai d ^-^ = sin aa cos a2(pF^)2 (7-75) 


and from conservation of momentum in the axial direction we get 

g j j (pF^) cos^ ai d ^“^ = gp2+cos^ a2(pF^)2 (7-76) 


Although these equations are subscripted for two-dimensional flow, they 
can also be applied to three-dimensional flow by integrating radially. 

If experimental survey data were available at station 1, the integrals in 
the above equations could be directly evaluated even with variations in 
static pressure and flow angle. These conservation equations could be 
written for any before-mix location at which data were available, and then 
used to evaluate the after-mix loss coefficient. In the case where the 
before-mix station is not station 1, it would not be possible to determine 
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the mixing loss completely by experimental means. In most cases, how- 
ever, it is only the final after-mix loss that is desired, and survey measure- 
ments are usually made a little farther downstream of the trailing edge, 
where angle and pressure variations have somewhat damped out. 

If static pressure and flow angle are constant across station 1, it is 
possible to express equations (7-74) to (7-76) in terms of the previously 
used boundary-layer parameters, as was done in reference 1. The analysis 
herein differs from that of reference 1 only in that the before-mix station 
used in reference 1 corresponds to station la herein. Equation (7-28) 
subscripted for station 1, where there is no trailing-edge void, can be 
written as 

i (7-77) 

Subscripting equation (7-29) for station 1 and combining it with equation 
(7-77) yields 

^ (pF2)id0^=(l-a,*-fli*)(pF*)/..i (7-78) 

Substituting equations (7-77) and (7-78) into equations (7-74) to 
(7-76) yields the following equations for conservation of mass and 
momentum in terms of the boundary-layer parameters previously deter- 
mined : 

cos ai ( 1 - «i*) (p7)/.,i = cos aj(pV)* (7-79) 

sin ai cos a,(l-5,*-fl,*) (pF*)/.,i = sin «, cos a»(pF*), (7-80) 

gpi+cos*ai(l-«i*-9,*) (pF*)/.,, = cos*a,(p7*),-f-grp, (7-81) 

These equations, along with the ideal gas law and the conservation-of- 
energy equation {T\ =Tt) , can be solved simultaneously as shown in 
reference 1 to obtain ej, the after-mix kinetic-energy loss coeflflcient, 
for both compressible and incompressible flow. 

For incompressible flow, the solution for St is 


sin* ai 


\ 1-Si* ) 


-f-COS*ai(l-«i*)* 


1-1-2 cos* a,[(l- 5 ,*)J- 


(7-82) 


For compressible flow, no explicit solution was found, and the following 
steps are required to obtain St : 
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( 1 ) The parameters C and D are computed from 


(I-A/.. 1 ) ^+cos2ai(l-5i* 
2y 




c=- 


cos ai(l 


-'•’(a. 




(2) The quantity ( VJV „)2 is obtained from 

(3) The density ratio {p/p')i is obtained from 

©■{-(SOMHr 

(4) The total pressure ratio P 2 VP 0 ' is obtained from 

cosai(l-5i*) 

p,' 

(pV.\ 

\p'vj, 

(5) The pressure ratio (p/p ')2 is obtained from 

(?).=e): 

(6) Finally, it is obtained from 

©y 


'\ (t-D/t 


-1 


^ = - 


(^) 


'\ (y-i)/7 


-1 


(7-83) 


(7-84) 


(7-85) 


(7-86) 


(7-87) 


(7-88) 


(7-89) 


Values of St include all the blade-row loss; that is, the frictional loss of 
the blade row, the trailing-edge loss, and the mixing loss. Values of Si 
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include all the blade-row losses except mixing loss. Therefore, 

= (7-90) 

where €miz is the fractional loss in available energy due to mixing. 

BLADE-ROW LOSS CHARACTERISTICS 

In this section, experimentally and analytically determined losses of 
the various types considered will be presented and compared, and the 
effect of blade-row geometry on losses will be discussed. 

Distribution and Comparison of Losses 

Figure 7-8, taken from reference 8, compares experimentally and 
analytically determined values of kinetic-energy loss coefficients at three 
different angle settings for a given stator and at three stations repre- 
senting different losses. The loss coefficient €io.m, obtained just within the 
blade trailing edge, represents the surface-friction loss at the mean 
(arithmetic mean radius) section; the coefficient Ci.m, obtained just 
beyond the trailing edge, represents the friction loss plus trailing-edge loss 
at the mean section; and the coefficient represents the total loss for 
the annulus including blade and end-wall fi;iction, trailing-edge drag, and 
mixing. In general, agreement between the experimental and analytical 
loss coefficients is reasonably good. 



Percent stator area setting 


Figure 7-8. — Comparison of experimental and analytical loss coefficients for different 
stator area settings. (Data from ref. 8.) 
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Mean-section ideal after-mix critical velocity, 

<V/Vcr'i.m.2 

Figure 7-9. — Variation of loss coefficients with velocity. (Data from ref. 9.) 


Figure 7-8 gives some idea of the distribution of losses in a stator blade 
row, but does not separate the mixing and end-wall losses. Figure 7-9, 
taken from reference 9, shows the variation in loss coefficient with ve- 
locity. Loss coefficient is seen to decrease slightly with increasing velocity. 
This figure also shows separately the mixing and end-wall losses, as well 
as the other blade-row losses. 

In this particular case, the friction loss was about 2 percent of the stator 
ideal energy and about one-half of the total stator loss. The t railing-edge 
loss w^as about one-quarter of the total loss. In general, the trailing-edge 
loss will vary with trailing-edge blockage as was shown in figure 7-7. 
The end- wall loss, which was about 15 percent of the total loss for this 
case, will vary with the design, depending primarily on radius ratio and 
spacing. The mixing loss made up the remaining 10 percent of the total 
loss. The loss breakdown will, of course, vary with the stator design, but 
the comparison does indicate that each of the losses may be of con- 
sequence. 

Effect of Blade-Row Geometry on Losses 

A study of the effect of turbine geometry on turbulent-flow boundary- 
layer loss is presented in reference 10. In that study, the assumption was 
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made that the momentum loss per unit blade surface varies as the inverse 
of the chord Reynolds number to the m power: 


— ^ReT (7-91) 

c 

where is a Reynolds number based on blade chord c. Expanding 
equation (7-91) by multiplying and dividing by like terms, substituting 
into equation (7-22), and then using equation (7-47) to express the 
three-dimensional effect yielded an equation of the form 



where Ren is a Reynolds number based on blade height h. As indicated, 
the three-dimensional momentum thickness parameter can be expressed 
as a function of the geometric variables — height-to-spacing ratio h/s, 
blade solidity c/s, and height Reynolds number Rch. The reference value 
of 0tot/c, as explained in reference 9, is based on the minimum loss for a 
given solidity and, therefore, becomes a function of solidity. The ex- 
ponent m is set equal to ^ in the analysis. 

In reference 10, the derivative of the dimensionless momentum thick- 
ness 0 * 2 ) with respect to each of the geometric variables was obtained in 
order to find the minimum-loss value of each variable in terms of the other 
variables (there is no minimum for height Reynolds number). With the 
optimum values known, the relative variations in momentum loss around 
the minimum values were then determined. The results of the analysis 
from reference 10 are shown in figures 7-10, 7-11, and 7-12. Also shown in 
each figure is the nature of the geometry variation associated with the 
change in each variable. 

Figure 7-10 shows that a wide variation (50 percent or more) in h/s 
value around the optimum causes little increase in momentum loss. This 
results from the two counteracting effects of changes in chord Reynolds 
number and end-wall area. Figure 7-11 shows that the solidity of a blade 
may be varied considerably from optimum with some, but not excessive, 
loss. Comparison of the results in figures 7-10 and 7-11 shows that the 
loss is more sensitive to solidity than to the height-to-spacing ratio. The 
curve shape of figure 7-11 reflects also the counteracting influences of 
chord Reynolds number and end-wall area. 
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HeighHo-spacing ratio relative to optimum ratio, 

Figure 7-10. — Variation of momentum-thickness ratio with variation in height-to- 
spacing ratio. (Data from ref. 10.) 



Figure 7-11. — Variation of momentum-thickness ratio with solidity ratio. (Data 

from ref. 10.) 
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Figure 7-12. — Variation of momentum-thickness ratio with height Reynolds number 

ratio. (Data from ref. 10.) 


Figure 7-12 shows the variation of momentum- thickness ratio with 
height Reynolds number ratio. While the figure indicates a change in 
Reynolds number due to change in geometry^ the change in Reynolds 
number could also result from change in inlet flow conditions. The curve 
shape, then, results from the loss being inversely proportional to the 
Reynolds number to the = power. These results show that an in- 
crease in height Reynolds number results in improved performance. The 
height Reynolds number is sometimes used in correlating the performance 
of different turbomachines. 
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SYMBOLS 

Ah surface area of one blade, m*; ft* 

Aft parameter defined by equation (7-15) 

Aw surface area of end walls for one passage, m*; ft* 

C parameter defined by equation (7^3) 

Cd drag coefficient 

c blade chord, m; ft 

Cp specific heat at constant pressure, J/ (kg) (K) ; Btu/(lb) (°R) 

^ fdrag, N/m;lb/ft 

\ parameter defined by equation (7-84) 

E energy factor 

g kinetic-energy loss coefficient 

g conversion constant, 1 ; 32.17 (Ibm) (ft)/(lbf) (sec*) 

H form factor 

^ f blade height, m; ft 

\height of body placed in boundary layer, m ; ft 
J conversion constant, 1 ; 778 (ft) (lb)/Btu 

I blade surface distance from forward to rear stagnation point, 

m; ft 

m exponent in equation (7-91) 

n turbulent boundary-layer velocity profile exponent 

p absolute pressure, N/m*; Ib/ft* 

q dynamic pressure, N/m*; Ib/ft* 

Rce chord Reynolds number 

Rch height Reynolds number 

r radius, m ; ft 

8 blade spacing, m ; ft 

T absolute temperature, K; ®R 

t trailing-edge thickness, m; ft 

u distance in tangential direction, m; ft 

V fluid velocity, m/sec; ft/sec 

X distance along blade surface from forward stagnation point, 
m; ft 

Y distance from surface normal to boundary layer, m; ft 

y distance from surface normal to boundary layer expressed as 
fraction of boundary-layer thickness 
a fluid flow angle from axial direction, deg 
a« blade stagger angle from axial direction, deg 
y ratio of specific heat at constant pressure to specific heat at 

constant volume 

5 boundary-layer displacement thickness, m ; ft 

boundary-layer thickness, m; ft 
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6 boundary-layer momentum thickness, m ; ft 

M viscosity, (N) (sec)/m2; lb/ (ft) (sec) 

p density, kg/m’; Ib/ft’ 

^ boundary-layer energy thickness, m ; ft 

Subscripts: 

cr critical flow conditions 

eff effective 

/ friction 

fs free stream 

h hub 

i ideal 

inc incompressible 

m mean section 

min minimum 

mix mixing 

opt optimum 

p pressure surface 

ref reference 

s suction surface 

t tip 

te trailing edge 

tot total 

X axial component 

0 blade-row inlet 

1 just beyond trailing edge of blade row 
la just within trailing edge of blade row 

2 downstream uniform state 

3D three dimensional 

Superscripts: 

' absolute total state 

* dimensionless value 
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CHAPTER g 


Miscellaneous Losses 

By Richard J. Roelke 


In the last chapter, the boundary-layer losses associated with the flow 
process in the blade channel were discussed. To determine the overall 
design-point efficiency of a turbine, other losses must also be considered; 
these include tip-clearance loss and disk-friction loss. In some instances, 
these losses represent a very small part of the turbine output and may be 
neglected; however, in other instances, these losses can be of such mag- 
nitude as to influence the selection of the turbine design point. The sum 
of these losses normally comprises all the losses that are considered in the 
design of a full-admission axial-flow turbine. If, however, a partial- 
admission turbine is being considered, there are additional losses that 
must be included. The partial-admission losses usually considered are the 
pumping loss in the inactive blade channels and the filling-and-emptying 
loss in the blade passages as they pass through the admission arc. Finally, 
a loss that occurs at off-design operation of any turbine is the incidence 
loss, which will also be covered herein. 


TIP-CLEARANCE LOSS 

Because a turbine must operate with some clearance between the tips 
of the rotor blades and the casing, some fraction of the fluid leaks across 
the tips, thus causing a reduction in turbine work output. This leakage 
loss is affected, first of all, by the nature of the tip geometry; that is, by 
the amount of radial clearance, by recesses in the casing, and by tip 
shrouds. For a given tip geometry, the amount of blade reaction affects 
the leakage loss, since a large pressure difference across the tip (high 
reaction) causes more higher-kinetic-energy flow to leak through the tip 
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gap from the pressure side to the suction side of the blade. With an un- 
shrouded blade, this leakage flow not only causes a loss due to its own 
reduced work, but also causes an unloading of the blade, primarily in the 
tip region. Analytical evaluation of the drop in turbine efficiency caused 
by tip-clearance leakage is inherently difficult because of the complex flow 
problem. Several empirical expressions for clearance loss have been 
developed, and some of these are summarized in reference 1; however, 
they are rather complicated, and the author states that none is entirely 
satisfactory, 

A number of tests have been made at the NASA Lewis Research Center 
to determine the effect of tip clearance and tip geometry on axiaLflow 
impulse and reaction turbines. An examination of some of the results of 
these tests helps to obtain a better understanding of the tip-clearance 
loss. Figure 8-1 shows the angle traces at the blade exit of a 5-inch single- 
stage turbine (ref. 2). Two things to be noted from the angle traces are 
that the flow in the clearance space and near the tip was not fully turned, 



Figure 8^1 . Variation of exit flow angle with radius ratio for four rotor tip clearances. 

(Data from ref. 2,) 


126 



MISCELLANEOUS LOSSES 



even at the smallest clearance tested, and that underturning of the flow 
increased with increasing tip clearance, and this effect occurred all the 
way down to the hub. This undertuming of the flow unloads the blade 
aerodynamically and results in lower turbine output and efficiency. The 
decrease in efficiency for this turbine, as well as for two others, is shown 
in figure 8-2. 

The solid lines in figure 8-2 represent test results from single-stage 
(ref. 2) and two-stage (ref. 3) reaction turbines and from a single-stage 
impulse turbine (ref. 4). All turbines were unshrouded. The importance 
that the level of reaction plays in the clearance loss is clearly evident from 
the figure. For the same ratio of tip clearance to blade height, the losses in 
efficiency for the reaction turbines were about double that for the impulse 
turbine. 

The dashed lines in figure 8-2 are estimates of the efficiency losses for 
the two single-stage turbines (refs. 2 and 4) as obtained from the curves 
published in reference 5 (as fig. 1.6) and reproduced here as figure 8-3. 
Extrapolation of the experimental data of figure 8-2 shows that figure 8-3 
gives satisfactory estimates of tip-leakage loss for small tip clearances. 
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Figure 8-3. — Tip-clearance correlation for unshrouded blades. (Data from ref. 5.) 


Reviewing the results shown in figures 8-2 and 8-3, it is apparent that the 
loss in efficiency increases with increasing reaction, and for moderate to 
large ratios of tip clearance to blade height, the loss is appreciable. 

In addition to reducing the tip clearance, methods for reducing the tip- 
leakage losses include recessing the casing above the blade tip while in- 
creasing the blade height, and adding a tip shroud. These loss-reduction 
schemes can be used either individually or in combination. The single- 
stage impulse turbine of reference 4 was tested at several ratios of tip 
clearance to blade height, both without and with the recessed casing and 
the tip shroud. Figure 8-4 shows the three general configurations tested 
in reference 4, and the turbine-performance results are shown in figure 8-5. 
A clearer understanding of the performance characteristics is possible if 
the loss mechanisms are considered. The factors affecting turbine wwk 
for the reduced blade-height configuration as compared to a zero-clearance 
configuration consist of (1) reduced blade loading area, (2) clearance-gap 
leakage flow, (3) mixing of the leakage flow with channel throughflow, 
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Flow 


Tip shroud-\ 



(c) 


(a) Reduced blade height (relative to (b) Recessed casing, 

zero-clearance blade height). 

(c) Shrouded rotor. 

Figure 8—4. — Tip-clearance configurations investigated for impulse turbine (ref. 4). 


and (4) blade unloading (as a result of flow going from the pressure side 
to the suction side). With the recessed-casing configuration, the blade 
extended to the passage outer radius and was of constant height as the 
clearance gap was changed by varying the amount of casing recess. There- 
fore, the reduced blade loading area was eliminated, and the leakage flow 
was reduced because of the indirect leakage path. With the shroud added 
to the blade, the blade unloading was eliminated, and the leakage flow was 
further reduced. Note from figure 8—5, however, that at tip-clearance 
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ref. 4.) 


fractions below some value, about 0.035 in this instance, the shroud no 
longer provides an increase in efficiency. This can be attributed to an 
increasing friction loss between shroud and casing as the clearance gap 
is decreased. 

The comparative results shown in figure 8-5 for different blade-tip 
geometries are dependent upon that particular design and may not apply 
to other turbines. This is particularly true of the shrouded blade, since 
the leakage flow depends not only on the clearance span and pressure 
difference but also on the number of seals used. With respect to the 
recessed-casing configuration, it should be noted that the blade should not 
extend into the recess. If it does, the overlapping section will just be 
churning stagnant fluid and creating additional losses. 

In summary, tip-clearance loss presents a complicated flow problem 
influenced by many factors and is not easily predicted with consistent 
accuracy. The clearance gap required for a turbine depends primarily on 
diameter (larger clearance for larger diameter) and, as seen previously, 
the loss increases as the ratio of clearance gap to blade height increases. 
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For any given diameter, therefore, the tip-clearance loss increases with 
increasing hub- to tip-radius ratio. It becomes increasingly difficult to 
maintain a desired small ratio of clearance gap to blade height as the 
turbine, and hence the blade height, becomes smaller. For a given radius 
ratio, therefore, the loss is more severe for small turbines and less severe 
for larger turbines. If tip leakage is considered to be a problem in a 
particular case, it might be worthwhile to carry out tests to evaluate the 
leakage effects. 

DISK-FRICTION LOSS 

The disk-friction loss (or windage loss) is due to the skin friction and 
circulation of fluid between the rotating disk and the stationary casing. In 
addition, some turbines for hot applications, for example aircraft engines, 
have a small steady stream of lower-temperature gas that bathes and cools 
the rotor disk. This cooling gas flows along the rotor-disk surface from 
near the engine centerline outward to the base of the blades. The qualita- 
tive nature of the flow patterns around rotor disks without and with 
throughflow of cooling gas are shown in figure 8-6. Equations for es- 
timating the associated losses are presented herein. 

No Throughflow 

For the case with no throughflow, as in figure 8-6 (a) , the thin layer of 
fluid close to the rotating surface is thrown outward by centrifugal action 



(a) Without throughflow. (b) With throughflow. 

Figure 8-6.— Flow patterns for rotating disks. 
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and returns via the stationary wall to the inner radius, thereby building 
up a continuous circulatory effect. Consider an element of area on one 
side of the disk 

dA—2Tcrdr (8-1) 

where A is the area, in m* or ft^, of one side of the disk, and r is the radius, 
in m or ft, of the area element dA. The fluid shear stress r, in N/m^ or 
Ib/ft^, acting over this area at the radius r produces a resisting torque to 
the disk rotation of 


UiYl 0 

—^~T2wr^dr ( 8 - 2 ) 

where Mo is the resisting torque, in N-m or Ib-ft, for both sides of a disk 
in the case of no throughflow. The shear stress can be expressed as 

T=^pVu^ (8-3) 

where 

Cf fluid shear-stress coefficient 

g conversion constant, 1 ; 32.17 (Ibm) (ft)/(lbf) (sec^) 

p density, kg/m^; Ib/ft’* 

Vu tangential component of fluid absolute velocity, m/sec; ft/sec 

At the disk surface, the fluid tangential velocity is 


= roj 


(8^) 


where w is the angular velocity, in rad/sec. By substituting equation (8-4) 
into (8-3), the total torque for both sides of the disk can be written as 

O-TT 

Mo= —Cfpe^^r^dr (8-5) 

•'o 9 

where a is the disk rim radius, in m or ft. Performing the int(^gration yields 


Mo — Ca/.oP — (8—6) 

2(7 

where Cm,o is a tonpie coefficient for the case of no throughflow. The disk- 
friction loss expressed as power is th(Mi th(‘ torqu(‘ tim(‘s the angular 
velocity: 


MoO} pcj^a-^ 


(^7) 


where Pd/ is the disk-friction power loss, in W or Btu/sec, and J is a 
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conversion constant (equal to 1, or 778 (ft) (lb)/Btu). The form of 
equation (8-7) that is found in most handbooks is 


where 

Pdf=^KdfpN^Dr^ 

(8^8) 


disk-friction power-loss coefficient 


N 

rotative speed 


Dr 

disk rim diameter 



A number of investigators have published values of the constant Kd/ in 
equation (8-8) to be used for different circumstances, while others have 
made small changes to the exponents to better fit the available data. The 
wide assortment of semiempirical equations used to predict this loss is, no 
doubt, due to variations of the test-apparatus geometry, the somewhat 
oversimplified model from which equation (8-7) is derived, and the 
existence of different types of flow that can occur in the space between the 
rotor and the casing. One thing that can be noted from equation (8-7) 
or (8-8) is that for a given blade speed, lower loss is obtained by having 
a smaller diameter and a higher rotative speed. 

An extensive investigation has been conducted (refs. 6 and 7) to deter- 
mine the effect of chamber proportions on disk friction and to present a 
clearer picture of the several modes of flow that may exist. In general, 
four modes of flow, or flow regimes, can exist in the axial space between the 
casing and the rotating disk, depending on the chamber dimensions and 
the flow Reynolds number. The torque coefficient Cm,o was evaluated 
both theoretically and experimentally in each regime. A description of 
each regime and the associated equations for the torque coefficient are as 
follows: 

Regime I : Laminar Flow, Small Clearance. Boundary layers on the 
rotor disk and casing are merged, so that a continuous variation in 
velocity exists across the axial gap s. Figure 8-7 (a) indicates the nature 
of the variations in the radial and tangential components of fluid velocity 
at a given radius in the gap. The best equation for torque coefficient, both 
theoretically and empirically, is 


(s/a)R 


(8-9) 


where s is the axial distance, in m or ft, between disk and casing, and R is 
the Reynolds number defined as 


R^ 






( 8 - 10 ) 


where p is the dynamic viscosity, in (N) (sec)/m^ or lb/ (ft) (sec). 


133 



TURBINE DESIGN AND APPLICATION 


Radial Tangential 

component component 

of velocity of velocity 



(a) 


Radial Tangential 

component component 

of velocity of velocity 



(a) Flow regimes I and III. 

(b) Flow regimes II and IV. 

Figure 8-7— Velocity patterns around rotating disks without throughflow. 


Regime II. Laminar Flow, Large Clearance. The combined thick- 
ness of the boundary layers on the rotor and on the casing is less than the 
axial gap, and between these boundary layers there exists a core of rotat- 
ing fluid in which no change in velocity occurs. Figure 8-7 (b) shows 
the variations in the radial and tangential velocity components for this 
case. The best theoretical and empirical equations for torque coefficient 
are 


C 




( 8 - 11 ) 


where Cn is a function of (s/a), as shown in figure 8-8(a), and 
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(a) Flow regime II. 

(b) Flow regime IV. 

Figure 8-8. — Evaluation of torque coefficients. (Data from ref. 6.) 


CAf — 


3.70 


( 8 - 12 ) 


respectively. 

Regime III: Turbulent Flow, Small Clearance. The turbulent 
counterpart of Regime I. The best theoretical and empirical equations 
for torque coefficient are 


0.0622 


( 8 - 13 ) 


and 


C 


M.a — 


0.080 


( 8 - 14 ) 


respectively. 
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Regime IV; Turbulent Flow, Large Clearance 
counterpart of Regime II. The best theoretical and 
tions for torque coefficient are 


The turbulent 
empirical equa- 


Cm,o— 


Civ 


Jivh (8--15) 

where Civ is a function of {s/a)^ as shown in figure 8-8 (b), and 


Cm.o = 


0.102 (5/a) 


(8~16) 


respectively. 

The particular flow regime that exists at any Reynolds number can be 
determined by plotting torque coefficient against Reynolds number from 
equations (8-9), (8-11), (8-13), and (8-15), as shown in figure 8-9 for 
several values of s/a. The discontinuities (changes in slope) in the lines 
01 figure 8-9 indicate transition from one regime to another. In this figure 
the flow regimes are determined by matching the slopes of the lines with 



Figure 8 9. Theoretical variation of torque coefficient with Reynolds number for no 
throughfiow. (Data from ref. 6.) 
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those shown in the insert in the figure. Torque-coefficient values (ref. 6) 
determined experimentally with a 50.8-centimeter (20-in.) disk rotated in 
water and oil for several values of s/a verify the theory. 

Throughflow 

For the case of the rotating disk with throughflow, as in figure 8-6(b), 
the friction torque increases with the throughflow. This problem has been 
analyzed for low values of throughflow with regime-IV flow. In this case, 
it is assumed that the fluid enters the chamber near the centerline with no 
angular velocity and leaves at the rim with some angular velocity K„om. 
The symbol K„ represents the ratio of the angular velocity of the rotating 
core of gas to the angular velocity of the disk. The increase in torque, 
AAf, over that without throughflow is the rate of change of angular 
momentum of the fluid flowing through the system: 

Q 0 

AAf = 2p — (/Co0>a)a = 2p — (8-17) 

9 9 

where Q is the volumetric throughflow rate, in mVsec or ft’/sec, in the 
clearance space on one side of the disk. The total torque for the through- 
flow case is then 


M = Mo+AM^^~^+2~QKmm^ (8-18) 

^9 9 

The value of Ko is approximately 0.45 for s/a ratios from 0.025 to 0.12. 

An assessment of the power loss can be obtained by calculating the 
friction torque of the throughflow case compared to that of the no- 
throughflow case: 


M ^ ^ 2pQK^^ _ 1 I 
Afo 

- Cm,(.p<oW 


(8-19) 


Substituting equation (8-16) for Cm.® yields 

Q 


M _ 

0.0255 (s/a) 


= 1+39.2 


Ko 


(s/a)*/'® 

where T is a dimensionless throughflow number defined as 


Q 

T = — 
oxi’ 


(8-20) 


( 8 - 21 ) 


According to the data of reference 7, equation (8-20) predicts values 
that are somewhat high; moreover, the effect of s/a is not accurately 
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Ratio of axial gap 
to disk rim radius, 



Figure 8-10. — Empirical variation of torque with throughflow number. (Data from 

ref. 7.) 


given by Empirically, the test data are represented to within 

±5 percent by the relation 


M 

M, 


l + lS.9Ko 


T 


( 8 - 22 ) 


Equation (8-22) is plotted in figure 8-10 for several s/a values. 


PARTIAL-ADMISSION LOSSES 

Full-admission axial-flow turbines are used for most applications; 
however, unusual conditions sometimes arise for which a partial-admission 
turbine may be a better choice. If, for example, the design mass-flow^ rate is 
so small that a normal full-admission design would give very-small blade 
heights, then it may be advantageous to use partial admission. The losses 
due to partial admission with long blades may be less than the leakage 
and low Reynolds-number losses of the full-admission turbine having short 
blades. In addition, for a given rotative speed, partial admission allows the 
freedom of larger diameter and higher blade-jet speed ratios. Also, the 
use of partial admission may be a convenient way to reduce power output 


138 



MISCELLANEOUS LOSSES 


of an existing full-admission turbine (physically block some of the stator 
passages). In general, partial-admission turbines have high specific-work 
output and low volumetric-flow rates. 

As mentioned previously in this chapter, the partial-admission losses 
are the pumping loss in the inactive blade channels and the filling-and- 
emptying loss encountered as the blades pass through the active sector. 
This latter loss has been referred to as expansion, scavenging, or sector 
loss. The mechanisms of partial-admission losses are not clearly or fully 
understood, but they do result in a decrease in output power and efficiency 
when compared to the same turbine operating with full admission. 

The pumping loss is that loss caused by the inactive blades rotating in 
a fluid-filled casing, and expressions for it are somewhat similar in form to, 
and often combined with, the expression for the disk-friction loss. These 
expressions all seem to trace back to reference 8, where the results of 
several experimental investigations are summarized. The equations for 
estimating pumping-power loss that resulted from these investigations 
showed that the effects of blade height and diameter on the pumping- 
power loss are quite uncertain, as evidenced by variations in the exponents 
on these terms. Further, the nature and location of obstructions (adjacent 
blade rows, casing wall, etc.) or lack of such in the vicinity of the three 
open sides of the blade channel were accounted for only by differences in 
the empirical loss coefficient. Therefore, it appears that a generally 
applicable expression for pumping-power loss is yet to be found. 

The one equation perhaps most often used is 

P, = K^p W (1-6) (8^23 ) 

where 

Pp pumping-power loss, W ; (ft) (lb)/sec 

Kp pumping-power loss coefficient, 1/m*^^; (Ibf) (sec2)/(lbm) (ft*^*^ 

Um blade mean-section speed, m/sec; ft/sec 
I blade height, m; ft 

Dm blade mean-section diameter, m; ft 

€ active fraction of stator-exit area 

The value of the coefficient Kp as reported in reference 8 and converted 
to the units used herein is 3.63 1/m*^^, or 0.0105 (Ibf) (sec*)/(lbm) (ft*^*), 
for an unenclosed rotor. For the same rotors enclosed, the coefficient 
values were one-quarter to one-half of the above values. More recently, 
the combined disk-friction and pumping losses for a single-stage rotor 
enclosed by the turbine housing were reported in reference 9. If a disk- 
friction loss estimated by equation (8~7) is subtracted from the combined 
losses of reference 9 and if the remaining loss is converted to the form of 
equation (8-23), the coefficient Kp is found to be 5.92 1/m*^^ or 0.0171 
(Ibf) (sec*)/ (Ibm) (ft*^*). This is significantly higher than the coefficients 
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reported in reference 8, and the difference is attributed to primarily the 
lack of an adequate loss model. 

The other partial-admission loss shall herein be called the sector loss. 
Imagine a blade channel as it just starts to enter the active sector. It is 
filled with relatively stagnant fluid that must be pushed out by the high- 
momentum fluid leaving the nozzle. This scavenging will continue until 
the blade channel is completely within the active sector. As the blade 
channel passes out of the active sector, a second sector loss occurs. As the 
inlet to the blade channel is cut off from the nozzle active arc, less and 
less high-momentum fluid enters the channel. Since this fluid has the 
entire blade channel area to flow into, it is rapidly diffused as it flows 
through the rotor. These losses cause an overall decrease in the momentum 
of the fluid passing through the rotor, thus decreasing the available energy 
of the fluid. It was reported in reference 10 that this decrease in momen- 
tum may be found by multiplying the rotor-exit momentum by a loss 
coefficient 

K..(l-|) (S^-24) 

where p is the rotor-blade pitch, in m or ft, and f is the nozzle active arc 
length, in m or ft. Effectively, K, is a rotor velocity coefficient that ac- 
counts for the sector loss. 

The effect of the sector loss on turbine efficiency is determined as 
follows. With the use of equations (2-6) and (2-14), from volume 1, and 
the associated velocity diagram geometry, we can express the specific 
work of an axial-flow turbine as 


= % {Wu,i-WuA)=—r (WismPi-WiSmfit) (8-25) 

gJ gJ 


where 

AA' turbine specific work, J/kg; Btu/lb 

Wu tangential component of relative velocity, m/ sec ; ft/ sec 

W relative velocity, m/ sec ; ft/ sec 

P fluid relative angle measured from axial direction, deg 

The subscripts 1 and 2 refer to the rotor inlet and exit, respectively. For 
an impulse turbine (which most partial-admission turbines are), where 


A/i' = % IT, sin 0i(l+K^) (8-26) 

gJ 

where is the rotor relative-velocity ratio Wi/W\ for the full-admission 
turbine. For the partial-admission turbine, applying the sector loss 
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coefficient yields 


W2 = K^K.Wx 


So, for the partial-admission turbine, 


sin fi^{\+K^K.) 


Since efficiency is 



(8^27) 


(8~28) 


(S-29) 


whore Lhid is the turbine ideal specific work, in J/kg or Btu/lb, the ef- 
ficiency of the partial-admission turbine with respect to that of the full- 
admission turbine is 


Vya _ 

T) Ah' 


(8^30) 


Substituting equations (8-26) and (8-28) into equation (8-30) then 
yields 


Vpa = v 


l+K^K, 


(8^31) 


The efficiency penalty expressed by equation (8-31) accounts for the 
sector loss only ; the pumping loss discussed earlier will reduce the overall 
efficiency further. Equation (8-24) indicates that a partial-admission 
turbine rotor should have closely spaced blades to reduce the sector loss ; 
however, as more blades are added to the rotor, the blade profile loss will 
increase. Also, the effect of the number of rotor blades on the pumping 
loss is not known. Therefore, the complete optimization of a partial- 
admission design cannot be done analytically at present. 

In the study of reference 9, the efficiency of a small axial-flow turbine 
was determined over a range of admissions from 12 to 100 percent. The 
total loss due to partial-admission operation was taken as the difference 
between the full- and the partial-admission efficiencies. The blade pumping 
and disk-friction losses were measured separately and were subtracted 
from the total partial-admission loss to give what was called other partial- 
admission losses. These other losses include the sector loss and any loss 
due to leakage from the active sector to the inactive sector. The partial- 
admission losses of reference 9 are plotted against admission-arc fraction 
in figure 8-11. The combined pumping and disk-friction loss increased 
with decreasing arc fraction, while the other losses remained nearly con- 
stant over the range of arcs tested. 

Predicted efficiencies (from ref. 10) are plotted against blade-jet speed 
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Fiquk£ 8r-ll. — Variation of partiaLadmission losses with active fraction of stator 

area. (Data from ref. 9.) 



Figure 8-12. — Design-point performance of partial- and full-admission turbines. 

(Data from ref. 10.) 


ratio (see discussion in vol. 1, ch. 2) in figure 8-12 for a particular turbine 
operating with full admission and with three different amounts of partial 
admission. The expected reduction in peak efficiency with reduced arc of 
admission is seen. The important thing to note from this figure is the 
reduction in optimum blade-jet speed ratio as admission arc is reduced. 
Aerodynamic efficiency is a maximum at a blade-jet speed ratio of 0.5, 
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irrespective of admission arc, and decreases with decreasing blade speed. 
Blade-pumping and disk-friction losses, which decrease with decreasing 
blade speed, become a larger part of the gross aerodynamic power as 
admission arc decreases. Therefore, as admission arc is reduced, the 
maximum net output power (aerodynamic power minus blade-pumping 
and disk-friction power) is obtained at lower blade speeds. Thus, for the 
design of a partial-admission turbine, the partial-admission losses must be 
factored into the design before an optimum or near-optimum blade-jet 
speed ratio can be selected. 


INCIDENCE LOSS 

The incidence loss is that loss which occurs when the gas enters a blade 
row (either stator or rotor) at some angle other than the optimum flow 
angle. Flow incidence would normally only occur at ofT-design conditions, 
since, theoretically at least, all gas and blade angles are matched at the 
design condition. The nomenclature used when speaking of incidence is 
shown in figure 8-13. The dashed line running through the blade profile 
is the camber line and defines the blade inlet angle. The incidence angle 
is defined as 

i = a — oh (8-32) 

where 

i incidence angle, deg 
a fluid flow angle from axial direction, deg 
«6 blade inlet angle from axial direction, deg 

The fluid flow angle must be the absolute angle for stators and the relative 
angle for rotors. The incidence angle may be positive or negative, as 
indicated in figure 8-13. The sign of the incidence angle is important 
because cascade tests have shown that the variation of loss with incidence 
angle is different for positive and negative angles. 



Figure 8-13. — Blade incidence nomenclature. 
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Figurk 8 - 14 . — Characteristics of l)lade incidence loss. 



The general nature of the variation of incidence loss with incidence angle 
is shown by figure 8-14, which represents a summary of cascade test 
results. The loss curve is not symmetrical about the zero incidence angle, 
but shows a loss that is larger for positive incidence than for negative 
incidence. This may be due to some local separation on the suction surface 
at large positive incidence, as indicated in figure 8-15, and the lack, or 
smaller area, of separation at the same value of negative incidence. Also, 
blades in which the mean acceleration of the gas flow is large (high- 
reaction blades) have a wide range of incidence over which loss is low, 
whereas low-reaction blades have higher losses for the same incidence 
range. 

Another thing to be noted from figure 8-14 is that the minimum loss 
does not occur at zero incidence, but at some small amount of negative 
incidence. This may be explained by the sketch of figure 8-16. The stag- 
nation streamlines for two inlet flow angles are shown; one at zero in- 
cidence and the other at some small negative incidence with respect to 
the blade inlet angle. Both tests and theory show that the stagnation 
streamline curves upward as the flow impacts on the blade leading edge, 
and the true zero incidence occurs when there is some negative incidence 
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Figure 8 - 16 . — Curvature of stagnation streamline at blade inlet. 


relative to the free-stream flow. The incidence angle at minimum loss is 
usually —4® to -“8°. Because of this, some turbine designers design their 
blades with a small amount of negative incidence, while others do not 
because of the small difference involved. 

The magnitude of the incidence loss takes on importance when the 
off-design performance of a turbine must be predicted. A method for 
determining incidence loss based on test data is described in reference 11. 
An analytical method is described here with the aid of figure 8-13. The 
inlet velocity Vi can be resolved into a component T„ normal to, and a 
component Vp parallel to the blade inlet direction (camber line at inlet). 
If it is assumed that the parallel component passes through the blade row 
without any entry loss and that the normal component is entirely lost, 
the recovered kinetic energy is 


2gJ 2gJ\Vj 2gj'"''^ " 

and the kinetic-energy loss due to incidence is 

Li = — (l-cos^t) 


(8-33) 


(a-34) 


In order to account for the differences in loss variation with positive 
and negative incidence, the effect of blade-row reaction, and the minimum 
loss not occurring at zero incidence, equation (8-34) has been generalized 
to 



[1 — cos’' 


{i '^opi)~\ 


(8-35) 


where iopi is the optimum (minimum-loss) incidence angle. This type of 
equation has proved satisfactory when used in off-design performance 
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prediction methods such as that of reference 12. Where specific incidence- 
loss data are lacking, values of n = 2 for negative incidence and n = 3 for 
positive incidence have been used satisfactorily, 
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SYMBOLS 

A area on one side of rotor disk, m*; ft* 

a disk rim radius, m ; ft 

C II coefficient used to evaluate C.v.o in regime II by equation (8-11) 

C IV coefficient used to evaluate in regime IV by equation (8-15) 

Cf fluid shear-stress coefficient 

Cm,o torque coefficient with no throughflow 

D diameter, m; ft 

/ nozzle active arc length, m; ft 

g conversion constant, 1 ; 32.17 (Ibm) (ft)/(lbf) (sec*) 

A/i' turbine specific work, J/kg; Btu/Ib 

A/i,<i turbine ideal specific work based on ratio of inlet-total pressure 
to exit-static pressure, J/kg; Btu/lb 
i incidence angle, deg 

J conversion constant, 1 ; 778 (ft) (lb) /Btu 

Kdf disk-friction power-loss coefficient 

Ko ratio of rotating-core angular velocity to disk angular velocity 
Kj, pumping power loss coefficient, l/m^^*; (Ibf) (sec*)/(lbm) (ft=*^*) 
Kt sector loss coefficient 

rotor velocity coefficient for full-admission impulse turbine 
Li incidence loss, J/kg; Btu/lb 

I blade height, m; ft 

M frictional resistance torque for both sides of rotor disk, N-m; 
Ib-ft 

N rotative speed, rad/sec; rev/min 

n exponent in equation (8-35) 

Pdf disk-friction power loss, W ; Btu/sec 

Pp pumping power loss, W ; Btu/sec 

;> rotor-blade pitch, m ; ft 

Q volumetric throughflow rate, mVscc ; ftVsec 

R Reynolds number 

r radius, m; ft 

5 axial distance between rotor disk and casing, m; ft 

U blade speed, m/sec; ft/sec 

V absolute velocity, m/sec; ft/sec 

W relative velocity, m/sec; ft/sec 

a fluid flow angle from axial direction, deg 

ab blade inlet angle from axial direction, deg 

fluid relative angle measured from axial direction, deg 
€ active fraction of stator exit area 

7} turbine static efficiency 

fji dynamic viscosity, (N) (sec) /m*; lb/(ft) (sec) 

p density, kg/m*; Ib/ft^ 
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T fluid shear stress, N/m^; Ib/ft^ 

T throughflow number defined by equation (8-21) 

CO angular velocity, rad/sec 

Subscripts : 
m mean section 

n component normal to blade inlet direction 

o no throughflow 

opt optimum 

p component parallel to blade inlet direction 

pa partial admission 

r disk rim 

u tangential component 

^ f rotor inlet 

\blade-row inlet 
2 rotor exit 
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